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"A large cloud chamber containing eight lead plates was operated without counter control 
at 10,000 feet. In 8500 photographs, 1090 cascade showers and 58 nuclear disintegrations were 
observed. The differential energy spectrum for the shower-producing rays, which was obtained 
from the relative frequency of occurrence for showers of a given size, is proportional to E+” 
for E=2X10* ev and to E-*“ for E = 10° ev. The shower-producing rays with E>2 X 108 ev con- 
stituted 6.5 percent of the radiation observed in the cloud chamber; one-third of the shower- 
producing rays were photons. Approximately one-third of the observed nuclear disintegrations 
contained particles that penetrated at least 0.7 cm of lead and four contained particles that 
penetrated at least 2.8 cm of lead. Most of the disintegration particles were protons or meso- 
trons. The initiating particles were neutrons and protons, with the former predominant for the 
lower energy disintegrations. There appears to be no sharp line of distinction between low 
energy and high energy disintegrations. No correlation was observed between the disintegrations 


and the cascade showers. 





ARIOUS types of cloud-chamber experi- 
ments have been performed in order to 
analyze the cosmic radiation at altitudes of ten 
to fifteen thousand feet.' In the present experi- 
ment the details of the interaction of cosmic rays 
with lead were studied at an altitude of 10,000 ft. 
A large cylindrical cloud chamber (30 by 30 cm) 
was operated at Tioga Pass, Yosemite National 
Park, without counter control. The chamber con- 
tained eight lead plates each consisting of a 
0.6-cm lead plate and two reflecting steel sheets 
plated with chromium. The equivalent lead 
thickness, calculated on the basis of densities, 
was 0.7 cm, and since the steel constituted such 
a small fraction of the absorber, it is assumed in 
all subsequent calculations that the plates were 
simply 0.7-cm lead sheets. Since there was suf- 


1 Several of these experiments will be mentioned later in 
the paper. 
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ficient illumination to obtain images of individual 
drops, diffuse tracks could be photographed, and 
these are included in the statistics in some cases. 
The chamber was filled with air and the saturated 
vapor from a 1 : 3 liquid mixture of water and 
ethyl alcohol to a total pressure of ~70 cm when 
compressed and ~60 cm when expanded. Other 
details of the apparatus have been given pre- 
viously.” 

The lead plates made possible the classification 
of the ionizing particles in two groups: particles 
of electronic mass, and particles of many times 
electronic mass. (mesotrons and protons). Fur- 
thermore, the details of the interaction of each 
type of particle with lead could be observed. 
Essentially all of the photons initiated cascade 
showers in the lead and thus could be detected, 


?W. E. Hazen, Phys. Rev. 64, 7 (1943). 
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Fic. 1. Three time-associated showers were photographed 
in this expansion. The initiating particles are probably part 
of an air shower. There are 30 particles in the largest shower 
at its maximum, which occurs under the sixth lead plate (4.2 
cm of lead). The energy of the particle that initiated this 
shower was 1.5 10° ev and the predicted position of the 
maximum is 4.5 cm of lead. Reflected images of the tracks 
can be seen in the plate surfaces, most clearly in the top and 
bottom plates. 


whereas only a small fraction of the non-ionizing 
particles interacted with the lead. The analysis 
of the data pertaining to the production of col- 
lision electrons by the heavier particles has 
already been given.’ 


CASCADE SHOWERS 


If an electron or photon has an energy greater 
than 10 ev, it is very unlikely to pass through a 
series of 0.7-cm lead plates without producing an 
observable cascade shower. In the cases where 
the initiation of the shower is observed, it is 
possible to distinguish between electrons and 


.than that for the larger showers, which nearly 
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photons as the originating rays; where the 
maximum of the shower is observed, the energy 
of the originating ray can be estimated. 

In the 8500 useful photographs, 1090 cascade 
showers with maxima of four or more particles 
were observed ; 731 were classified as giving rige 
to sharp tracks, the rest as giving rise to diffuse 
tracks. Examples of the photographs are repro. 
duced in Fig. 1. In the same series of photo. 
graphs, 8678 sharp tracks from penetrating 
particles were observed. A size-frequency block 
diagram for the showers is given in Fig. 2. The 
relative number of four-particle showers is too 
low because many showers of this size originated 
in one plate and stopped in the next plate. This 
situation makes the chance of observation less 


















always passed through several plates. 

The energy values shown along the scale of 
abscissae were obtained from the number of 
particles at the shower maximum by the use of 
Eq. (2.104) of Rossi and Greisen.* The number 
of particles at the shower maximum is not appre. 
ciably different for showers initiated by electrons 
and those initiated by photons of the same 
energy. When plotted with logarithmic scales 
the data are only reasonably well represented by 
a straight line (wth a slope of —2.4) ; see Fig. 3. 
If a simple power law is used to describe portions 
of the differential energy spectrum, the exponent 
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Fic. 2. Frequency of occurrence of cascade showers it 
lead as a function of the size of the shower. The size of the 
shower is expressed in terms of the number of particles at 
the maximum of the shower. The energy scale is obtained 
from shower theory [Eq. (2.104) of Rossi and Greisen, 
reference 3]. 








2B. Rossi and K. Greisen, Rev. Mod. Phys. 13, 2@ 
(1941). 
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varies from —1.7 for the lowest energies (2  10* 
ev) to —3.0 for the highest energies (10° ev). 

The intensity of the shower-producing par- 
ticles relative to the penetrating particles was 
obtained by considering sharp tracks only (the 
uncertainty in the equivalence of the total times 
of observation is thereby reduced). Since the 
total number of sharp tracks caused by showers 
with maxima greater than four particles was 605, 
while in the same series of photographs there were 
8678 sharp penetrating tracks, the relative 
number of high energy electrons and photons 
(E>2X 10%) is seven percent at 10,000 ft. From 
the cases where it was believed that the beginning 
of the shower was observed, the ratio of photons 
to electrons was found to be 1 : 2.5. 

Auger has compared the absorption of cosmic 
rays in lead and in aluminum, at about the same 
altitude as in the present experiment. Since 
Auger used a counter telescope with the absorber 
placed between the counters, the number of 
photons that would have been recorded is neg- 
ligible. He obtained a ratio of 8 : 54 or 15 percent 
for the number of shower-producing particles 
with E>10* ev compared with the number of 
penetrating particles. If the spectrum from the 
present data is extrapolated to an energy of 
10* ev, the value for the relative number of 
shower-producing rays (photons and electrons) 
with E>10*® would not be much greater than 15 
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Fic. 3. The energy spectrum for photons and electrons 
(shower-producing rays). The ordinates represent the 
number of rays in an energy band of width 45 Mev. The 
experimental point with the smallest abscissa, which is 
represented by +, is relatively too low because showers of 
this size seldom penetrate more than one lead plate and 
therefore are less likely to be observed than larger showers; 
the other experimental observations are represented by 
lines whose lengths indicate the probable errors. 





*P. V. Auger, Phys. Rev. 61, 684 (1942). 
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percent, even if the largest reasonable slope for 
the lower end of the spectrum is used in the 
extrapolation. A comparison of the results from 
the two experiments would thus lead one to the 
conclusion that the number of photons is neg- 
ligible; however, the experimental uncertainties 
are fairly large in both cases and the result of 
the comparison cannot be considered as in 
conflict with the above-mentioned value of 
1 : 2.5 for the ratio of photons to electrons. Most 
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Fic. 4. The observed number of star particles as a func- 
tion of their range in lead. The auxiliary scales of abscissae 
give the energies for protons and mesotrons. The dotted 
lines indicate the results of the correction for the finite size 
of the lead plates. 


of the shower-producing rays originate as elec- 
trons,> and hence one would expect to observe 
more electrons than photons in a region where the 
rays are still being produced. 


STAR SHOWERS OR NUCLEAR 
DISINTEGRATION 





A few stars have been observed in the course 
of many cloud-chamber experiments and in one 
recent investigation 156 stars were observed.® In 
the present experiments a total of 58 stars con- 
sisting of two or more particles has been found 
in the photographs. The large area and close 
spacing of the lead plates resulted in a reasonably 
large chance that a particle would find several 
lead plates in its path if its range was several 
centimeters of lead. The interaction of a particle 
with the lead plates provides information con- 
cerning the nature of the particle. 

Thirty-one of the stars consisted entirely of 


5K, Greisen, Phys. Rev. 63, 323 (1943). 
SW. M. Powell, Phys. Rev. 61, 670 (1942). 
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the particles with range less than one lead plate 
(0.7 cm) and greater than 2.5 cm of air. Included 
in these thirty-one stars are the two which 
originated in the gas rather than in one of the 
lead plates. The rest of the stars contained 
particles with range of one or more lead plates, 
in addition to particles with range fess than one 
lead plate. The average number of observed par- 
ticles per star in the former group was 3.3, and 
in the latter group 5.5. 
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Fic. 5. Curves for 5-ray frequency, ionization, maximum 
range of é-rays, and residual range of mesotrons or protons 
as a function of the velocity of the primary particle. N/No 
and I/I¢ are, respectively, the frequency of é-ray produc- 
tion and the ionization, both of them relative to the value 
of these quantities for a primary particle with velocity =—c. 
Rmazx is the maximum possible range for a 5-ray. R is the 
residual range of the primary particle. For an a-particle, 
N/Noand I/I, are four times as large as indicated, Rmax is 
unchanged, and R is the same as for a proton. 


A block diagram showing the observed fre- 
quency-range distribution for particles originat- 
ing in nuclear disintegrations is given in Fig. 4. 
The relative numbers of short-range particles are 
too high since the solid angle subtended by a 
plate decreases as we proceed away from the 
origin of the star; the dotted lines of Fig. 4 
indicate the results of the corrections for this 
effect. The original range distribution of the star 
particles is distorted by their passage through the 
lead plate in which they originated. The observed 
distribution is approximately the same as the 
original distribution for particles of range greater 
than one lead plate, but the observed distribution 
gives no idea of the original distribution for par- 
ticles with range less than one lead plate. The 
total number of particles in the latter group can 
be estimated from the fact that 2/58 of the stars 
originated in the air between plates whereas the 
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7H. J. Bhabha, Proc. Roy. Soc. 164, 257 (1938). 
y 









stopping power of the air relative to the plates ip 
g/cm? was only 2/5800. We therefore conclude 
that there were roughly 100 stars (containing 
only short-range particles) that were not ob. 
served for each one that was observed ; Powell® 
calculated a value of 50 for one-cm lead plates, 
which would reduce to 35 for 0.7-cm lead plates, 
in contrast to the value of 100 found in the 
present experiment. However, his calculation was 
made by dividing the lead thickness by the aver- 
age range of the star particles; if Powell’s data 
are analyzed in terms of stopping power of the 
gas relative to that of the lead, we again obtain 
a value of 100 for the number of unobserved stars 
per observed star. The discrepancy is attributable 
partly to experimental uncertainty and partly to 
the assumption that the frequency of star pro. 
duction is proportional to the nuclear mass. Thus 
we conclude that the true range distribution of 
star particles includes 35 to 100 times the ob. 
served number for particles with range less than 
one centimeter of lead ; most of these unobserved 
short-range particles have ranges much less than 
one millimeter of lead. 

In many cases the nature of the particles can 
be inferred from their ionization, change in 
ionization, scattering, or production of 6-rays, 
The first three methods have been discussed 
many times, but quantitative information for the 
latter method has not been found in the litera- 
ture. Since low energy transfers predominate, the 
interaction due to spin is negligible and the col- 
lision probability for heavy particles’ with 
electrons is proportional to 


(1/B8?)(1—B°E/Emax)dE/E?, (1) 





























where E is the energy of the ejected electron, 
8 the velocity of the primary particle in units of 
c, and Emex is the maximum energy that can be 
transferred to the electron. For the type of col- 
lisions considered here, Emax is given by Eq. 
(1.5a) of Rossi and Greisen* with sufficient 
accuracy, 











Emax = 2moc*6?/(1—6*), (2) 


where myc? is the rest energy of the electron. The 
expression for the number of 5-rays per unit path 
with energy between Emin and Emax found by 
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integrating (1) and expressing Emax in terms of 


8 from Eq. (2) is 
oe i—s* 


B? Emin B? 








2myc?- B? 
—(1— £7) In | (3) 


Emin(1 — B?) 


In the present experiment 6-rays with ranges 
as short as one mm could be observed, and this 
fact leads to a value of 7 kev for Emin.* If we 
define Np as the value of N for B~1, then N/No 
represents the frequency of occurrence of 6-rays 
relative to the minimum frequency of occurrence 
(primary particle with 81). A plot of N/No as 
a function of 8 is given in Fig. 5. The ionization, 
the maximum possible range for a 6-ray, and the 
residual range of the primary are also plotted in 
Fig. 5. The value for No obtained by Hornbeck 
and Howell® is 0.011 per cm of air at 76 cm and 
15°C for Emin=12 kev. For Enmin=7 kev the 
value for No would be 0.011 12/7 or 0.019 per 
cm, since we see from Eq. (3) that No (the value 
of N for B™1) varies inversely with Ein. 

Of the 166 tracks of particles with range less 
than 0.7 cm of lead, 70 showed ionization corre- 
sponding to singly-charged particles with velocity 
greater than 0.75 c; such particles could be either 
electrons or mesotrons. The remaining 96 short- 
range particles ionized at least several times 
more strongly and therefore could not be elec- 
trons. Three of these particles were probably 
a-particles because of the large number of 6-rays 
occurring along their paths. Each of them had 
four or more 6-rays with range at least one mm 
in a track length of ~2.5 cm. Thus N/Nois ~85, 
which is a highly unlikely value for a proton or 
mesotron (Fig. 5). The rest of the particles were 
probably not a-particles since an a-particle with 
range <0.5 cm of lead would produce at least 
one 6-ray per cm. Shapiro!® has found that 90 
percent of the star particles registered by photo- 
graphic emulsions are protons, and the remainder 
probably are a-particles. Since such particles lie 
near the short-range group under discussion, 


*G. Fr. von Droste, Zeits. f. Physik 84, 17 (1933). 
sai} Hornbeck and I. Howell, Proc. Am. Phil. Soc. 84, 33 
 M. M. Shapiro, Phys. Rev. 61, 115 (1942). 
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there is good agreement with the present results 
in respect to the small number of a-particles. 
Among the particles with range greater than 
0.7 cm of lead, there were twelve that might 
have been electrons, as judged from scattering 
and ionization ; three of these twelve also multi- 
plied and this indicates that the three were 
probably electrons. However, the identity of the 
twelve is not certain even in the case of the 
latter group of three. Three of the long-range 
particles might be interpreted as a-particles, but 
they might also be protons. The remaining 65 
long-range particles were either protons or 
mesotrons. It is difficult to estimate how many 
of these long-range star particles are mesotrons 
rather than protons. Nielsen and Powell" photo- 
graphed the tracks of slow particles at 14,000 ft. 
with a technique that left virtually no uncer- 
tainty in the identification of mesotrons. They 
found that only about one percent of the heavy 
tracks observed at high altitudes was due to 
mesotrons. However, this does not necessarily 
indicate that most of the star particles are not 
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Fic. 6. Frequency of occurrence of stars in an energy 
band of width 20 Mev as a function of the total observed 
energy. The energy, which is calculated on the assumption 
that all the particles are protons, would be approximately 
one-half as great if the particles are mesotrons. 


mesotrons. In the first place, the range of a slow 
proton is ten times that of a mesotron moving 
with the same velocity, and hence the observa- 
tions of Nielsen and Powell indicate that the 
relative frequency of formation might be ten 
percent. In the second place, a large fraction of 
the slow particles occurs singly; that is, such 


uC, E. Nielsen and W. M. Powell, Phys. Rev. 63, 384 
(1943). 
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particles did not originate as part of what we 
have termed a star shower (although the dis- 
tinction may be based on energy alone). It is 
difficult to estimate the number of single slow 
particles from the present series of photographs 
because contamination a-particles traversed the 
full distance between the plates even when the 
path was inclined 45 degrees from the vertical. 
A lower limit, which may be several times too 
small, is 400 single particles that certainly were 
not a-particles. Anderson and Neddermeyer™ 
observed on Pikes Peak about twenty times as 
many single heavily-ionizing particles as stars. 
Consequently, it is quite likely that Nielsen and 
Powell observed very few star particles, and 
their results do not necessarily indicate that most 
of the star particles are not mesotrons. Anderson 
and Neddermeyer” identified two particles in one 
star (Fig. 10, reference 12) as probably protons; 
one particle in each of two stars (Figs. 12 and 13) 
is definitely a mesotron; and four other particles 
in one star (Fig. 12) are either mesotrons or 
electrons. Crussard and Leprince-Ringuet™ pho- 
tographed a star (Fig. 3, reference 13) containing 
two particles that are certainly protons. There 
are two other cases that have been reported in 
which the multiple production of identifiable 
heavy particles has been observed.:'5 The par- 
ticles have been identified as mesotrons in both 
cases. In several other reports penetrating par- 
ticles produced in multiples have been interpreted 
as mesotrons without any real justification for 
such an interpretation. 

Only one of the thirty-one stars containing 
merely short-range particles was observed to 
have been caused by an ionizing ray, but in 
several cases such a ray could not have been 
observed. A _ time-coincident cascade shower 
accompanied one of the stars; in this case a non- 
ionizing ray gave rise to the star. Eight of the 
stars containing long-range particles were initi- 
ated by ionizing rays and thirteen by non- 
ionizing rays. In two photographs time-coincident 
cascade showers were also present. In one of 


122C, D. Anderson and S. H. Neddermeyer, Phys. Rev. 
50, 263 (1936). 

13 J. Crussard and L. Leprince-Ringuet, J. de phys. et 
rad. 8, 213 (1937). 

4G. Herzog and W. H. Bostick, Phys. Rev. 59, 122 
(1941). 
15D). J. Hughes, Phys. Rev. 60, 414 (1941). 
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these cases the star was initiated by a pene. 
trating ionizing particle, which also appeared to 
initiate the shower through the medium of a 
knock-on electron; in the other case the nature 
of the ray that initiated the star could not be 
ascertained. 

Twenty of the stars consisting of short-range” 
particles contained only particles that were’ 
projected downward whereas the other eleven” 
contained both upward and downward projected” 
particles, the number of upward and downward 
projected particles being about equal in the 
latter case. On the other hand, only six of the 
stars containing long-range particles included 
upward projected particles, and the total 
number of such particles was only eight. 

The energies of the particles can be estimated 
from their range limits in lead® (see Fig. 4). The 
most energetic particles had energies greater than 
100 Mev if they were mesotrons or 200 Mev, if 
they were protons. A block diagram of the fre- 
quency versus total energy actually observed in 
the star is given in Fig. 6. The energies were cal-_ 
culated by assuming that all the star particles” 
were protons. Examples of the photographs of 
stars are reproduced in Fig. 7. 

In the above analysis, heavily-ionizing par- 
ticles that occurred singly and originated in the 
lead plates have been omitted because, in the 
first place, it was impossible in many cases to 
distinguish between a-particles from radioactive 
sources and heavily-ionizing particles ejected by 
cosmic rays, and, in the second place, it is not 
known that the single particles result from the 
same type of event. 





















SINGLE HEAVILY-IONIZING PARTICLES 





As mentioned in the preceding section, there 
were observed at least 400 heavily-ionizing par- 
ticles that were projected from the lead plates 
and that could not have been contamination 
a-particles. Most of these particles evidently 
were projected from lead nuclei by neutrons, i.e, 
the incident ray was non-ionizing and the event 
was not associated with a cascade shower. Eight 
of the heavily-ionizing particles were projected 
from the lead by penetrating particles, which 
lost only a portion of their energies in the process 
and then continued through the remaining lead 
plates. 
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Fic. 7. Examples of nuclear disintegrations or star showers. Upper left: the most frequent type of star. The particles 
have short range and heavy ionization. Upper right: a star containing long- and short-range particles, probably initiated 
by a neutron. Lower left: a similar star initiated by a charged penetrating ray. Lower right: a star containing several long- 
range particles. The track doubling in one section results from a residual clearing field. The right long-range particle 
produces secondaries, but it is probably not an electron. 
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In some cases heavily-ionizing particles, which 
originated outside the chamber or in the glass 
cylinder, passed nearly horizontally through the 
chamber, and twenty or more centimeters of 
track length could be observed. In several cases 
eight or ten 6-rays occurred, a number that cor- 
responds to about twenty-five times the number 
for B~1. These particles could not have been 
mesotrons since the range of such a mesotron 
would be only twenty centimeters (Fig. 5). They 
might be protons or a-particles. 

A number of penetrating particles were ob- 
served that terminated their ranges in the lead 
plates but originated outside the chamber. In the 
last 4000 of the pictures, 5 inches of paraffin were 
placed above a portion of the chamber in order 
to determine whether the particles were protons 
projected from hydrogenous material by neu- 
trons, as suggested by Powell.* A more detailed 
analysis of the photographs will be required 
before any conclusion can be drawn. 


CONCLUSIONS CONCERNING NUCLEAR 
DISINTEGRATIONS 


(1) The initiating particles are very seldom, if 
ever, electrons or photons. If they were electrons 
or photons with sufficient energy to produce the 
observed disintegrations, they would produce 
cascade showers in the lead plates above the star 
origin or they would be accompanied by other 
electrons or photons. In previous cloud-chamber 
experiments that seem to indicate a correlation 
between stars and cascade showers the cioud 
chamber was tripped by counters, in some cases 
with counter arrangements that selected multiple 
events only. For example, Anderson and Nedder- 
meyer” observed cascade showers in more than 
one-third of their photographs taken at Pikes 
Peak. Hence their observation of 30 cases of 
time-coincident showers in 113 photographs with 
heavily-ionizing particles does not indicate asso- 
ciation between the two types of events. The 
other published results also indicate no more 
association between cascade showers and stars 
or heavily-ionizing particles than would be 
expected of independent events. In the present 
work expansions were not counter controlled, and 
hence there was much less distortion in the 
relative frequencies of occurrence of the various 
events. Since eight percent of the photographs 
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showed cascade showers of four or more particles, 
one would expect four or five of the 58 photo. 
graphs of nuclear disintegrations to contain 
cascade showers ; the observed number was three. 
Powell® also noted little if any association be. 
tween showers and stars. On the other hand, 
Auger'* observed a real correlation between 
heavily-ionizing particles and extensive (Auger) 
showers. However, this does not necessarily 
indicate production of penetrating rays by the 
soft component since it is known that penetrating 
rays are also present in Auger showers. 

(2) The stars are probably produced by neu. 
trons and protons. Since the non-ionizing rays 
that initiate stars are not photons, we are left 
with neutrons as the reasonable alternative 
choice. The ionizing rays that initiate showers 
are neither electrons nor a-particles, but in this 
case there are two further alternatives: protons 
or mesotrons. However, the frequency of oc. 
currence of stars increases several times more 
rapidly than the intensity of the penetrating 
particles as we go to higher altitudes from sea 
level. Since most of the penetrating particles 
are mesotrons, it would be possible for the rela. 
tive number of protons to increase several times 
in going from sea level to 14,000 feet, whereas 
the relative number of mesotrons could not. 
Hence it seems most likely that the star-initiating 
ionizing particles are protons. 

(3) There may be no sharp line of distinction 
between low energy stars with particles projected 
in all directions and high energy stars with long. 
range particles projected in the general direction 
of the incident particle. Rather there appears to 
be a gradual change from one type to the other 
with increasing energy of the initiating particle. 
Examples of all types of intermediate cases were 
photographed (Fig. 7). The single heavily: 
ionizing particles originating in the lead plates 
are perhaps merely the limiting case for low 
energies. 

(4) The initiating particles are predominately 
neutrons for the lower energy stars and about 
equally divided between ionizing and non- 
ionizing particles (protons and neutrons) for the 
higher energy stars. Nearly all of the single 
















































16 P, Auger, R. Maze, P. Ehrenfest, and A. Freon, J. de 
phys. et rad. 10, 39 (1939). 
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heavily-ionizing particles were produced by 
neutrons. 

(5) The star particles are mostly mesotrons 
and protons. In the present work very few elec- 
trons or a-particles were observed in the high 
energy stars and less than half of the short-range 
particles could have been electrons. Examples 
have been quoted from other papers that gave 
proof of the presence of both mesotrons and 
protons. 

Two other conclusions, which are related to 
reported results from other investigations, should 
be mentioned. (1) It seems rather presumptuous 
to assume that an observation of the multiple 
production of long-range particles with either 
counters or a cloud chamber is necessarily an 
observation of the production of mesotrons, in 
cases where there is no evidence that the par- 
ticles are not protons. (2) Detailed analyses of 
stars observed in photographic emulsions that 
make use of conservation of momentum or of 
energy are likely to be meaningless with present 
techniques since particles with minimum ioniza- 
tion occur in a few of the stars and are not 
observed in photographic emulsions. 

The observed characteristics of stars are satis- 
factorily described in terms of the semiquanti- 
tative theory of nuclear collisions as discussed 
by Heitler.'? He not only suggests that all of the 
transferred energy of the incident particle might 
be dissipated in “‘heating’’ the entire nucleus, 
with subsequent emission of many particles, but 
he also suggests the possibility that the particles 


17 W. Heitler, Phys. Rev. 54, 873 (1938). 


near the edge of the nucleus receive the energy 
and escape before sharing their energy with the 
rest of the nucleus. Heitler’s estimates give 
2x10’ ev for the maximum energy of the 
emitted particles in the former process, when the 
energy of the incident particle is 210* ev. In 
the present as well as in previous work energies 
at least five or ten times larger than 2X10’ ev 
have been observed for star particles. This result 
might indicate that the initiating particles had 
energies ~10" ev or else that the latter of the 
processes suggested by Heitler was occurring. 

There appears to be a gradual transition from 
the low energy uncollimated type of star pro- 
duced mostly by neutrons to the high energy 
collimated type initiated by neutrons and pro- 
tons. It may be that stars of the latter type 
constitute examples of the production of meso- 
trons by protons such as is indicated by the 
experiments of Schein, Jesse, and Wollan.'* If 
this is so, it seems likely that protons are also 
produced in the same events. 

Hamilton, Heitler, and Peng’® have suggested 
a cascade process for the mutiple production of 
mesotrons. None of the present photographs 
gives evidence for such a process. 

The author welcomes this opportunity to 
express again his gratitude to Superintendent 
Frank A. Kittredge and other members of the 
staff of Yosemite National Park for their helpful- 
ness. 


18 M. Schein, W. P. Jesse, and E. O. Wollan, Phys. Rev. 
59, 615 (1941). 

1” J. Hamilton, W. Heitler, and H. W. Peng, Phys. Rev. 
64, 78 (1943). 
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¢ 
given by the empirical formula 


, 1 1 n 3 i 
v(n, j) = [2 —oo— o(Z) P+ aX(Z—o2)( =F) ++ 2 DZ, 
where oo and o2 are constants for a given level and 
¢(Z) = —0.098— 1.80 x 10°°Z + 1.40 x 10-42? 


for the K and L levels and equal to zéro for the M level. The values of the screening constants 
best suited to the above formula are also given in the paper. 


NUMBERS 3 AND 4 


On the Term Values of X-Ray Spectra 
Ts1-Minc Hu 
Department of Physics, National University of Chekiang, Meitan, Kweichow, China 
(Received October 8, 1943) 


It is shown in this paper that the term values of the K, L, and M levels of x-ray spectra can be 
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INCE Moseley’s discovery of x-ray spectra, 

many attempts have been made to calculate 
the term values of different levels. So far as the 
author knows, the most satisfactory formula is 
the one given by Pauling.’ The formula is 





1 1 n 3 
Hn, j= (Zon) + (2 —01)*{ : --) 
= n' j+4 4 


Z-S; 


n;* 


(1) 





+++--2D Z, 


By means of this formula and the S; values 
calculated by him, Pauling has shown that the 
values of a are nearly constant for a given level 
(independent of Z). He has also calculated the 
theoretical values of o) and a2 and found that 
they agree in general with the empirical values. 
But detailed calculations based upon his formula 
show that the agreement is still not quite satis- 
factory. The discrepancy is most prominent in 
the case of the K level. In order to see what cor- 
rection must be introduced into (1), the term 
values of the K, LZ, and M levels have been cal- 
culated by means of this formula with the screen- 
ing constants given by Pauling. The experimental 
values chosen for comparison are taken from the 
Handbuch der Experimental Physik.’ Since all the 
empirical values were determined by the crystal 
spectrometer, they have been multiplied by a 
numerical factor (1—0.002) on account of the fact 


(1930), pp. 180-191. 

2 Handbuch der Experimental Physik, Vol. 24, Part 2. 
A. B. Lindh, Réntgen Spektroskopie (1930), pp. 124-174, 
227-231. 





1L. Pauling and Goudsmit, The Structure of Line Spectra 


that the value of the electronic charge is 
4.803 X 10° e.s.u. instead of 4.774 10-"® e.s.u, 
The value of the fine structure constant used is 
1/136.9.3 

If the term values thus calculated are denoted 
by v° and the experimental values by v*, then 
the deviations Ay are given by v°—v*. Plotting 
Av against Z, we find the points obtained do not 
lie on smooth curves. But the deviations are 
certainly systematic in case of K and L levels. 
The precision of the calculations is limited by 
the presence of uncertain errors in the values of 
the size screening constants, especially for elec- 
trons in the incomplete shell. To overcome this 
difficulty the values of the wave numbers (in 
Rydberg units) for K and L lines have been 
calculated also and compared with empirical 
values. In this case, only the screening constants 
for electrons of a single complete shell enter into 
the calculations. The values of Av thus obtained 
change continuously with Z and can be joined 
with smooth curves (see Figs. 1 and 2). The 
similarity between the two curves naturaily 
suggests that there may be some simple relation 
between them. Further investigations show that 
they are connected by the relation AvxK,’=3Avly 
+0.15Z.4 The screening constants used in the 
calculations are given in Tables I and II. 

Encouraged by the simplicity of the relation 
and the smoothness of the curves, attempts have 

























’ Compton and Allison, X-Rays in Theory and Experi- 
ment, pp. 697-703. 

* According to Sommerfeld’s nomenclature, the Ka’ line 
corresponds to the transition K-—>Ly and the Lg line 
corresponds to the transition Ly;—~>Mty. 
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Fic. 1. Avx-, against Z. Equation of the solid curve is given by (3). 
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Fic. 2. Avis against Z. Equation of the solid curve is given by (4). 


been made first to smooth up the correction 
curves for different levels by making some 
alterations of the size screening constants® 
and then to introduce some simple corrections 
into (1) to bring it in better agreement with 
experiment. The results of such attempts may 
be summed up as follows: 

(a) The values of different terms may be given 


’ The alterations of the values of size screening constants 
are quite arbitrary; they are given only to illustrate that 
the seeming irregularities of the values of Av for different 
levels can be eliminated by using suitable size screening 
constants. 


by the formula 


1 
v(n, j) =—[Z—o0— 9(Z) P 


n* 


1 n ; 
+—08(Z—o)'(——— , 


n' j+i 4 


where 


¢(Z) = —0.098 — 1.80 X 10-*Z + 1.40 X 10-*Z 
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Fic. 3. S against Z. Dotted lines indicate the values given by Pauling. 


for the K and L levels and equal to zero for the 
M level. 

(b) The values of oo and a2 best suiting the 
above formula are given in Table III and that 
of S may be read from Fig. 3. 

Assuming the validity of (2) the values of Avx,’ 
and Avzg can be taken as the difference between 


TABLE I. Screening constants. 








2p 3s 3p 





4.52 10.9 13.2 








TABLE II. Screening constants. 
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the values of the wave numbers given by (1) 
and (2). We get, therefore, 


Avka’ = (1.5Z2+2.2)¢(Z)+0.005Z+1.2 (3) 


and 
Avrg=(0.5Z —3.4)¢(Z) —0.005Z—1.2. (4) 


The constant terms in these formulas are intro- 
duced by the difference of the values of the size 
screening constants for L and M shells used in 
(1) and (2) (see Fig. 3). These relations are repre- 
sented by the solid curves in Figs. 1 and 2. All 
the points lie close to these curves. As a further 
illustration of the validity of (2), the term values 


TABLE III. Values for formula (2). 








oo 0.6252 6.80 
C2 0.17 3.6 








* This value was assigned by considering the structure of the two K 
electrons and the nucleus to be helium-like. 
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5.89 6.79 20.0 
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TERM VALUES OF X-RAY SPECTRA 


TABLE IV. Calculated and observed term values. 














* 





Zz VK VK" vK* PL YL YL ¥My ¥My ¥My 
|< 

30 711.7 710.3 711.2 77.1 76.9 77.8 0.9 0.9 2.0 
35 992.6 990.6 990.6 117.7 117.5 117.3 5.2 5.2 5.2 
40 1325.7 1323.1 1322.4 170.0 169.7 168.5 13.4 13.4 12.3 
45 1709.6 1706.2 1707.5 231.2 230.7 231.3 22.1 22.1 23.6 
52 2345.0 2340.3 2340.1 340.3 339.6 340.2 42.4 42.3 42.6 
56 2756.4 2750.9 2752.6 414.3 413.5 413.7 57.6 57.5 56.7 
60 3214.2 3207.8 3203.1 495.5 494.5 494.7 72.5 72.4 71.8 
66 3972.5 3964.6 3972.2 632.2 630.9 630.6 95.8 95.6 94.8 
70 -- — — 735.4 733.9 733.7 112.8 112.6 112.3 
74 5113.8 5103.6 5109.9 850.6 848.9 849.2 133.7 133.4 133.4 
78 5764.0 5752.5 5763.1 978.7 976.7 977.1 156.4 156.1 157.0 
82 6463.0 6450.1 6471.1 1121.9 1119.7 1118.1 183.0 182.6 183.1 
90 8073.5 8057.4 8066.8 1451.5 1448.6 1447.3 244.9 244.4 244.0 














of K, L, and M levels have been calculated from 
(2) and given in Table IV together with the 
empirical values chosen for comparison. From 
this table it is easily seen that the term values of 
Ln and My levels given by (2) are quite close 
to the empirical values. Since either (1) or (2) 
can account for the doublet structure of x-ray 
spectra successfully, it follows naturally that (2) 
can be applied to calculate the term values of 
any of the Z and M levels with results agreeing 
with experiment. In the case of the K level the 
agreement is less satisfactory for heavy elements. 
Since (2) gives the values of wave numbers for 
K lines and L levels close to the empirical values, 
it may be inferred from the combination principle 
that the term values of the K level given by (2) 
must be close to the actual values, and the 
deviations may arise from some systematic 
errors of experiment.In Table IV vx, vz, and 
vMy are the experimental term values. 


ve’ =(1-— 0.002) vx, 


and similar equation for y111" and yMy’. vL11°, yMy* 
and vxk* are the values calculated from (2). 

The above corrections may arise from manifold 
sources, but the electronic interaction is the most 
probable one as is usually the case in spectros- 
copy. It is interesting to note that the correction 
term ¢(Z) is the same for the K and L level. 
Attempts have been made to use the same func- 





tion for the M level by changing the value of 
a» and S, but the result is less satisfactory. 
Evidently, the required correction may be intro- 
duced into (1) in other ways. We may, for 
example, write the correction formula in the form 





1 a? n 3 
Hm, j)=—[Z—ooP+= 2-01} ) 
n n 





i+} 4 
Z-S; 
+---—-230 2; . —Av(n, Z), 
i nN; 


where Av(n, Z) may be taken as the difference 
between (1) and (2) or expressed in any other 
form. But whatever may be the form of the cor- 
rection term, two facts are certain. Firstly, 
Av(i,Z) and Av(2,Z) are connected by the 
relation 


Av(1, Z) =4Av(2, Z)+aZ+b. 


The values of a and b are affected by the screen- 
ing constants. They can be made to vanish if 
suitable screening constants are used. Secondly, 
the corrections for the M level are small and 
masked by the errors in the values of the screen- 
ing constants. 

In conclusion the author wishes to express his 
sincere thanks to Professor K. C. Wang. This 
work was assigned by him and has been done 
under his direction. 
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The Angular Distribution of Alphas from Li’(), a)« 
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NUMBERS 3 AND 4 


The tracks formed in photographic emulsions by alphas from the reaction Li’(p, a)a at 
various angles from the beam were observed at bombarding energies E between 400 and 900 kev. 
The angular distribution proportional to 1+ A(Z£) cos? @, as found by Young, Ellett, and Plain 
at lower voltages, was also observed here, and a maximum of A(£) was found at 675 kev. 
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NE of the simplest reactions yielding disin- 
tegration products whose angular distribu- 
tion is of some theoretical interest is: 


Li’+H!'=He'+He'+17 Mev. (1) 


The alphas from this reaction have a range of 
about 8 cm in standard air. This reaction has 
been studied by several observers,! and especially 
carefully by Young, Ellett, and Plain,? who 
employed proton energies up to 400 kev and 
showed that in this energy range the intensity 
as a function of angle 7(@) was very well repre- 


sented by 
I(@)~1+A(E£) cos? (6). (2) 


Here @ is the angle between the incident proton 
and the observed alpha in the center-of-mass 
coordinate system, and A(£) is a factor depend- 
ing only on the energy E of the proton. Their 
values of A(E) are represented by the small 
circles in Fig. 2. Although the statistical uncer- 
tainties were too great to allow a definite con- 
clusion on this point, the curve seemed to level 
off as it approached their maximum bombarding 
voltage, 400 kev, which suggested the existence 
of a maximum nearby. 

Critchfield and Teller,* in a theoretical discus- 
sion of this reaction, correlated this apparent 
maximum with a resonance at about 400 kev 
indicated by the yield curve‘ together with a 


* Now at S. A. M. Laboratories, Columbia University, 
New York, New York. 

t Now at the Western Electric Company, Baltimore, 
Maryland. 

1F, Kirchner, Physik. Zeits. 34, 785 (1933); J. Giar- 
ratana and C. G. Brennecke, Phys. Rev. 49, 35 (1936); 
H. Neuert, Ann. d. Physik 36, 437 (1939). 

2 Young, Ellett, and Plain, Phys. Rev. 58, 498 (1940). 

3C. L. Critchfield and E. Teller, Phys. Rev. 60, 10 (1941). 

* Rumbaugh, Roberts, and Hafstad, Phys. Rev. 54, 651 
(1938); Herb, Parkinson, and Kerst, Phys. Rev. 48, 118 
936° Heydenburg, Zahn, and King, Phys. Rev. 49, 100 
1936). 
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theoretical treatment of the penetration factor, 
They showed that both of these observations 
could be explained by simple assumptions involy. 
ing only a broad S and a narrow P state of the 
compound nucleus Be’. 

The protons with energies up to 900 kev were 
provided by the pressure insulated Van de Graaff 
generator in the Johns Hopkins University 
Physics Department. The target chamber as. 
sembly is shown diagrammatically in Fig. 1. The 
beam entered the chamber through a slit 2 mm 
high by 6 mm wide. The inclination of the target 
block made the target spot about 5 by 6 mm, 
The thin lithium targets used were made by 
evaporating metallic lithium onto a hexagonal 
target block made of copper, which could be 
turned to present any one of its faces to the beam 
or to the furnace. The evaporation was done in 
vacuum with a simple lithium stove as shown. 
The metallic lithium had to be packed very 
tightly into the cavity of the furnace to insure a 
steady and slow evaporation. With small pieces 
of lithium merely dropped into the stove, much 
more power was needed, and the pieces disin- 
tegrated in an explosive fashion. A rough check 
by means of the 440-kev gamma-ray resonance 
indicated that the target thickness was usually 
equivalent to about 20 to 60 kev. 

The bombarding voltage was measured by 
magnetic deflection of the proton beam. The 
deflecting field was held constant and _ voltage 
was indicated by the angular displacement of the 
target chamber necessary to put the proton spot 
on the entrance slit. The deflection was calibrated 
by the 440-kev gamma-resonance® of Li and 


































5 A description will be published in The Review of Sciew 





* Hafstad, Heydenburg, and Tuve, Phys. Rev. 50, 5 
(1936). : 
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Fic. 1. The angular distribution of alphas from Li? (p, a)a. Target chamber. 








verified by the 334-kev and 660-kev resonances*’ 
of F. Two values of the magnetic field were used 
to secure convenient deflections over the whole 
range of energies used. The accuracy of this 
method of measurement was limited by erratic 
fluctuations in the position of the undeflected 
beam position. The uncertainties in voltage 
arising from voltage fluctuation and inaccurate 
measurement are expressed by the width of the 
circles in Fig. 2. 

The alphas were detected simultaneously at 
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Fic. 2. A(E) as a function of bombarding energy. Small 
circles and broken curve: data and curve of Young, Ellett, 
and Plain. 


7 Bennett, Herb, and Parkinson, Phys. Rev. 54, 398 
(1938). 
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five different angles by means of photographic 
plates placed opposite each of the five thin 
aluminum windows as shown in Fig. 1. Since the 
range of these alphas in air is about 8 cm and the 
air equivalent of the 0.0005” aluminum foil is 
about 2.5 cm, there is a residual range in the 
emulsion quite sufficient to produce a track of 
satisfactory length. 

The angle of incidence was made 45° as a 
mean between normal incidence, in which the 
track would appear as just one point, and grazing 
incidence, in which the track would be lost in the 
surface fog. The plateholder was designed to fit 
around the target chamber with an aperture 
opposite each thin-foil window. The apertures 
were closed by a shutter to allow removal of the 
plateholder for loading. The plates used were 
Eastman fine grain Alpha Particle Plates 1” by 
3’’. It was found necessary to store the plates in 
a dry and fairly cool place to avoid fogging 
which, although still invisible without the micro- 
scope, would interfere with the counting of 
alpha-tracks. After exposure the plates were 
developed in a specially constructed rack of 
glass rods, and after drying the exposed parts 
were protected with a cover glass. 

The tracks were counted under a microscope 
with dark-field illumination and 400-fold mag- 
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Fic. 3. Angular distribution of alphas from Li’(p, a)He‘ at 
several bombarding voltages. 


nification. The ocular had a diaphragm with a 
rectangular opening traversed by a cross-hair. 
The plate was moved by a mechanical stage, and 
the number of tracks passing the cross-hair in a 
swath of fixed width was counted as the stage 
moved through unit distance. About a thousand 
tracks were counted on each plate. 

The photographic method was chosen in 
preference to electrical counting methods because 
of its simplicity and because the simultaneous 
exposure of five plates required less operating 
time of the generator than the customary elec- 
trical counting at only two windows simul- 
taneously. This advantage is, however, partially 
offset by the tedious task of counting the tracks 
at leisure. 

The data so obtained were transformed? to the 
center-of-mass coordinate system. This intro- 
duces modifications up to almost 20 percent in 
spite of the great energy liberated. 


§ Haxby, Allen, and Williams, Phys. Rev. 55, 140 (1939). 


Figure 3 shows the intensity J(@) plotted 
against cos? @. The results are shown for several 
different bombarding energies. In most cases it 
was not possible to draw a straight line through 
all five of the circles (whose radii equal the 
standard deviation) which represent the data, 
In some cases the best straight line passes rather 
far outside of some of the circles, and there is 
only a small probability that, at the same time, 
Eq. (2) applies and no experimental error has 
occurred other than statistical fluctuations, 
There is a tendency for one or two of the circles 
to fall below a straight line through the other 
four or three, but there is no apparent regularity 
in these deviations. We prefer to attribute them 
to experimental error (such as non-uniformity of 
plates or of development which might possibly 
persist in spite of precautions), rather than to 
deviations from (2). In particular, the lowest 
circle in the set at the highest voltage falls 
seriously below the line determined by the other 
four, and we have arbitrarily ignored it. In all 
cases, the slopes in Fig. 3 were determined by 
personal judgment. Through each set appear 
three lines to represent what we consider the 
most probable and the two extreme admissible 
slopes. 

Figure 2 shows A(£) as a function of bom- 
barding energy E. The circles represent values of 
A (the slopes of the lines in Fig. 3) which were 
regarded as the most probable. The length of the 
vertical lines represents the estimated range of 
admissible slopes. It can be seen that the curve 
passes through a moderately sharp maximum and 
then comes down rather rapidly as predicted.’ 
However, the maximum occurs not around 400 
kev, but at about 675 kev (which corresponds to 
about 590 kev in the center-of-mass coordinate 
system). Our data are not sufficiently detailed to 
confirm or reject the inflection at 400 kev sug- 
gested by the data of Young, Ellett, and Plain, 
and indicated by the dotted line in Fig. 2. 

Weare very grateful to Professor J. A. Bearden 
for his advice and direction during the early part 
of this work. 











tted 
eral 
2s it 
ugh 
the 
ata, 
ther 
e is 
ime, 
has 
ons, 
‘cles 
ther 
rity 
hem 
y of 
ibly 
1 to 
vest 
falls 
ther 
| all 
| by 
ear 
the 
ible 


om- 
s of 
vere 
the 
e of 
irve 
and 
ed? 


s to 
late 
1 to 
sug- 
ain, 


den 
art 












PHYSICAL REVIEW 


VOLUME 65, NUMBERS 3 AND 4 


The Angular Distribution of Alphas from F'*(p, «)O" 


H. H. Rosst anp C. D. Swartz* 





FEBRUARY 1 AND 15, 1944 





Johns Hopkins University, Baltimore, Maryland 


(Received November 12, 1943) 


The angular distribution of alphas from the reaction F"(p, a)O™ has been investigated at 


bombarding voltages between 500 and 900 kev. The distribution was found to depend markedly 
on bombarding energy, contrary to the result found by McLean, Ellett, and Jacobs at lower 
voltages. 


HE angular distribution of alphas from the 
reaction"? 


F19+ H!=O0'*+ He!+7.95 Mev (1) 


has been studied by McLean, Ellet, and Jacobs 
in the energy range 330 to 435 kev. The variation 
with voltage observed in this range was so slight 
that they interpreted it as hardly significant in 
view of the statistical error. Gerjuoy* has shown 
theoretically that these results could arise from 
three overlapping levels of Ne?’ with J=0, 1, 
and 2. These states would lead to no powers 
higher than the third when the intensity is 
expanded in powers of cos 06. 

We have investigated the reaction (1) in the 
energy range 500 to 900 kev. This was done in 
conjunction with the experiments reported in the 
preceding paper, and the apparatus and the 
methods employed were essentially the same. 
The targets used were slabs of fluorite (CaFe2) 
about 3” thick mounted on the copper target 
support. Attempts to use a thin target of NiFs, 
made by electrolizing HF with a nickel electrode, 
were unsuccessful with proton currents through 
the defining slit of the order of one microampere 
available. 

The use of a thick target seemed justified at 
first by the combination of two circumstances: 
the very rapid rise of the total yield with voltage 
and the slow variation of the distribution which 
was anticipated by extrapolation from the earlier 
experiments.’ Our results show that this antici- 
pation was false, however, and the rapid vari- 


* Now at S. A. M. Laboratories, Columbia University, 
New York, New York. 
1 W. E. Burcham and C. L. Smith, Proc. Roy. Soc. A168, 
176 (1938). 
1 + Fowler, and Lauritsen, Phys. Rev. 59, 253 
* McLean, Ellett, and Jacobs, Phys. Rev. 58, 500 (1940). 
‘E. Gerjuoy, Phys. Rev. 58, 503 (1940). 






ation of yield is not by itself sufficient to justify 
substitution of thick for thin targets, except for 
the sake of a rough indication of trends. A thin- 
target result may, of course, be obtained from 
sufficiently complete thick-target data by dif- 
ferentiation with respect to voltage. But our 
thick-target data are not complete enough to 
allow this to be done accurately. 

Between 500 and 900 kev the yield curve? of 
this reaction is not a steadily and rapidly in- 
creasing function, but has two pronounced 
maxima. Therefore, for certain values of the 
bombarding energy V, the production of alphas 
will be less abundant at the surface of the target 
than at some distance below it, where the 
protons have a lower energy V. The number of 
collisions of energy between V and V+dV will 
be proportional to the thickness of the layer in 
which the corresponding deceleration takes place. 
This distance is equal to dV divided by V’, the 
slope of the energy-range curve at energy V. 
This must be multiplied by y(V), the yield at 
V and g(V), the charge decelerated from V+dV 
to V, to get a measure of the number of alphas 
produced between V+dV and V. In the experi- 
ment the bombarding voltage varied during an 
exposure, and a record was made of the incident 
charge as a function of V;. By integrating over 
all charges incident at values of V; greater than 
V, one obtains the total charge g(V). Thus the 
number of alphas arising from collisions of energy 
between V and V+dV during an entire exposure 
is proportional to dV multiplied by 


I(V)=q(V)y(V)/V’. (2) 


The results of three exposures at the lower ranges 
of bombarding voltages are shown in Fig. 1. The 
relative numbers of alphas produced by col- 
lisions of various energies are indicated by the 
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Fic. 1. Angular distribution of alphas from F'*(p, a)O"* 
at low bombarding voltages. The broken line represents 
the data of McLean, Ellett, and Jacobs between 330 and 
435 kev. The inserts give in arbitrary units the relative 
number of a’s as a function of voltage for each set. 


functions J(V) shown in the inserts to the figure. 
In order to verify that the striking variation of 
the angular distribution with voltage was not 
due to some erratic experimental error, three 
exposures were made at the higher voltage range 
with two different fluorite crystals as target. The 
results are in satisfactory agreement, as shown 
in Fig. 2. The insert in this figure gives 7(V) for 
one of the exposures, the other two having quite 
similar curves of J(V). 

The ‘‘points” are plotted as circles, squares, 
and triangles (the dimensions of these figures are 
approximately equal to the standard deviation). 
We were unable to check the data of McLean, 
Ellett, and Jacobs. Every attempt to secure data 
at voltages below 500 kev produced very heavy 
general fogging of the region of the plates exposed 
to the alphas. The cause of this fogging appears 
to have been some soft radiation from the target, 
since it could penetrate the 0.0005” aluminum 
window but not the 7,” brass wall of the target 
chamber. It should be pointed out, however, that 
the data of McLean, Ellett, and Jacobs do not 
extend beyond cos @=0.70, so it is possible that 
even at their voltages the curve turns downward 
beyond that point. 

Particular pains were used to check possible 
sources of error. Evidence that the photographic 


method of determining relative intensities is not 
at fault is furnished by the consistency of the 
three sets of data at higher voltages, and also by 
the good internal consistency of the data on the 
angular distribution of alphas from the Li’(p, a)q 
reaction.’ The possibility of extraneous alphas 
produced by contaminations in the target ma- 
terial must be considered. In the arrangement 
used, an alpha can enter the emulsion only after 
traversing 0.0005” of aluminum and about 1 cm 
of air. Since the air equivalent of aluminum is 
about 1800, this means that only alphas whose 
range in air is at least 3.3 cm can enter the 
emulsion at all. 

We consider that a range of at least 2-cm air 
equivalent in the emulsion is necessary in order 
to assure recognition of a track, so this particular 
apparatus will detect only alphas at least 5.5-cm 
range in air. Fluorine alphas, whose range is 
about 6 cm in air, are already difficult to count. 
Lithium and fluorine are the only light elements 
taking part in a (p, a) reaction which produce 
alphas of over 5-cm range. Precautions were 
taken to eliminate lithium contamination and, in 
addition, any large proportion of the long range 
lithium alphas would have been quite noticeable 
in examining the plates. The only other reaction 
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Fic. 2. Angular distribution of alphas from F'*(p, a)O" 
at high bombarding voltages. The broken line represents 
the data of McLean, Ellett, and Jacobs between 330 and 
435 kev. The insert gives in arbitrary units the relative 
number of alphas as a function of voltage for one set; the 
other two are quite similar. 


’ Swartz, Rossi, Jennings, and Inglis, Phys. Rev. 65, 80 
(1944). 
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offering any reason for suspicion is B*(p, a) Be’. 
Although these alphas have a range of only 4.4 
cm in air, the cross section for this reaction is 
larger than that for fluorine by a factor of the 
order of 1000. Enough long range stragglers 
might possibly occur to produce an appreciable 
effect if a large amount of boron were present. 
A rough spectroscopic test revealed no boron. 
The great strength of a boron line on a com- 
parison spectrum made us believe that the pro- 
portion of boron in the fluorite was almost cer- 
tainly less than one part in a thousand, and 
probably not more than one part in ten thousand, 
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but we have no definite basis for stating an upper 
limit. 

The fact that the targets from two separate 
pieces of fluorite (one of which was definitely 
pink, the other white) which would not be 
assumed to contain equal boron contamination, 
gave results identical within experimental error, 
seems to indicate that the distribution observed 
is indeed due to fluorine alphas only. 

We are very grateful to Professors D. R. Inglis 
and J. A. Bearden for direction and advice, and 
to Mr. Burridge Jennings, who constructed 
almost all the apparatus used. 
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Proton penetration factors for Li7(p; a)a are calculated under several assumptions of nuclear 


potential. The status of the Critchfield-Teller theory is briefly discussed in the light of these 
calculations and the angular distribution found by Jennings, Swartz, and Rossi. Measurements 


at higher energies are called for. 


DETAILED knowledge of Li’(p; a)a, one 
of the most widely studied of nuclear reac- 
tions, is quite important for nuclear theory. 
Since it is one of the simplest reactions, calcula- 
tions based on various detailed assumptions of 
nuclear forces can be carried through and 
checked against observation. Although experi- 
ments do not give directly the energy levels of 
light nuclei, it is convenient to divide the cal- 
culations into two stages, viz., the calculation of 
energy levels based on detailed assumptions of 
nuclear structure, and the calculation of excita- 
tion and angular distribution functions, based 
on these levels, but not involving very detailed 
assumptions of nuclear model. Hence a check of 
the second part of the calculation with experi- 
mental results can be used to infer the actual 
energy levels. 

Critchfield and Teller! have treated the theory 


* Now at the National Bureau of Standards, Washing- 
ton, D. C. 

'C. L. Critchfield and E. Teller, Phys. Rev. 60, 10 
(1941). 










of the angular distribution and excitation func- 
tions in the Li’(p; a)@ reaction by means of the 
usual perturbation theory formulation. More 
recently Swartz, Rossi, Jennings, and Inglis® 
have extended the range of angular distribution 
measurement past 400 kev to about 900 kev. 
They found that the results could be fitted by a 
factor of the form 1+A(E£) cos?@ over the 
entire range of bombarding energies thus far 
measured. The maximum value of A(E£) was 
found to occur at 675 kev rather than near 400 
kev as anticipated on the basis of Young, Ellett, 
and Plain’s data.* The later experiments were 
not sufficiently good in the neighborhood of 400 
kev to exclude the possibility of a secondary 
peak there, though they made this seem unlikely. 
A secondary peak there would be incompatible 
with the assumptions of Critchfield and Teller. 

Most of the energy dependence of the matrix 

?C. D. Swartz, H. H. Rossi, B. Jennings, and D. R. 
Inglis, Phys. Rev. 65, 80 (1944). 


’V. J. Young, A. Ellett, and G. Plain, Phys. Rev. 58, 
498 (1940). 
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elements involved in the perturbation calcula- 
tion is due to the Coulomb repulsion of the 
bombarding proton. Sigce only P protons are 
involved in this reaction,‘ it seems legitimate to 
factor the matrix elements into the so-called 
penetration factor for P protons impinging on 
Li’, which contains most of the energy de- 
pendence, and the “internal yield” factor, which 
is relatively constant. Critchfield and Teller! 
pointed out that the maximum in A(E£) would 
be expected to occur at lower bombarding energy 
than the maximum in the internal yield at right 
angles, as deduced in reference 1 from the data 
of Rumbaugh, Roberts, and Hafstad® and 
others. Actually the opposite was found, and the 
discrepancy is now accentuated in the light of 
the more recent data, by the shift of the maxi- 
mum of A(E£) to a higher energy than antici- 
pated. 

Since the internal yield is equal to the 
measured yield divided by the calculated pene- 
tration factor, or Y/P, it seems reasonable to 
examine the calculation of the penetration 
factor for stability of the result under reasonable 
changes in the details of the model used. The 
dispersion formula will be most accurate if the 
unperturbed wave functions used cause rapid 
convergence of the perturbation calculation. A 
proper choice will make the matrix elements 
V(r;m, p) occurring in the second-order terms 
as small as possible. Now 


V(r; n, p) = S ver VarpWan¥p,4r, 


where yc, is the wave function of the compound 
nucleus as calculated from the Hamiltonian of 
Be®, except that it reduces rapidly to zero at 
about r=7fo. Wan is the wave function of the Li’ 
nucleus as calculated from its Hamiltonian, H,. 
¥p, is the wave function of the incident proton, 
calculated with the Hamiltonian 7,4p+ Up. The 
Hamiltonian for the complete system may be 
written H=H4+T74pt+Upt+Vap and Vap is 
the perturbation potential. 

The choice of Up which minimizes V(r; n, p) 
is not the average potential of a proton in the 
Li’ field, which is attractive at small distance, 
but rather a repulsive potential. This has a 


*E. D. Courant, Phys. Rev. 63, 219A (1943). 
5 L. H. Rumbaugh, R. B. Roberts, and L. R. Hafstad, 
Phys. Rev. 54, 657 (1938). 


reasonable interpretation, viz., that the proton 
does not exist as such inside the nucleus, and 
hence its wave function should die out there.* 
Since H is fixed, Up+Vap is not subject to 
choice. If Up is chosen as a large repulsive poten- 
tial, Vap must be chosen large to compensate, 
There is thus an optimum Up which makes 
V(r; m, p) smallest. For this Up, yp, is small 
inside r=ro. Since wc, reduces to zero rapidly 
at about r=7o, the important contribution comes 
from the region about 79, and hence the de- 
pendence of V(r;m,p) on bombarding energy 
is closely the same as that of (|Wpp|),-r, since 
the other factors are independent of bombarding 
energy. 

By use of Yost, Wheeler, and Breit’s’ tables 
of Coulomb wave functions, values of the relative 
penetration factor P~(F;*/p*),.,, were cal- 
culated for several different values of ro and for 
several different choices of Up. 

(A) With Up a constant repulsive potential 
U> for r<ro, and a Coulomb potential for r>r, 
(cut off to zero at large distance for convenience 
of calculation), the result was not significantly 
dependent on the choice of Uo for the several 
trial values 785, 1500, 3000, and 15,000 kev. 

(B) Calculations were also made by using for 
Up a Coulomb potential down to r=0, again cut 
off at large distance. 

(C) Penetration factors arising from Kapur 
and Peierls’ * treatment of the dispersion for- 
mula, i.e., 1/pG,? in the notation of reference 7, 
were also calculated for the case of Up a Coulomb 
potential down to r=0. As Breit pointed out,* 
1/pG,’ is almost the same as F;*/p* as long as 
¢$;6;=1, which is true in this region of bombard- 
ing energies. 

With ro=4x10-" cm, a maximum was ob- 
tained in all cases, below A (EZ) max., independent 
of the details of the nuclear model used. For ry 
in this neighborhood the values for the right 
angle yield’ Y would have to be badly in error 
to make the maximum of Y/P occur at higher 
energy than the maximum of A(Z£). As f¢ is 
assumed larger, differences appear in the shape 


*H. A. Bethe, Rev. Mod. Phys. 9, 91 (1937). 

7™F. L. Yost, J. A. Wheeler, and G. Breit, Terr. Mag. 
40, 443 (1935). 

*P. L. Kapur and R. Peierls, Proc. Roy. Soc. 166, 277 
(1938). 

*G. Breit, Phy. Rev. 58, 506 (1940). 
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TABLE I. Internal yield at right angles, calculated from 
the observed data for various nuclear radii by several 
methods of approximation. Values listed are log Y/P plus 
a constant selected to make yield and penetration factor 
match at E=225 kev. 


—— 


barding energy (kev) 
——— | hf 
» (K.E. in center of mass system) 
— 083 622 494 393 








Calculation (A) 
I 0.40 0. 
II 0.38 0.36 
Ill 0.32 0.33 
IV 0.29 0.32 


38 


Iculation (B) 
ry 35 


0 

I 0.33 ~— (0.34 

ll 0.30 = 0.32 
IV 0.28 0.31 


Calculation (C) 
I 0.50 
Il 0.47 


III 0.38 
IV 0.34 


0.36 
0.36 
0.33 
0.33 


0.48 
0.45 
0.38 
0.35 


I. ro=6.17X10-" cm 
Il. ro=5.50X 10-8 cm IV. 


III. r9=3.80X 10-" cm 
r9=3.09X 10-8 cm 











of Y/P depending on the type of calculation 
used. (See Table I.) 

Calculation (B) apparently still shows a 
maximum in the same region as before, but in 
(A) and (C) the maximum, if it exists, is moved 
out of the range of energies thus far used. One 
may take these results to indicate either that ro 
must be in the neighborhood of 4X 10—" cm, as 
is customarily assumed, or that this perturbation 
calculation is not particularly trustworthy, as 
Breit® has remarked. 

A closer examination also reveals that the 
first choice above is not necessarily at odds with 
the Critchfield-Teller theory. It! leads to the 
result: 


m4 


where 





2(€—€o)8 cos 6 cos & =] 
1+é 


68 cos 6 cos mene) 
2 
+cos ite? ‘ 





€9= —tan 6+(u/2) cos 2 
and u=8/cos 6 cos £, and also 
68 cos 6 cos £(€9—e) 
~ 1+e?+2(e—e0)B cos 5 cos §-+ 26" 





A(E) 


with A(¢€)=0. Here é is the angle denoting ad- 
mixture of triplet and 6 exp (74) is a complex 
number depending on the wave function of the 
compound nucleus, v1z., 


B exp (15) = 104aeI"/2aol's. 


From this result it follows that A(Z£) max. 
=A (€m) = —36*/(6?+ue,) and o max.=o(e,) 
~1+?/ue,, where €,,is the position of A(£) max. 
and ¢, is the position of o max., and 


e=(E—E2)/3T 2. 
For u>0, €m=€o— (1+26?+ €?)#<0 and 
én= (€o—u) +[(eo—u)?+1 ]'>0, 


leading to €,>€», contrary to observation. 

On the other hand, the sign of u is merely a 
question of the phase of the triplet admixture 
in the wave function of the compound nucleus, 
and for “<0, é€m=eo+[1+26+e?]!>0 and 
€n=(€o9—u) —[(e9—u)?+1]'<0. In this case 
ém > €, as Observed. Also 


Em — €9 = (14+28?+ &?)*>0 


and hence € is below ¢,,. Therefore the observed 
zero in A(£) in the neighborhood of zero bom- 
barding energy occurs at ¢€. There is little dif- 
ficulty in finding sets of consistent values for 
t, 8B, 6, and Ez fitting the observed separations 
€m — €9 and €,—€,. Their direct calculation, how- 
ever, hardly seems feasible, since detailed con- 
siderations of nuclear model are involved. It 
does seem reasonable to assume, however, that 
X2=X, sin +X, cos —, the combination of 
incident triplet and singlet states with J=2 
which corresponds most strongly to the com- 
pound state involved in the alpha-decay, will 
not be more triplet than singlet, so that | £| <45° 
and cos 2é>0. 

Measurements at higher energies are needed. 
The present expectation is that A(E) goes to 
zero asymptotically at higher energy. Negative 
A(E), at least below the region where higher 
levels of Be** become involved, would force a 
revision of the theory. 

I gratefully acknowledge the constant en- 
couragement given and the kind interest shown 
in this work by Dr. D. R. Inglis, under whose 
direction it was done. 
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The differential equations governing the probability distribution of events distributed over 
a multidimensional domain are derived. They are a generalization of the equations governing 
the probability distribution in a time series. Let f,(x, y)dxdy be the chance that an event will 
occur in the domain dxdy when r events have already occurred in the domain (0, x; 0, y). The 
equations which are derived then govern the set of functions W,(x, y), where W, is the chance 
that » events occur in the domain (0, x; 0, y). The problem of solving them reduces to the 
solution of a sequence of first-order ordinary equations, but these are exact, so the solution is 
merely a matter of quadratures. The general solution is written, and a few simple illustrations 


are given, 





I. INTRODUCTION 


N spite of the frequent use of differential 

equations in special problems of probability 
theory, it appears that there is a need for explicit 
statement of the general equations which govern 
the probability of having m events in a given 
domain described in terms of variables x, y, etc. 
The purpose of this paper is to derive these equa- 
tions, and to give their solutions for some 
important cases. The problems in which such 
probabilities are encountered are exceedingly 
diversified and numerous. They range from 
astronomy to zoology, from studies of the dis- 
tribution of stars and nebulae to problems con- 
cerning the distribution of bacteria in culture 
dishes. The general features of such problems can 
be well understood from a simple example. 
Suppose we have a piece of cloth, of length x and 
width y; and let f(x, y)dxdy be the chance that 
a flaw occurs in the element dxdy. Then the 
chance that there are m flaws in the area xy may 
be called W,(x, y) and we desire a set of differ- 
ential equations from which we can find the 
functions Wo, W, etc.. 

In 1910 Professor Bateman! solved a problem 
of this type involving a time series. Considering 
the emission of alpha-particles from a radio- 
active source of constant strength, he assumed 
that the chance an alpha-particle will be pro- 
duced in time dt is fdt, where f is a constant, 
namely, the average number emitted per second. 
Bateman showed that the chance of getting m 


* Presented before the Southeastern Section, American 
Physical Society, April, 1943. 
1H. Bateman, Phil. Mag. 20, 698 (1910). 
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particles in time ¢ is 
W,= (ft)"e"'/n! (1) 


In 1935 Ruark and Devol? formulated the differ. 
ential equations for a more general case. Suppose 
we know f,(t)dt, the probability that an event 
will occur in dt when r have already occurred in 
time ¢. Then the chance W,(t) that n events will 
occur in the time interval (0, ¢) is governed by 
the equation: 


dW,,/dt=fn_1Wr-1—faWn; 


When 1 is 0, we omit the first term on the right. 

The fact that the chance of the occurrence of 
an event now depends on the number which have 
preceded it represents a considerable increase in 
scope. It is easily possible to go further, dealing 
with the situation in which f depends on the 
detailed arrangement of the previous events in 
time. Dependence of f on the occurrence of still 
other types of events can also be included, lead- 
ing to systems of coupled equations describing 
the chance of each type; but problems which 
require these refinements have not yet come to 
the writer’s attention. 

It has become customary to speak of the events 
as ‘‘random”’ when f is a constant, and we shall 
adhere to this nomenclature. 


n=1,2,-++. (2) 


II. EQUATIONS FOR PROBABILITIES WHICH 
DEPEND ON SEVERAL VARIABLES 
For brevity’s sake, we shall deal with two 
variables since the generalization to more than 
two is simple. 
2 A. E. Ruark and L. Devol, Phys. Rev. 49, 355 (1936). 
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1. Rectangular Coordinates 


Referring to the above illustration of’ the 
piece of cloth, let fadxdy be the chance that an 
event occurs in the element dxdy when n have 
occurred in the rectangle whose opposite corners 
have the coordinates 0, 0 and x, y. We seek the 
value of W,(x+dx, y). Holding y constant, we 
can get m events in the rectangle of base x+dx 
and height y, in two ways convenient for our 
purpose. Writing the sum of their probabilities, 


we have: 


W,(x, y) + (OW,,/dx)dx 


=Was{ f forty a+ wi -( f tay) as} 


A similar equation holds for dW,,/dy. Following 
a suggestion of Professor John W. Lasley, we 
shall let /“fdy denote an integral in which x is 
held constant. Writing 


y zx 
f tied. f fod = Yu (3) 
0 0 


we finally have 
aW,/dx = —XoWo; dW) /dy= — YoWo; 
OW,,/dx =X n1Wn1—XnWn) 


; n=1,2,etc. (4) 
dW,/dy= Y,-1W,-1- } .W,) 


The boundary conditions are: Wo=1 and all 
other W’s are zero when x=0 and y=0; all W’s 
are zero when x= and y is finite, etc. The 
notation in (3) is chosen for the sake of general- 
ization to more dimensions. In three dimensions, 


for example, 
z u 
X= ff faye. 
0 0 


2. Cylindrical Coordinates 


A second example will help to show how the 
method should be handled. For good practical 
reasons we have defined W,, to be the chance of 
n events in the particular rectangle mentioned. 
Of course, we could seek for equations governing 
the chance of having m events in any contour, in 
terms of any coordinates, but usually this quan- 
tity would not obey a manageable set of equa- 
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tions; to get tractable results, the contour should 
in general follow coordinate lines. Let us consider 
the case of axial symmetry in cylindrical coor- 
dinates. The reader may think of the problem of 
flaws in a length z of a cast metal bar of radius R, 
whose properties do not depend on the azimuth 
in a plane cross section. Let f,(r, z)2ardrdz be 
the chance of a flaw in the range drdz when n 
have already occurred in the range r, z. Then, 
if we let the area zr* be a new variable a, and 


write 
ff t05=A., ff t.00=2., (S) 
0 0 


we obtain, 
OW,,/da =A — —A aW, 


n=1,2,etc., (6) 
OW,,/d2= Zn-1 n—1— Z,W,, 


which are of the same form as (4). Here again we 
omit the first term on the right in case n=0. 


Ill. INTEGRALS OF THE EQUATIONS IN 
RECTANGULAR COORDINATES 


The general theory of simultaneous linear 
partial equations is old, but I have been unable 
to find a reference to this particular set. Equa- 
tions (3) and (4) have some interesting properties 
connected with the special character of the coef- 
ficients which are not obvious at first sight. They 
can always be integrated as follows: 


1. Calculation of W, 
We have 
0 log Wo/dx = — Xo, etc., 
so that 
Wo=exp (—So), (7) 
where 


dSo=X odx+ Yody. (8) 
Here dS» is a complete differential because of 
(3), and we obtain So by integrating it from the 
point (0, 0) to the point (x, y). 
2. Calculation of W,, 
We have 
0W,,/dx+0W)/dx= —X,W,, 


with a similar relation involving Y,W;. Thus 
X,W, and Y,W, are components of a gradient 
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the same is true of X2W2, Y2We:; and so on. 
Therefore, 


0(X,W,,)/dy=0(Y,W,)/dx, 
and noting that 


0X ,,/dy=dY,/dx, 
we get ° 


X,0W,,/dy= Y,0W,,/dx. (9) 
By Lagrange’s method it follows that 
W,=F,(Sz), (10) 


where S,=c is a solution of the auxiliary equa- 
tion, dx/Y,=dy/X,. Thus we have a practical 
method for finding how the variables x, y are 
contained in W,, for X,dx+Y,dy is always a 
perfect differential dS,,. 

After getting the S’s we need only to note 
that exp (S,) is an integrating factor for Eqs. (4). 
Multiplying each one by this factor, making use 
of (9), and adding, we get 


d(exp (S,)W,) =exp (Sn) Wad Sn-1 


and remembering the boundary conditions, 


zy 
W,=exp (-sy) f exp (Sn) Wa—-1dSn-1; 
0,0 


n=1,2,etc., (11) 


which can be generalized to any number of 
dimensions. For the case of one dimension, with 
constant values of the f’s, these are the equations 
governing the disintegrations in a radioactive 
series, for which Bateman! gave the explicit 
solution. There is no difficulty about writing out 
the explicit result of (11) when the f,’s are all 


ARTHUR E. RUARK 





and it follows from the equations for W2 that 





constants, but the results are lengthy, so we 
omit them. 


3. The Case of Independent Events 


In this case f does not depend on n. We can 
write, X, Y, Sin place of X,, Yn, Sn. All the W’s 
are functions of S and can be found as follows: 
Substituting (10) in (4), we get 


F,/0S/dx=X(Fr1— F,), 


but X =0S/dx, so F, obeys an equation of the 
Bateman form. The complete solution of the 


problem is 
W,=S"e-5/n! (12) 


This result can be restated in an instructive 
form. If S is considered as a new independent 
variable, the distribution is random in terms of S. 
On a given surface S=c, each W is constant. 
The surfaces S=c may be thought of as one set 
of coordinate surfaces in a curvilinear system, 
the other sets being chosen in any convenient 
fashion. 


4. Random Events 


When f is constant, we have S=fxy. From 
general considerations we see that in this case 
the value of W, associated with any area A is 


W,=(fA)re“4/n! (13) 


We refrain from further applications, simply 
noting that most of the useful theorems and 
formulas which have been developed in studies 
of fluctuations in one dimension can easily be 
generalized with the aid of the results presented 
here. 

I wish to thank Dr. Nathan Rosen for reading 


this paper. 
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The following problem is treated: Electrons enter the 
space between two infinite parallel planes with uniform 
velocity at right angles to the planes. The current is studied 
for all possible values (positive and negative) of the po- 
tential difference between the planes. The complete solution 
can be obtained if this boundary condition for the current 
is accepted: The number of electrons entering the discharge 
space must be equal to or smaller than a given number No 
per cm? per sec., and is for each potential as high as the 
potential permits. The problem depends on two dimension- 
less parameters, a reduced plate distance £ (involving the 
current) and no= Vo/Eo, where Vo is the impressed poten- 
tial difference and eZ,» the energy of the electrons. For each 
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value of £ the current is space-charge limited below a 
critical value of 9. The limitation, however, is not due to 
the appearance of a potential minimum low enough to stop 
the electrons, but is due to the fact that the solution does 
not exist unless the current is limited. The space-charge 
limited characteristic is a simple generalization of Child 
and Langmuir’s well-known formula [see Eqs. (28) or 
(35) ]. The deviations from this formula are very consider- 
able if Eo~ Vo (see Fig. 2). The behavior of the potential 
distribution is graphed for some typical cases of positive, 
zero, and negative potentials Vo. Currents can pass with 
zero or negative potentials only below critical values of 
the reduced plate distance. 





I. INTRODUCTION 


T has been understood, for almost 30 years, 
that the well-known space charge equation of 
Child and Langmuir represents a limiting law 
which is valid only when the initial velocities of 
the electrons approach zero. The influence of the 
initial velocities, first recognized by Schottky, 
has been studied theoretically in a very complete 
way under the assumption of a Maxwellian 
distribution. However, there remains a gap in the 
theory. The case where all electrons entering the 
discharge space have a uniform velocity has not 
been treated in a satisfactory way, not even in 
the case of parallel planes where the initial 
velocity is at right angles to the planes. 

This latter problem has been treated in an ap- 
proximate way by Langmuir and Compton,? but 
their solution is subject to the same criticism as 
the original Child-Langmuir solution, inasmuch 
as it leads to an infinite space-charge density (and 
velocity zero) at the point where the potential is 
minimum. The following analysis will show, 
furthermore, that one essential assumption in the 
derivation of Langmuir and Compton is too 
restricted. The problem itself has practical impli- 


1W. Schottky, Physik. Zeits. 15, 526, 624, 656, 872 
(1914); Ann. d. Physik 44, 1011 (1914); Verh. d. D. Phys. 
Ges. 16, 490 (1914). P. S. Epstein, Verh. d. D. Phys. Ges. 
21, 85 (1919). T. C. ig Phys. Rev. 17, 441 (1921); 22, 
445 (1923). I. Langmuir, Phys. Rev. 21, 419 (1923). E. L. E. 


Wheatcroft, J. I. E. E. 86, 473 (1940). 
* 1. Langmuir and K. T. Compton, Rev. Mod. Phys. 3, 
191 (1931), see p. 239. 


cations, and it seemed desirable to the author to 
attempt its complete solution. 

It may appear strange that the more complex 
problem of Maxwellian distribution has been 
solved in an entirely satisfactory way, whereas 
the simpler problem of uniform velocity remains 
open to discussion. However, there is an essential 
difference between the two cases. If all velocities 
are represented, e.g., in the Maxwellian distri- 
bution, each electron may be considered as 
subject to an impressed field treated as continu- 
ous and created by electrons of other velocities. 
Then there is no ambiguity as to how any one 
electron is going to behave in a field so impressed. 
This supposition has been clearly recognized and 
formulated by Epstein.* If, however, all electrons 
have the same velocity, they have to travel in a 
field created by electrons of the same velocity, 
and it remains an open question whether these 
electrons are able to create a minimum of po- 
tential which is low enough to stop them. It will 
be seen that this occurs exclusively in a limiting 
(unstable) case. In general, the minimum will not 
be low enough to stop electrons which have 
entered the discharge space. In spite of this fact 
the current will be space-charge limited (for 
sufficiently low values of the potential) because 
the solution of the problem does not exist for any 
given value of the plate distance unless the 
current is below a critical value. Thus the space- 


3 P. S. Epstein, reference 1, see p. 86. 
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charge limitation is caused, in this case, not by 
the fact that some of the electrons return to the 
cathode, but by the fact that only a limited 
number of them can get into the discharge space. 
The analysis of the problem is quite straight- 
forward. The gist of it is the determination of the 
limits of validity of the solution. If these are 
known, it is easy to determine the proper form of 
boundary condition for the current which will 
have to replace the artificial boundary condition 
[i.e., (€E/dx)o=0] of Child and Langmuir. 


Il. THE EQUATIONS AND THEIR SOLUTION 


We shall treat here exclusively the case where 
the electrodes are infinite parallel planes and 
where all electrons are supposed to enter the 
space between the two planes with the same 
velocity vo at right angles to the planes. Let x be 
measured in this normal direction and let the 
emitter plane, at x=0, be kept at the potential 
E=0, and the receiver plane, at x=d, at the po- 
tential E= V» (not necessarily positive). 

The equations of the problem are given by 
Poisson’s equation 

@’?E/dx? =4rne, (1) 


by the equation of continuity, and by the equa- 
tion of motion for the individual electron. We 
write the latter two equations in the integrated 
form. If m and v represent the number per cm‘ 
and the velocity of the electrons, respectively, we 


have 
evn =1 (2) 
and 


v=09(1+E/E»)}. (3) 


Here i represents the current density‘ and Eo the 
potential which is characteristic of the initial 


velocity 
eo = mvo"/ 2. (4) 


Substituting from (2) and (3) into (1) and 
integrating once, we obtain 


(dE/dx)? = (16miEo/vo)(1+E/Eo)'+Ci. (5) 
Here we introduce dimensionless variables 


n=E/Eo, §=(16ri/(voEo))*x, (6) 


* We count, in the usual way, ¢ as positive if an electron 
current is flowing in the direction of the positive x axis. 


and obtain 
dn/d§= +((1+n)!+a,)', (7) 
where a; is a dimensionless constant of integra- 
tion replacing Ci. 
Finally, we set 


u=(1+n)! (8) 
and find, after a second integration, 
(4/3) (u—2a1)(u+ai)'= +é+a2, (9) 
with a2 as second constant of integration. 

It follows from (7) that dn/dt, or dE/dx, may 
be positive or negative. Limiting ourselves at 
first to the case Vp >0, we have to distinguish two 
possibilities. Either dn/d§>0, and 7 will increase 
in a monotonic way, or dn/dé <0 at first, and 7 will 
present a minimum. 


First Cdse, Monotonic Increase 
Let =» correspond to the distance d and let 
uo=(1+0)'= (1+ Vo/Eo), (10) 


then the conditions imposed on the potential 
demand u=1 for ¢=0, and u=wupo for &=§£,, 
Hence, from (9) (with the plus sign) 


(4/3)(1—2a1)(1+a1)'=ae (11) 
and 
(4/3) (uo — 2a1)(uo+ ay)! =Fo+ az. (12) 
These two equations determine the constants a, 
and a provided 7 is known. The field distribution 
then follows from (9) (with the plus sign). 


Second Case, Negative Potential Minimum 


Let 6 be the value of — for which the minimum 
occurs. The value of 6 is obtained by integrating 
(7) from §=0 to §=6 and bearing in mind that 
dn/dt = —(u+a;)'=0 for §=6. Thus it is found 
that 

56= (4/3)(1—2a:)(1+.a;)!. (13) 

The solution now becomes 

(4/3)(u—2a1)(u+a:)'=—E+a3, O<E<S, (14) 


and 
(4/3)(u—2e1)(u+ai)'=Et+as, §26, (15) 
where a; and a, are new constants of integration. 


They are determined from the boundary con- 
ditions for E, and from the condition of continuity 





at 
an 


W 
bu 


th 


gr 
of 
sic 
is 
int 
of 
ter 


CURRENTS CARRIED 


BY ELECTRONS 




















“, 


Fic. 1. The curves represent the dependence of the reduced plate distance 
£>) on the constant a; for various values of the external potential (wo =(1-+-n0)! 
= (1+ Vo/Eo)*=const.). Solutions of the problem for prescribed values of £ 
and mo are given by the intersection of the horizontal £)=const. with the 


curve #9=const. 


at §&=6. We find 
(16) 


az=6, ag=—, 


and, for the determination of a, 
(4/3) (uo —2a1) (uot+ar)*=fo— 6. (17) 


When the constants are known, the field distri- 
bution follows from (14) and (15). 


Il. THE BOUNDARY CONDITION FOR THE 
CURRENT 


It will be noticed that the Eqs. (11) and (12) in 
the first case, and Eqs. (13), (16), and (17) in-the 
second case, determine the constants of inte- 
gration as functions of 9 and &9. Thus the solution 
of the problem depends on these two dimen- 
sionless parameters only. However, the solution 
is not complete since £) involves the knowledge 
of the current density 7. Now 7 is a constant of 
integration itself, originating in the integration 
of the equation of continuity, and has to be de- 
termined by a boundary condition. 

We shall assume conditions under which a true 


saturation current exists. Our problem does not 
correspond to thermionic emission, but rather to 
the case where electrons of uniform velocity are 
shot into the discharge space. Therefore, it is not 
necessary to include the possibility of a Schottky 
effect. Preliminarily we assume that not only the 
velocity of the electrons is given, but also their 
number. Let No=movo signify the number of elec- 
trons which enter the discharge space per cm? per 
sec. Then the saturation current will be 


(18) 


1=19 = ENV. 


With this boundary condition the current will 
be constant, i.e., independent of the potential Vo, 
within the limits of validity of our solutions. 

It should be noticed that Eq. (9) will give no 
physical meaning unless the variable wu is real. 
Consequently, from (8), 7 is subject to the con- 
dition »> —1, and the potential minimum cannot 
become as low as qmin= —1, and eventually stop 
the electrons, except in a limiting case. We are 
going to discuss this limiting case later. 
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It remains to be seen within what limits our 
solutions are valid. The practical situation will 
usually be that 70 and & are given and that the 
constants a have to be determined. However, the 
discussion becomes simpler if we assume that, 
besides yo, a: is given and that, consequently, £ 
and the other a’s have to be determined. The 
equations available for the determination of £ 
are (11) and (12) in the first case and (13) and 
(17) in the second case. 

It will be seen, from (11), that the constant ae 
in the first case is identical in form with 6 in the 
second case, though it does not have the same 
physical meaning. Therefore, we can write (12) 
and (17) in the joint form 


£ot6= (4/3) (uo — 2a) (uo+a)!, (19) 


where the plus sign corresponds to the first case. 
We shall refer to the two cases as “branch |” and 
“branch II” of the solution. 

Because of the last factor in (13), the constant 
a, is subject to the limitation 


a,>-—1, for branch I and II. (20) 


Furthermore, on branch II, a is limited to 
negative values as otherwise dyn/dé could not 
become zero, according to (7). Hence 


0>a,>—1, for branch II. (21) 


If a; is known, the value of the reduced potential 
n at the minimum nmin is given, from (7), as 


Nmin = — (1 — a’). (22) 


In Fig. 1 the two top curves represent £ as 
functions of a; for the two values up = 2 and up=3 
which correspond to the positive values yo = 3 and 
no=8. The curves are calculated from (13) and 
(19). Branch I is shown by a full line and branch 
II by a broken line. Branch II has been calculated 
for the interval (21) only as there is no solution 
of our problem beyond that interval. A point on 
branch I represents a possible solution without a 
potential minimum, a point on branch II a 
possible solution with a minimum [given by 
(22)]. A point where the two branches join 
represents the limiting case where the potential 
curve starts with a horizontal tangent at the 
cathode. : 

It is easy to show that the curves for all values 
of uo>1, i.e., 70>0, represent the same general 


character as that exhibited by our two examples: 
Both branches begin at the value 


£o= (4/3) (wo+2)(uo—1)! (23) 


with a vertical tangent. Branch I decreases in a 
monotonic way and approaches £)=0 asymptot- 
ically. Branch II has a unique maximum at 


a= — (uo/(uo+1)), (24) 
which is given by 

Emax = (4/3)(1+u0)!, (25) 
[from substitution of (24) into (19) with the 


lower sign ]. 

It is the existence of this maximum which is 
decisive for space-charge limitation. To show this 
we now return to the usual state of affairs where 
no and & are prescribed. Then the possible solu- 
tions will be given by the intersections of the 
curve %=const. with the horizontal £)=const. 
Evidently that leads to three possibilities— 
namely, that there is one solution, that there are 
two solutions, or that there is no solution at all. 

We discuss the last of these possibilities first: 
it will arise when the given value of £o is larger 
than the maximum value (25) corresponding to 
the given value of uo. Thus, under our preliminary 
boundary condition, there will be saturation 
current for all values of no larger than the critical 
value mo* given by 


fo = (4/3)(1+u0*)!, (26) 


and for smaller values there will be no solution. 
This is a result which cannot be admitted from 
the physical point of view, and the boundary 
condition must be modified. 

From the definition of the variable & [Eq. (6) ], 
it follows that £ can be lowered by decreasing the 
value of 7 involved and that means by decreasing 
No Since Vp is supposed to be constant. Thus, when 
the critical value of uo* (26) is reached (by 
decreasing the potential), a solution for lower 
values of » becomes possible if mo is suitably 
reduced. Therefore, we relax our boundary con- 
dition to the following statement : The number of 
electrons entering the discharge space has to be 
equal to or smaller than the given number No=nm% 
per cm? per sec. and is for each potential as 
high as the potential permits. The last clause has 
to be added because otherwise the solution 
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Fic. 2. The curves represent the dependence of the re- 
duced current j =9xid*?/(2e/m)* on the external potential 
Vo (in volts) for different values of the initial velocity 
(Eo=e volts). 


would not be determined. It would be possible to 
obtain solutions for all values lower than ip 
before the critical mo* is reached. 

With this boundary condition the current will 
decrease, from the critical value no* downward in 
such a way that (26) remains fulfilled since this 
relation gives the highest value of 7 which can 
pass for given values of d and 7. It should, 
however, be stressed that the limitation is not 
due to the appearance of a minimum which makes 
the electrons (or some of them) return; for the 
minimum of potential, as given by (22), is less 
low than nmin=—1 for all relevant values of a; 
given by (24). Therefore, the space-charge limi- 
tation, as stated in the introduction, is due in our 
case to the fact that, for a given potential, 
electrons of given velocity can enter the dis- 
charge space only in a number limited for each 
plate distance. 


IV. DISCUSSION OF THE SOLUTION 


Our Eq. (26) contains the space-charge limited 
branch of the characteristic. If we bear in mind 
the definition of u» (10) and & (6) and introduce 
the “reduced current density j’’ by 


j =9wid?/(2e/m)}, (27) 


CARRIED BY ELECTRONS 


the relation (26) can be written in the form 
j= Voll (1+ Eo/Vo0)'+(Eo/Vo)*}. 


This is evidently the generalization of the 
Child-Langmuir law to the case of electrons with 
uniform initial velocities, and reduces to it in the 
limit Eo/Vo—0. It need hardly be pointed out 
that our solution does not lead to infinite values 
of the space-charge density, neither at the 
cathode nor at the potential minimum. Further- 
more, it is automatically limited to the values 
i <1». As soon as the value 1 =i» has been reached 
by increasing the potential the space-charge 
limited characteristic, (26) or (28), has to be 
continued by the horizontal saturation branch. 

For small values of Eo/V> we obtain 


j= Voll1+3(Eo/Vo)'+- ++), 
and for large values of Eo/ Vo 
j =8Eo!(1 + (3/4)(Vo/Eo) +--+ J. 


Thus, for Ve=0 the current reduces to the 
constant value jo>=8E,!. However, the charac- 
teristic can be continued into the domain of 
negative V ’s (see Section VI). 

In Fig. 2 we have represented the reduced 
current density j as a function of Vo (in volts) for 
various values of Eo (also in volts). All curves 
ultimately (for Voe—~) approach the curve with 
E,=0, which is the Child-Langmuir law. It will 
be seen that, for Eo comparable with Vo, there 
are very large deviations from the $-power law, 
as was known from the study of Maxwellian 
distribution. 

We now revert to the case where the given 
values of & (with t=i9) and mo lead to two 
solutions of the problem. As an example, the 
horizontal &=10 in our Fig. 1 intersects the 
curve u%=3 in the two points B and D corre- 
sponding to a,;=—0.39 and a,;=-—0.975, re- 
spectively. The two solutions which are thus 
possible may be both on branch II, or one on 
branch I and the other on branch II; the latter 
case occurs only if w»>1.81. 

The two solutions belong to different values of 
a, and, therefore [according to (9), (14), and 
(15) ], have different potential distributions. It is 
easy to show that the solution with the smaller 
negative value has the lower minimum (if both 
solutions are on branch II) and that it leads to 
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Fic. 3. The curves represent the potential distribution 
in the case £9 = 10 for typical values of the ratio no= Vo/Eo. 
At no=m the potential minimum begins to appear and at 
no=no* the current becomes space-charge limited. The 
broken line corresponds to the unstable form of discharge. 


lower values of the potential for any value of ~ 
except £=0 and = & . It would require work to 
change the discharge with the higher potential 
curve to that with the lower potential curve, and 
consequently the former discharge is the stable 
one. Thus we arrive at the result that points to 
the right of the maximum of branch I], in Fig. 1, 
represent unstable forms of discharge. 

We obtain the following description of the be- 
havior of a discharge for a given value of & as the 
potential is decreased. For very large values of uo 
there will be saturation. The intersection between 
the curves £)=const. and u»=const. will be at 
large values of a, i.e., on branch I, and the 
potential will be almost linear (see next section). 
As up is decreased, the intersection remains on 
branch I until a; has reduced to the value —1. 
This will happen at the value u, given by® 


£o = (4/3) (ui +2)(u,—1)!. (31) 


At this point the potential curve has a horizontal 
tangent at §=0. As the potential is further re- 
duced the intersection will be on branch I], i.e., 
a maximum begins to make its appearance. The 
minimum gets lower and 6 gets larger (see the 


5 Equation (31) represents a cubic for u,; which is easy 
to solve. It has one real root only. 
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curve for 6 on Fig. 1) until the point is reached 
where the horizontal )>=const. is tangent to a 
curve %»=const. This will happen at the critica] 
value of up given by (26). Until this point is 
reached the current is saturated. If now the 
potential is further reduced, the representative 
point will follow downward the line joining the 
maxima of the various uo curves. 

We still have to say a few words regarding the 
limiting case where the potential at the minimum 
becomes as low aS qmin=—1, i.€., Vmin= —Ep. 
This happens exclusively for the value a,=0 
[see (22) ]. The event is, therefore, characterized 
by points of branch II located on the ordinate 
a,=0. We have just seen that such points repre- 
sent unstable forms of discharge. 

Apart from this fact there are strong reasons 
which speak against the physical admissibility of 
this limiting case. For any given value of uo the 
case a; =0 (on branch I1) occurs for a finite value 
of &. Therefore, we must draw the conclusion 
that the current density 7 would be finite, as 
given by this value of &). Now, the velocity of the 
electrons naturally becomes zero at a point where 
n=-—1 [see (3)], and therefore the density n 
would have to be infinite. We are of the opinion 
that no solution should be admitted which makes 
a physically measurable quantity infinite. We, 
therefore, conclude that the limiting case a; =0 
(on branch I1) should be excluded, even as a 
possible form of unstable discharge. 


V. THE POTENTIAL CURVE 


After the constants a have been determined as 
indicated above, the potential distribution be- 
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Fic. 4. The curves represent the potential distribution 
in the case Vo=0 for two values of the reduced plate 
distance. The broken line corresponds to the unstable form 
of discharge. 
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tween the plates becomes calculable from (9) or 
(14), (15) for the two branches. We first discuss 
the case of very high potentials (79>1). 

We obtain the value of a, corresponding to 
this case by developing (19), with the plus sign, 
and (13) in powers of uo/a. This yields, to a first 
approximation, 


a} wget (32) 


Substituting this value into (9) and (11) we find 
as limiting law 


f= (£0/no)n. (33) 


Thus, in the limit yo, the potential distri- 
bution becomes linear, as has always been as- 
sumed and as is evident from the fact that all 
electrons are swept out of the discharge space 
with increasing velocity as the potential is 
increased. 

It is not worth while to establish approximate 
formulae for any other cases as the exact expres- 
sions are easy to handle. In Fig. 3 we have 
graphed some potential curves for the case® 
t)=10. We have drawn the curves for the value 
n=m=8.79 where the minimum begins to 
appear, for no=8, and for the critical value 
no=no* =7.02 where the space-charge limitation 
sets in. For no=8 there are two solutions as we 
have seen above. The unstable form of discharge 
is given by a broken line. 

VI. ZERO POTENTIAL AND NEGATIVE 
POTENTIALS 

In our discussion we have limited ourselves so 
far to positive values of V» though our solutions 
are valid for all values of 79> —1. It is of interest 
to discuss zero and negative values since the 
equations developed under the assumption of 
Maxwellian distribution are not applicable to 
these cases.’ 

If Vo=0, or u)=1, the branch I of the curve 
uo= 1 degenerates into the horizontal axis and the 
branch II takes the form indicated in Fig. 1. It 
hasits maximum £max = (4/3)2!=3.771ata,;= —}. 
Therefore, currents can pass with Vo>=0 only if 
fo<3.771. For the interval. 2.667 < )<3.771 
there are two solutions, one of them unstable. 

* The case {)>=10 does not correspond to any practical 
case of thermionic emission. Much larger values of & would 
have to be chosen (order of 10° to 10‘). However, the 
minima are much more marked with comparatively low 


values of no and £p. 
71. Langmuir, reference 1, see p. 424. 
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In Fig. 4 we have drawn the potential distri- 
bution for Vo=0 and for two values of £o, 
namely, the limiting case )=3.771, and the 
intermediate case §)=3.33. The unstable distri- 
bution is given by a broken line. 

If, finally, Vo becomes negative we have again 
to distinguish two cases, namely, the case of a 
potential decreasing in a monotonic way : branch 
I, and the case of a negative potential minimum: 
branch II. By calculations which are exactly like 





j— @& 1 LS 2 
rh T T —" 
‘ (1) } 
\ y= -0.75 
i; ON (2) 
‘\ 
-05 a \ 
\. (") 
i (2) Pw 
at ie SO 





Fic. 5. The curves represent the potential distribution 
in the case of a negative potential Vo= —(3/4)E» for the 
value £=1.886. The broken line corresponds to the 
unstable form of discharge. 


those following Eq. (9), we find for the determi- 
nation of a, on branch I (negative potentials) 


5 — fo = (4/3) (uo — 2a) (uo +a;)!, (34) 


where 6 retains its definition (13). 

For branch II our previous formulae remain 
valid without any change. Thus a can be de- 
termined from (19) with the lower sign. Also the 
maximum, as defined by (24), (25), retains its 
significance and the current begins to decrease on 
the space-charge limited curve as soon as the 
critical value is reached. Therefore, the charac- 
teristic (28) remains valid though it has to be 
written now in the form 


J=((Eo+ Vo)*+Eot ¥. (35) 
The reduced current decreases according to 
this equation down to the value Vo= — Eo giving 


j=E£E,!. For still lower values of Vo, no solution 
exists and no current, therefore, can flow. We 
have continued the characteristics of Fig. 2 down 
to their limit. The tangents are vertical at the 
limits. 

As for the limits of validity of the solutions, the 
argument remains similar as for positive po- 
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tentials. It follows from (34) and (19) (lower 
sign) that solutions exist only for 


a,> =~ Me. (36) 


Both branches of the curves u»=const. begin 
with the value a;= —wo at the value é.=6 and 
with vertical tangent. They form together a loop 
which contracts around the point G on Fig. 1 
(a,=0, &=4/3) as uo decreases from up=1 to 
uo=0, and for uw»=0 there remains only the 
branch I. In Fig. 1 the curves for u)=} and 
uo =0 have been added. 

It is now easy to discuss the behavior of the 
solution. Three intervals have to be distinguished : 

(a) 0<&<4/3. The solution starts, with 
Vo=0, on branch II until, with decreasing po- 
tentials, the value ).=6 is reached. Then the 
solution changes to branch I (no minimum) and 
ends on the curve %#)=0. The current has the 
saturation value throughout. 

(b) 4/3<&<1.886. The upper limit of this 
interval is set by the maximum of 6 which occurs 
for a;= —0.5. The solution begins, with Vo=0, 
on branch II, changes to branch I when the 6 
curve is intersected for the first time, changes 
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back to branch II when the 6 curve is met for the 
second time, and finally becomes space-charge 
limited. 

(c) 1.886< &<3.771. In this interval the soly. 
tion is on branch II throughout. The current js 
saturated at first and becomes space-charge 
limited as soon as the critical value (26) is 
reached. All space-charge limited solutions [in 
cases (b) and (c) ] end at the point G of Fig. 1, 
i.e., at £5 =4/3 or j7=Eo!. 

‘In Fig. 5 we have represented the potential 
distribution in the case of a negative potential 
Vo. We have taken the value uo =} corresponding 
to Vo= —(3/4)Eo, and have chosen £9 on the § 
curve; i.e., §)= (4/3)2'= 1.886. Whenever, in the 
intervals (a) and (b), the limiting case )=8 is 
realized, the solution of branch I coincides with 
the (stable) solution of branch II. This case, 
therefore, corresponds to the case (dn/d£)o=0 for 
positive potentials, but the horizontal tangent is 
now at the receiver plate. 

The two solutions on branch II correspond to 
the values a; = — } and a, = —0.025, respectively, 
Both potential distributions are represented in 
Fig. 5, the unstable one by a broken line. 
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The vector wave equation in prolate spheroidal coordinates £, 9, ¢ is set up, the variables are 


separated, and the characteristic values (eigenvalues) and characteristic functions (eigenfunc- 
tions) of the resulting ordinary differential equations are obtained in series which converge 
rapidly in the neighborhood of resonance for spheroids of eccentricity near to unity. The 
coefficients of the two independent primitives in the linear combinations representing diverging 
and converging waves at infinity are calculated. The zeros of certain of the characteristic values 


are investigated. 


N a previous paper! both the free oscillations of a perfectly conducting prolate spheroid and the 
oscillations forced by a plane wave with the electric field parallel to the long axis of the spheroid 





were discussed. The forced oscillations were treated by an approximate method valid only for very 
eccentric spheroids. For the exact investigation of these oscillations the complete solutions of the 
vector wave equation in prolate spheroidal coordinates are needed. The object of the present paper is 


1. Page and N. I. Adams, Jr., Phys. Rev. 53, 819 (1938). This paper will be referred to as I and the present paper as II. 
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to obtain these solutions in the form of series which converge rapidly in the neighborhood of resonance. 
In a following paper these solutions will be applied to the problem of the antenna. 






1. FIELD EQUATIONS 





If the time factor is taken as e~*“’, where r is the time, the field equations for simple harmonic waves 
of angular frequency w in a medium of permittivity « and permeability u are 








VXE=i-H, VXH=-i—E, (II-1) 
c c 


in Heaviside-Lorentz symmetrical units, leading to the wave equations 









2 
re ~ H. (II-2) 









w Ku 
VXVXE=——E, VxXVxXH= 
c c 









In the right-handed orthogonal prolate spheroidal coordinates £, n, ¢ defined in I, either of these 
vector wave equations yields the three scalar equations 



















































{ot “|- ee ta OF + _ oF, (II-3) 
Onin? —& EJ Anln?—£ On J] (1—£€*)(n*—1) OG? °—1 OEdG 
1 oF, 1 aF, “alee . tile aFy adh 1. 
a dnd (1-2) (92-1) 0g? dtLn?—& dE n?—& On J 
7 = 
| "| ~i a +2 = =<¢? os Fy, (II-5) 
aes 064 ?—1 0& 1—£€ On? ner ae rs (1—£*)(n?—1) 
where a on 
¢= { (n?— &) (1 —&*) yl’ 
.= =» £8 a» $538 . 
{ (n? — &*) (n?—1)} a 
— _ 4 ° 
={(1—#)(n?—1)} ! 
and f asf 
c= (ula, 
Cc A 






f being the semi-interfocal distance, and \ the wave-length. 
Evidently the solutions are of the form 





nm, r=G I, F.=G,| om . 
—sin md 


=<, 


" 
cos m@ cos md 





where m is a positive integer, and G;, G,, Gs are functions of £ and 7 only, satisfying the equations 






























rs] 1— 2 = 2 
“| =~ |-s/ = + m a. we (11-6) 
Onin?—£? OE J Anln?—£ On | (1—£€*)(n?—1) n’—1 0€ 
m aG 2 af 1—£ 8G.) af n?—1 8G) 

-—_—_*+_~ a-s| its: —|=8Gy (II-7) 
1—£ dn (1—£*)(n?—1) OtLn?—£& OE] AELn?—E£? On J 

















a G 2 2 2_ ¢2 
| mG: - 1 9G, 1 Se, Of mG, _. at G.. (IL-8) 
dL (1—&)(n?—1) J n?—-1 OE) =1—E* On? = nL (1—£*)(m?—-1)] (1 —€*)(n?—1) 
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If we differentiate (11-6) with respect to £, (II-7) with respect to n, and combine the two equations go 


obtained with (II-8), we get 
ve Fite m 4 m i hr 
at an =L1i—g 2-1) ” 7 














which is just the solenoidal condition. 
Eliminating G, from (II-6) and (II-7), we find 

































~ 3 a 2E (AG: AG), 2n [Ge OG.) _ wil 
— ir” COi—ae el eRe as = 
a ‘. 28 [8G. 8Gi] 29 [8G, aGe) 

+ fee. 4NGI— oa wane FM | SY (11-11) 
Li—g* 9?—1) 1—EL OE On n?—1L On 8 € J 























where S and WN are the operators 





+e2(1—£), 






s=(1-)—- 
=(- 







N=(n? *(n?—1). 








Equations (II-10) and (II-11) admit two pairs of solutions in which the variables are separable. The 
first pair is obtained by putting G;= &G, G,= —nG and the second pair by putting G;=G, G,= — &G. 
In both cases G satisfies the equation 





2 


{i | m? c-¢e G42 (a? v|- m? 
ri) ‘ of] 1-# "= On ' dn} 7°?-1 


So, if we put G=u(£)v(n), we are left with the two ordinary differential equations 


0G =0. (II-12) 















il" -#F/-; = wtaute(1—)u=0, —1Sé<1, (II-13) 






dv m? 
No? -12]- av+e%(n?—1)v=0, 1<n<~, (1-14) 
dn dnj 1*—1 





which differ only in the range of the independent variable. 

Finally, after G; and G, have been obtained, Gy may be determined from (II-9). 

Evidently the method fails when m=0. In this case, however, F;, F,, Fy are not functions of 
@, and (II-5) reduces to an equation in F, alone in which the variables are separable. If we put 
F,= {(1—&*)(n?—1)}*u()v(m), the differential equations satisfied by u and v are found to be just 
(II-13) and (II-14) with m=1. The functions F; and F, corresponding to this solution for F, are both 
zero since the field component proportional to F, is itself solenoidal. A second solution for the case 
where m=0 is obtained by taking the curl of the field intensity corresponding to the first solution. 

For a given m, solutions of (II-13) and (11-14) exist only for a discrete set of values of the constant 
of separation a, corresponding to the various harmonics possible. In fact, we conclude from the 
similarity of these differential equations to that for the associated Legendrian functions that a isa 
function of two positive integral indices /] and m such that m <1. The solutions of (I]-13) and (II-14) 
for a given ajm we shall designate by uim(£) and vim(n). Evidently each group of functions corresponding 
to a specified m constitutes an orthogonal set. If, then, we put 


sin md sin md Es - cos mod 
faze ofl, | fame A 1 fReE | oe 



























Aim ' 
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with omission of the time factor, we have the following pairs of solutions, distinguished by single and 
double accents: For m=0, , 


. 1 
=10 =———_ Ty) un zret — en) 


> 
LK 
= sta é*)!un(€) Jon(n), 


— 
=0=0, | 


Hin=0, | (II-16) 


:' - 19 = uni(€)vn(n) ; 
or m™ é 


LU 


~"i-oe-ore 








n 
{ (n? — 1) (?— &*)}! 
4 eT onl) mnt" | 
a ” 
"(1 8)? — 8} 
n= ——_— Tae ttin()Pm()., (11-18) 
a” = {(1—&)(n?—1)} 4 dui 
m—#) de 
2. SOLUTION OF THE EQUATION FOR u;,,(é) 
The form of Eq. (II-13) for uim() suggests a solution in the associated Legendrian functions P;»,(£). 
Hence we put 
tim =1(L+1) + 2am term +E aim + + *, 
Uim = Pim(§) +€ ~ Ax’ Pem(£) +e ~. an” Pim(€) + ~. Oy” Pim(E) +++, 


Uim(€)Vim(n), (11-17) 





Hin = — 


heed 





Uim(£)Vim(n) ’ 














2——Vim(n) — gual)" | 


substitute in the differential equation,’ and equate to zero the coefficient of each separate power of ? 
after getting rid of the factor 1 — £* in the last term by means of the recurrence formula 


(1—£) Pin = 0 (—1—2) 0 (—1—1) [Pim — Pina, m JA+TU—1)T(D [Pim — Pie], 
where I'(x) =(x-+m)/(2x+1). 


In this way we find the following formulas for the first four terms following /(/+1) in the series 
for the characteristic value (eigenvalue) aim: 
L(i+1 2—] 
_ td+1)+m , (11-19) 
(2/—1)(2/+3) 
» H(-1)HQ) HU+1)H(4+2 
© (¢-1)H@) H¢+1)H( ) (11-20) 


2(21—1) 2(21+3) 


* This is a much simpler and easier method of solution than that given in I. 
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TABLE I. Characteristic values (eigenvalues). 

















+ 4 8 124 
a1, =2—-&— f é® 


+ - € 
5 53-7 3+ 58-7 58-73-11 


4 a 8 5420 
a2 =6—-é— P e+ --- =6—0.571,429e2 —0.003,887 4 +0.000,014«°+0.000,001 «8+ ---, 


4 , 
3-78 3-78-11. 3% 77-11-13 
8 152 115,568 34,094,936 
an= 12— + — «6 Pa 
3-5. 353-11. 37-5®-11-13  3-58-113-13 
40 7064 462,736 


a4, = 20 —-—— f+ 


é— € 
7-11 73-113-13 7§-115-13 









8+ --++ =2—0.800,000e —0.004,571 et +-0.000, 1226 —0.000,002«8+ ---, 
























+++ =12—0.533,3332+0.001,365¢4 
— 0.000, 118¢«*+-0.000,002 «8+... 










+-++=20—0.519,481e2+0.001, 190 —0.000,013«8+---, 































20 1108 301,160 
as, =30———é e— +--+ =30—0.512,8212+0.000,889«! —0.000,003 «8+ ---, 

3-13 34-7 -138 37-7-135-17 

6 2 : 44 ~~ 6 = 84 6—0.857,143e2—0.001,944«4+ 0.000,036«' —0.000,001 €® 
a22=6—-é— € = +++ =6—0.857,1432—0.001, .000,036«* — 0.000, tee, 
et ark re See ATT 7 “ ; sl 
2 4 62 

a32=12—-é— ‘ 6 +--+ =12—0.666,667 2 —0.002,245'+0.000,013«'+0.000,000«8+ - -- 














+ + 
3 #11 37-11-13. 3. 113-13 










































m____ 4m*—1_—f AG-1)H@M —_-H(+1)H+2) 
_ aI na Gorn Geran - 
av _HC- HO 4(4m?—1)? " Cin H(J—3)H(/—2) 
"~~ 2(21—1) | (27—5)2(27—1)4(27+3)2 2(27—-1) 2-4(27—1) (22-3) 
_HG+1)H0+2)f 4(4m?—1)? Pad aan as 11-29 
2(21+3)  [(2i—1)2(2/+3)*(27+7)2 2(2/+3) 2-4(2/+3)(2/+5)] 






where H(x) =I'(x)I'(—x) = (x?—m?*)/(4x?—1). 
In fact, we can obtain by this method an equation for the entire characteristic value aim. Putting 





; b= (I+s)(I-+s+1) -e[P(—1—s—2)F(—1—s—1) +°(0+-s—1) P+) ]—am, 
= _H(+s)H0+s+1) 
a Re Roi1 








’ 







this equation, involving two continued fractions, is 








ex_1 ex) 
+ = 1, (11-23) 


e’xX_3 €°X3 






1—ex_;/1— 1—e*x;/1— 









However, the first five terms of aim are sufficient for most purposes, and they are more easily calcu. 
lated from (II-19) to (II-22) than from (II-23), the apparently simple form of which is deceptive. 
We list in Table I characteristic values for m=1 and /=1, 2, 3, 4, 5, and for m=2 and /=2, 3. 
When the characteristic functions in “jm(£) are expanded as a series in the associated Legendrian 
functions Pim(), the coefficients are found to be much more complicated than when they are ex 
panded in a simple power series in s=(1—£*)}. So the latter expansion has been adopted. We have 
to consider separately the cases ]—m even and !—m odd. In both cases it is found convenient to put 
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TaBLeE II. Characteristic functions (eigenfunctions). 
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2m+1 dz 
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> » CpS”, the recurrence formula for the coefficients is 
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TABLE II.—Continued. 
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For /—m odd we put uim(t) = §s"Z, getting 


@Z 2m+1dZ dZ 
(i~*}t Bo OTe FG ETFS tie Hoe a6, (11-25) 
5 5 5 











which gives, when we put Z=)>_, c,s”, the recurrence formula 


b(p+2m)cy— {(b—l-+m—1)(p+1+m) — Bim} Cp2+e*cp_«=0, 







representing a series that starts with Co. 
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TABLE II.—Concluded. 
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The characteristic functions u;»,(£) obtained from (11-24) and (II-25) are listed in Table II together 
with the functions v;»(), the arbitrary coefficient co in the tabulated functions u;,(£) being taken as 





unity. 





3. SOLUTION OF THE EQUATION FOR vin(n) 





Since we are interested primarily in conducting spheroids of eccentricity very close to unity 
(antenna), the objective in solving Eq. (II-14) for vim(m) is to obtain series solutions representing 
diverging waves at »= © which converge rapidly in the neighborhood of resonance for » nearly equal 
to unity. 

By changing the independent variable from 7 to p= en the differential equation can be thrown into 
such a form that the zero-order terms are satisfied by the wave functions 


_ (-— yp - < |e. 
p dp 


Therefore, a solution can be constructed in the form of a series of these functions, which meets the 
first requirement stated above. Unfortunately, this series converges so slowly in the neighborhood of 
resonance for 7 nearly equal to unity as to be quite useless for the antenna problem. 

Another possibility is to construct two independent solutions in series of associated Legendrian 
functions of the first and second kinds, respectively. The coefficients, however, turn out to be com- 
plicated, and difficult to calculate in the second case, and there remains the problem of finding what 
linear combination of the two independent primitives represents a diverging wave at infinity. 

Finally, it was concluded that the most satisfactory method was that employed in I. First, we 
obtain two independent primitives in the independent variable* :=(n?— 1)! which goes to zero as 
goes to unity, and then determine the coefficients of the desired linear combination by comparison 
with the solution in p= et which represents a diverging wave at infinity. As in the case of the equation 
for Uim(€) we have to consider separately the cases /—m even and 1—m odd. 

For /—m even we put Uim(n) =?"y, getting 


d*y 2m+1dy 
(144 yo EES + 20+ 1) — Im) (l4+-m+1)y—€(Bim—#)y=0. (11-26) 

















The first primitive for small ¢ is the power series y:= >> py apt” v ixere 


b(p+2m)ay+ {(p—1+-m—2)(p+1+m—1) —€Bim}ap-2+e*ap_4=0. 


* Note that ¢ used in this paper is the square root of the ¢ in I. 
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A second primitive is obtained by putting 





ate 
*(+e)'— 


1 


y2=391 lo pt +eyie, 






where x satisfies the equation 


















m+1d ‘x ld 
(1+? y+ +2042) —(I—m—1)(-+-m-+2)x—¢ (Bim —t*)x = = 4 my | (11-27) 
For /—m odd we put vim(n).= nt" Y, getting 
d*?Y 2m+idY dY 
(+e yee sy + 2(m-+2)-—— (I—m—1)(L+-m+2)¥—E(Bin—2)V=0. (I-28) 





The power series Y= >°, a,t” supplies the first primitive for small ¢, where 


b(p+2m)a,+ | (p—1+m—1)(p+1+m) —eBim}ayp2+eay_4=0, 












and the second primitive is 


(++i 1 
(i+e)'-1  (1+2)1 








Y,= 3 Y; log 









where X satisfies the equation 












d?X Ph dahon’ dX 
(1+2?) an ay FAQ + I — C— me) C++ IX —e* gy 


m 


To obtain the solution for large ¢ for the case 1—m even we change the independent variable in 


(11-26) to p=et and put 
1 m*—1 
a= 
2 2? 
















y -_ pt D (pe, 






getting 





dv 





(11-30) 





ay dv ay 
22°— — aim V = e? 


22@—_— 
dz? dz 














where s=ip. The solution is of the form 






W= > ¢,5° 


p=0 







where 


2pcy = {(p+I)(p—1--1) —e(Bim +1) }epr+e*(2p— 3)cp-2—e | (b—2)?— mf ep», 





After calculating the coefficients c, we write the solution in the form 





even odd 


y=- YY (—e)2C, t-te 44 YY (—e2)(-V2C, pomt ni 
E€ an n 









for comparison with the primitives for ¢ small, where 
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II 





For 1—m odd we make the same change in independent variable in (II-28) and put 


Y= p-™+2) (pet, 


d?y dv a?*y 3\dv 3 m*—4 
2°—_+ 22° — am V =e s+ (2-- _—" i-—-— )wl (11-31) 
2 z/ dz 





getting 














dz? dz dz : 2? 





The recurrence formula for the coefficients in the solution 





VW= > cr” 


p=0 


2pcp= {(p+l)(p—l— 1) —€(Bim +1) }¢pirte(2p— 1)cp-2—€*{ (p— 1)?—m?*}cy_3, 







and Y may be written in the form 





even odd 


y=- > (—e*)"/*C,i-™ n—244 7 (— 6?) (°—-DI2C p-m ta? 


€ on 





where C, has the same significance as before. 

In calculating the coefficients a, of y; or Y; the manner in which successive coefficients become of 
higher and higher order in e? provides a perfect check on the arithmetic as well as on the previous 
calculation of aim. The same is true of the calculation of the coefficients c, in the solution for large ¢, 
although here we have an additional check in the vanishing of those C,’s corresponding to negative 
powers of ¢ in y;. In every case at least one such vanishing coefficient was calculated as an arithmetical 
check. The arbitrary coefficient in the solution for large ¢ was chosen in each instance so as to make 
Yim(n) equal to (1/p)e* at infinity, whereas that in the solution for small ¢ was determined by making 
ao equal to unity. 

We shall designate by pim(n) and gim(n) the two independent solutions of the equation for vim(n) for 
small ¢ corresponding, in the case /—m even, to y; and ye, respectively, and, in the case /—m odd, to 
Y, and F¥2, respectively. The function p;,,(y) remains finite at £=0, while gim(n) becomes infinite. For 
all values of the indices the latter has the form 

















; (1+e))+1 
Gim(n) = + Pim(n) log G+ey—i Jin 












Hence, in tabulating these functions it is sufficient to give pim(n) and fim(n). In Table II the three 
functions “im(£), Pim(n), and fim(n) are given for those values of the indices for which the characteristic 
values are listed in Table I. They are expressed in such a form that the first term in each of the 
series A , or B, which has been calculated is unity. 

Finally we designate the function vm(7) representing the diverging wave (1/p)e’* at infinity by 
Yim(n) and write 


(i)'*! (+m) le! 
2-4-6-+-2mL(2/+-1) {12-3?-52---(2/—1)?} 











rim(n) = 2 imPim(n) 





214-1) {12-32-52---(27—1)2 
+(-1)g = bat) (11-32) 


(l—m) 'e+! 











The coefficients ai», and bi are listed in Table III. Evidently the function representing the converging 
wave (1/p)e—* at infinity is just the complex conjugate of rim(n). 
It will be observed that, for m=1, the equation 


@nbu=an/l(l+1) 
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holds for all five cases. This relation provides a valuable check on the arithmetical computations, 
Also, for m=2, we have a check in the relation 


Ai2b12 = 1/B’, 


although here each coefficient was calculated by two independent methods to preclude the possibility 
of arithmetical error. 


4. ZEROS OF CERTAIN CHARACTERISTIC VALUES (EIGENVALUES) 


The zeros of the coefficients b;; are very important in the applications of the theory to the antenna 
problem. The series for 6;; and for the characteristic value a1, have been carried far enough to show 
that they have a common zero e= w/2 accurate to five significant figures. To a lesser degree of accuracy 
the given series indicates that 62; and a2; have a common zero e= 7 and b3; and a3; a common zero 
¢= 3/2. We shall now prove rigorously that 6;, and a;, have a common zero and shall find its exact 


value. 
If we make m=1 in Eq. (II-14) for v and put v=2/(n?—1)!, the resulting equation for z is 


2 


( » : . 
S$ am — “sg i— =(). 
n i € az 


If €9 is a zero of a, then, the complete solution of (II-14) for this value of the parameter is 





exp (teon) + exp (—1eon), (11-33) 





v= 
(n?—1)! (n?—1)! 


and the solution representing a diverging wave of amplitude 1/p at infinity, where p=e«(n?—1)}, is 





vi(n) = (ea) exp (teon), (11-34) 


P\€0 
exact for all values of 7. This function, therefore, must be identical with the function rj(7) given by 
(1I-32) for ¢=€o, for all valués of the index /. Consequently the coefficient of the logarithmic term in 
gu(n) must vanish for ¢=¢o, that is, 
bir(€o) = 0, 


proving that a, and b;; have a common zero. 
Furthermore, if we replace 7 by (1+/)! in (II-34) and use the power series expansion for the 


exponential factor, we find, on separating real and imaginary parts, 


1 
v1(n) ae COS €9— ft Sin €o+$f*{Sin €o—€0 COS €o} +--+ 
€ 


3 
Hite — sin mee. hile €9 — €9 COS €o} + Zs (3 —€0?) sin €9 — 3€o Cos €9} +: +} (I]-35) 
€ot 2 €& 8 € 
Now the function ~::(n) contains no term in negative powers of t. Therefore, comparing with (I1-32), 
we find that when / is odd, cos ¢) =0 and consequently €9 is equal to (n+4)z, whereas, when / is even, 
sin ¢)=0 and consequently ep) =z. The value of the integer m is easily ascertained from the series for 
the characteristic values given in Table I. We find, for / odd or even, that ¢9=/x/2. 
When / is odd the part of ri(n) containing pi(m) is real, and the first term in the power series for 
this function is ¢. Hence, comparing (II-32) with (II-35), we find that 


(20+ 1) {18-3%-5#- - -(2h—1)9} (11-36) 
(1+1) leo! 





a11(€0) = 
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TABLE III. Coefficients of diverging wave functions. 








1-—-——-é+ é— 6+ +++ =1—0,020,0002+-0.000,763<—0.000,026c°+ ---, 
2-5® 28-54-72 24-34-58. 7? 





a= 


19 2609 32,593 
bu=1l- é P 


_ e+ € 
2-52 28-54-72 24-34-55.732 


6+ +++ =1—0.380,000e — 0.010,649e' + 0.000, 164«°+ ---; 





3 389 
a a v +++ = 1—0.030,6122+0.000,750e+---, 
2.72 28.33.74 


19 1751 


2 me 





— +++ = 1—0,064,6262—0.003,376e+ ---; 


€ € 
2-3-7? 28.33.78 
23 113,549 
en=1- é- é4 +++ =1—0.017,0372—0,000,257¢+---, 
2-38-52 28.36.54. 412 





37 42,233 
—s ¢- f+ ++ =1—0,027,4072—0.000,096+ - --; 
2-33-5? 28-36-54. 112 





181 
a4; = 1————— +--+» = 1—0.010,7102+--+, bywi— @+-+-=1-0.015,2642+4-+-; 
2-72-4112 2-78-11? 


67 89 
ay; =1—--———# + - -- = 1—0.007,3422+ +++, du =i— 24... =1—0,009,7522+-++; 
2-33-13? 2-38-13? 


1 121 
ax2= 1—-——e+ e+ ++ =1—0.010,204e+-0.000,233e+ +++, 
2.72 28.33.74 


53 14,381 
be =1— e+ ef -++ =1—0,180,2722+0.027,730e+---; 
2-3-72 23.33.74 


1 223 
ay2=1— e+—— + +>» = 1—0.018,5192+0.000, 3164+ ---, 
2-38 23-38-11 


19 5903 
bsp = 1 -———¢ + ————_-  - - - = 1— 0.070, 3702+. 0.001,673<'+--- 
2-38-5 23. 38-5-112 











When | is even, the part of r1:() containing p1:(y) is imaginary, and the first term in the series for this 
function is (1+2*)#. Therefore, comparing again with (II-35), we find that (11-36) holds for / even as 
well as / odd. 

When / is odd, it is evident from Table II that the first term in fi:(n) is — {2(1+#°)#}/ {ant}. Hence 


bir (€0) ™ (l— 1) leg! 
cin(€o) (20-1) {12-32-5?---(27—1)?} 


On the other hand, when / is even, the first term in fi:(y) is —2/ {ant}. Therefore, (11-37) holds for / 
even as well as for / odd. 

Similarly we can find exact values of the other coefficients in p1:(y) and fi:() for «=o. Thus we have 
a criterion of the rapidity with which the series given in the tables converge in the neighborhood of ¢o. 
This is illustrated in Table IV, where the exact values and the series values of some of the coefficients 
are given for e=¢o. The figures in parentheses indicate the number of terms used in the series. Of 
course the agreement becomes worse the larger /, partly because ¢) becomes larger and partly because 
fewer terms of the series have been calculated. In fact, the functions for ]>1 were not computed 


(1-37) 
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TABLE IV. Coefficients at zero of an. 














Series value 





Coefficient Exact value 


1 a/2 an 0.95493 --- 0.95491 











A; 1.013212--- 1.013210 (4) 
bisr/anr 0.5236: -- 0.5241 (4) 
Bz 2.026: -- 2.014 (3) 
B, 3.06+ ++ 3.26 (2) 


2 . an 0.76-+- 0.77 (3) 


A, 1.0639: -- 1.0637 (4) 

bei /an 0.219--- 0.220 (3) 

1 0 0.5 (3) 

3 3n/2 @31 0.63--- 0.71 (4) 
A, —0.200 —0.195 (4) 

bs:1/asi 0.133- ilies 0.137 (3) 

1 —0.07--- 0.05 (3) 





————=—=====_[—apO»O)")COoO*COw*"_—===—_—_———————~— ——— 





primarily for evaluation at the zeros of their characteristic values, but rather as correction terms in 
the neighborhood of the important first resonance (e= 2/2) for ¢>0. Here the convergence is entirely 
satisfactory. Therefore, the coefficients for which it is important that the series should give accurate 
values in the neighborhood of ¢9 are those for which /=1. In judging these, it must be remembered 
that the terms in f1;(7) involving Bz and B;, as compared with that in B), have numerical coefficients 
1/10 and 1/280, respectively. Hence, for small ¢, the accuracy is not less than that of the series for 
by;/a11, the first four terms of which show an error of only one part in one thousand. 

Turning now to Eq. (I1-13) for u for the case m = 1, we note that, if we put u=2/(1—&)4, we get for 


z the equation 
(1 |S +es]+ 0 
— )| —+es |+az=0, 
dé? 


which yields a solution, when a=0, of the same form as that previously obtained in yn. Hence we con- 
clude, when 7 is odd, that 


fault howe (= pelt) {12.32-52---(2/—1)*} cos ecg 


5 


(14-1) !eo'+4a11(€0) (1—e2)!" 





and, when 7 is even, that 


[un(€) ] ( pare = eS — 
u(§) je=eo= (— /2 


(1+-1) leo'+4ar1(€0) (1—#2)) 


Using the value of ai(€o) specified by (II-36) these reduce to 











2 cos €ok 
e=eg = (— 1) 2— — 11-38 
peli [11() Je=eo= (—1) = ae) (11-38) 
or d, anc 
2 si 
[uii(&) Je=eo= (— 1)? 2— = = (11-39) 
€0 =~ 


for / even. 
The method followed in this section does not succeed when m>1. 
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The straight wire antenna is approximated by a perfectly conducting prolate spheroid of 
eccentricity very close to unity. (1) The first three harmonics of the axially symmetrical free 
oscillations are discussed, with calculation of the wave-length \ and the logarithmic decrement 6. 
(2) The tangential electric field component of a plane wave with E parallel to the axis of the 
spheroid is expanded in the prolate spheroidal wave functions of Paper II. (3) The oscillations 
forced by such a plane wave are discussed for frequencies in the neighborhood of resonance with 
reference to current distribution, radiation resistance, and scattered field. 














1. FREE OSCILLATIONS 






E shall discuss briefly the axially symmetrical free oscillations of a perfectly conducting prolate 
spheroid so thin that the squares and higher powers of the value t of =(n?—1)! at the 

surface! can be neglected as compared with unity. Under these circumstances /o represents effectively 

the ratio of the minor semi-axis b of the spheroid to the major semi-axis a. We are concerned with the 

wave functions (11-15), in which 9,; stands for the function 7, representing a diverging wave as given 

by (1-32). 

The sole boundary condition is E;=0 for t= to, that is, 








d 
qa) =(Q), =e. 









a. Fundamental or First Harmonic 


For /=1, Eq. (III-1) becomes 





bu=(4, S)tel ai, (111-2) 





where 








(1+¢t2)'+1 i 9 886 | 
ELLE ee 4 64 ees 
5 7 


li:=] lo —— 
’ S i+t?)'—1 53.7 3458.7" 





if we neglect squares and higher powers of fo. 

For the limiting case to=0 of an infinitely thin antenna /,,;=0 and hence 6,;;=0. Hence, as shown 
in II, e=2/2. Consequently, if \ is the wave-length and 6 the logarithmic decrement of the free 
oscillations, 
















A\=4a, 6=0. (111-3) 






To discuss antennas of finite thickness we put 49=(10)~-" and make e= 2/2 everywhere except in 
the factor of b;, which vanishes for this value of «. This approximation yields results accurate to 
within one percent for n>3 and to within two or three percent for n=2. Then we find from (III-2) 


that 
‘| al 
e=— 1 —-————__ I. 
2 n—0.228 







‘See Paper II for definitions of the symbols appearing in this paper. References to Paper II will be indicated by II 
and to Paper I by I. 
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In terms of the wave-length and the logarithmic decrement 


aE | III 
=e =§-— i. - 
’ 2X 2r : * 
Therefore 

A\=4a, 6=1.059/(n—0.228). (111-5) 


b. Second Harmonic 


For /=2 we have the boundary condition 


4 
ba ogi tat, (III-6) 


where 
= 


(4+t)'-1 37° 387 





1+¢7)!+1 1 95 - 
In=| tog +to?)*+ 3 Heo | 


For a thin antenna this gives approximately 
1 0.1077 
” n—0.5)' 


A\=2a, 6=0.67/(n—0.5). (111-7) 





from which it follows that 


c. Third Harmonic 


For /=3 the boundary condition becomes 


2 
baa ca. 7aht ata (III-8) 


; [ (1+t0?)'+1 i ae 59 ‘i | 
=| log —____ - — ¢- e+ee-]. 
Le i tte)i—1. 3 2-3-5) 2-3-5811 


where 





As the series for 53; converges too slowly for e=32/2, we have used (II-37). For a thin antenna we 


get, approximately, 
| 0.0862 | 
e=—1 1-— 
2 n—0.4 


A= (40/3), 5=0.54/(n—0.4). (111-9) 


The series given in II have not been carried far enough to calculate even rough values of the 
logarithmic decrement for higher harmonics than the third. However, the wave-length for a thin 
antenna is given by \=4a// for the /th harmonic, however great the integer / may be. 





from which it follows that 


2. INCIDENT WAVE 


Preparatory to the investigation of the electrical oscillations of a perfectly conducting prolate 
spheroid forced by a plane electromagnetic wave of angular frequency w with the electric vector 
parallel to the long axis of the spheroid, it is necessary to expand the field components of the incident 
wave in the prolate spheroidal wave-functions specified by Eqs. (II-15) to (II-18). Evidently 
Vim() =Pim(7) since gim(n) becomes infinite at the origin. 

Designating unit vectors in the directions of increasing &, 7, @ by &1, m, ¢1, respectively, and 
putting s=(1—£#)!, t=(n?—1)!, e=(xu)*(wf/c) =2nf/A, as before, we find the electric and magnetic 


of 


Ey 


whe 
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field intensities in the incident plane wave are 


z_| ba itm kM sin, al1-10) 
L (W—-#)' (9? —#)! 

wn Base ag 108 # 
. 6 (PF) (n? - FI 
for a wave of unit amplitude. As usual we have omitted the time factor e-". 


Expanding the exponential as a power series we find for the two components of the electric field 
of the incident plane wave: 


Hi= 








— ¢1 sin a sais (III-11) 




















n ef? 14. et! . €®t6 14 + —e of: , e 4 
‘=x —__—_—_| s——s5 s5— S'+-+-++4e sin $j ts?— si... 
Ee (n?— &)! . 22 22.42 22.42.62 . 2.22 
é et é 
tH cos 26| Ast a vee i sin 3@{#st—..-}+4.. | 
| id a 4 et? . roe et! ’ 
EY =(a_ Bp) ‘75 awry ~ 92.42.62 + +++ -#e sin {Ass  diiale 
é 15 . 
te COs 29 Ast — ot oes +i sin 3p{fs?—.--}4.. | 


Since the boundary conditions involve E;‘ only, we need concern ourselves with this component 
alone. In order to obtain the coefficients with the desired accuracy, it is necessary to calculate the 
first few terms of u7,(£) and p7:() in addition to the functions given in Table II of II. 

We find for the desired expansion 


3 
72.92 


i 








"ui ; 42 2M) : ‘Mi : 
{= roman ($)7 lou) } +e sotan(§)— tpar(n)} + é isn (toa) ++ --| 


_ 1 
~ (?— 2)! 


E; 
° . n 4? 2 | ,? gE , 
+e sin 5 Man )Pu(n) + oo Maitar(€)Pai(n) + tee —Maimar(£)pa1(0) 


2 g , é 7. ne 
—— ~OM gittai(E)Par(n) +: - |-5 cos 26f *Msiuas(8) a(n) + tee 
7? 5 2? Ss 


—*aiauss(®)Pue(n) - v + vee , (III-12) 


where 
—_ 2, 93 ; 446 . 
ee Tye gage tT 





. 32 742 
My=1 —_ é 
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Mu=i+---, Masi+---, Mimi+---, Migmi+... 















LEIGH PAGE 






3. THIN RECEIVING ANTENNA 


In this section we shall investigate the current produced in a very eccentric perfectly conducting 
prolate spheroid by an incident plane electromagnetic wave of angular frequency w and unit ampli. 
tude, the electric field of which is parallel to the long axis of the spheroid. We shall be particularly 
interested in the frequency of resonance and frequencies close to resonance. Since the eccentricity 
of the spheroidal conductor under consideration is close to unity, the value to of the parameter ; 
at its surface is effectively equal to the ratio of the minor semi-axis } of the conducting spheroid to 
the major semi-axis a. We shall confine our discussion to conductors for which tp<(10)~*, with par. 
ticular attention to the limiting case attained as to goes to zero. Therefore we can neglect terms of 
order ¢o? as compared with terms of order unity. 

Remembering that the first term in p;;(n) for / odd is #, it is clear, then, that we need retain only 
this first term in the part of the expression (III-12) for E;* which is independent of ¢. As regards the 
part of this expression involving sin ¢, terms in the first power of ¢ appear, but this part of the ex. 
pression determines only the distribution of the current in azimuth and has no effect on the total 
current through any cross section of the conductor. Therefore we shall neglect it. The same is trye 
of the part of the expression involving cos 2¢, which, however, contains no terms of order less than 
# and is negligible on that account. Similar considerations apply a fortiori to the remaining parts of 
the expression. 

Hence we write as a sufficiently good approximation 





















, 2 , 
B= | Mit) +—é M30ttai(&) + 
(n?— &)4 “ll 


The corresponding component of the electric field intensity in the scattered wave must be 





we gat Msomsi(8) + . + (ITT-13) 






Bro. —1 4 +E oot 8) btra(n) 
F -~"..| om(S)> a rii(m)} HE Ngousi(é , rsi(n)} 


(n?— &)! id 
+e Neosi(8)- qr }+-:: | (IIT-14) 





in order to satisfy the boundary condition E;'+E;"=0 for t=t%9. Consequently, if terms in # and 
higher powers are neglected, 






qf ludu—tbn 











, 2 , 
Nio= ehuMio : , ’ 
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where /,; and /3; are defined in Section 1 of this paper, and 


1+¢9?)'+1 137 1 ‘ta 
inm| alld = ete] ‘ 
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The magnetic field strength in the scattered wave corresponding to (II1-14) consists of the single 


component 


} 
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and the current 1s 
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with inclusion of the time factor. 
We consider first the limiting case of an antenna of length 2a=2f and zero thickness. Then 


hy, =/31=/51=0, and the current is zero except for those wave-lengths which make one of the coef- 
ficients bi, 631, bs: vanish. It was shown in II that 6,; vanishes for e= 2/2, bs; vanishes for e= 32/2, 
and 65; vanishes for e=5x/2. Hence, as «=22f/ and a=f when to=0, it follows that the first three 
resonances occur at wave-lengths equal to twice, two-thirds, and two-fifths the length of the antenna. 
The resonance is infinitely sharp but not infinitely high. 


a. First Resonance 


We shall discuss the first resonance in some detail, as it is by far the most important. First consider 
an infinitely thin antenna. Making use of Eqs. (II-34) and (II-38), we find for the current at first 


resonance 


(111-17) 
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in exact agreement with the result of the approximate method employed in I. 
Since the element of length dd; corresponding to the increment dé is 


e—e 3 
n= s(“—) ae 


the applied electromotive force along d); is 


, adg , 2 ’ ’ 
E;‘dd: = “| M jot (€) +—e M gots: (£) + €'M sous (€) + snail COS WT. 
— - be a 


72.92 
Remembering that the functions u,,(&) are orthogonal, we find that the mean time rate of absorption 


of energy is 
} 1 
r(*) a Mio(e) Ff [ 111(€) Pdé. 
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The radiation resistance R is defined as the ratio of this quantity to one-half the square of the current 
amplitude at the center of the antenna. Hence, 












R (111-18) 
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If R, is the radiation resistance in ohms 





i 
Ry = 4xc4(10)-*R=73.08(“) ohms, (IIT-19) 
K 





in exact agreement with the result of the approximate method employed in I. 
Finally we obtain, with the aid of (II-34) and (II-38), the following simple expressions for the 
three non-vanishing field components in the wave scattered by an infinitely thin antenna at resonance: 
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where Mjo(eo)=0.8206 and, for the limiting case of an infinitely thin antenna, 7.=1. These 
expressions are valid at resonance for all »>1. 

Next we pass to the case of an antenna of small but finite thickness. As in Section 1 we put 
to=(10)-*. If n»>3 we do not incur an error as great as one percent if we retain the expressions for 
the wave-length of resonance and for the radiation resistance at resonance which were found for 
the infinitely thin antenna, and even for »=2 the error is only about two percent. But now the 
resonance peak in the current is no longer infinitely sharp. To determine its form we shall use only 
the first term in the general expression (III-16) for the current, ignoring a small but almost neg- 


ligible asymmetry due to succeeding terms. We find 
rT 
1.045 cos { —é 
x\? ] 2 
I=|2xrc{-—}a 
uJ J 


6A 
1+35.2(n— 0.228)*(— 
for the current, where \, = 4a is the wave-length of resonance and 6d the deviation from \,. The lead 
@ of the current ahead of the electromotive-force is given by 
6A 
i. 
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tan ¢=5.93(nm—0.228) 






These expressions are valid, of course, only for values of the ratio 6/A, small compared with unity. 
They represent quite typical resonance curves. 











b. Second Resonance 





In discussing the second resonance we shall consider only the case of an infinitely thin antenna, 
and, as the series involved converges more slowly in this region, we can do no more than calculate 
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a rough value of the current amplitude at resonance. We find for the current at resonance 


«\? 3x 
[= -| 2xc(*) « o.04 cos (Fs) COS wr. 
m 2 


The current is in phase with thé impressed electromotive force in the two extreme thirds of the 
antenna, but out of phase in the middle third. As the current amplitude at the center of the antenna 
is only some 4 percent of that at first resonance, the second and higher order resonances are evidently 
of little importance as compared with the first resonance. 


(III-25) 
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Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition 
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The partition function of a two-dimensional “ferro- 
magnetic” with scalar “spins” (Ising model) is computed 
rigorously for the case of vanishing field. The eigenwert 
problem involved in the corresponding computation for a 
long strip crystal of finite width (m atoms), joined straight 
to itself around a cylinder, is solved by direct product 
decomposition; in the special case n= an integral 
replaces a sum. The choice of different interaction energies 
(+J,+J’) in the (01) and (10) directions does not 
complicate the problem. The two-way infinite crystal has 
an order-disorder transition at a temperature T= 7, given 
by the condition 


sinh(2J/kT.) sinh(2J’/kT,) = 1. 


INTRODUCTION 


HE statistical theory of phase changes in 
solids and liquids involves formidable 
mathematical problems. 

In dealing with transitions of the first order, 
computation of the partition functions of both 
phases by successive approximation may be 
adequate. In such cases it is to be expected that 
both functions will be analytic functions of the 
temperature, capable of extension beyond the 
transition point, so that good methods of ap- 
proximating the functions may be expected to 
yield good results for their derivatives as well, 
and the heat of transition can be obtained from 
the difference of the latter. In this case, allowing 
the continuation of at least one phase into its 
metastable range, the heat of transition, the 
most appropriate measure of the discontinuity, 


The energy is a continuous function of 7; but the specific 
heat becomes infinite as —log |T—7.|. For strips of 
finite width, the maximum of the specific heat, increases 
linearly with log. The order-converting dual transfor- 
mation invented by Kramers and Wannier effects a simple 
automorphism of the basis of the quaternion algebra which 
is natural to the problem in hand. In addition to the 
thermodynamic properties of the massive crystal, the free 
energy of a (0 1) boundary between areas of opposite order 
is computed; on this basis the mean ordered length of a 
strip crystal is 


(exp (2J/kT) tanh(2J'/kT))*. 


may be considered to exist over a range of 
temperatures. 

It is quite otherwise with the more subtle 
transitions which take place without the release 
of latent heat. These transitions are usually 
marked by the vanishing of a physical variable, 
often an asymmetry, which ceases to exist 
beyond the transition point. By definition, the 
strongest possible discontinuity involves the 
specific heat. Experimentally, several types are 
known. In the a—§8 quartz transition,! the 
specific heat becomes infinite as (T.—7)~; this 
may be the rule for a great many structural 
transformations in crystals. On the other hand, 
supraconductors exhibit a clear-cut finite discon- 
tinuity of the specific heat, and the normal state 
can be continued at will below the transition 


1H. Moser, Physik. Zeits. 37, 737 (1936). 
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point by the application of a magnetic field.? 
Ferromagnetic Curie points’ and the transition 
of liquid helium to its suprafluid state‘ are 
apparently marked by essential singularities of 
the specific heats; in either case it is difficult to 
decide whether the specific heat itself or its 
temperature coefficient becomes infinite. Theo- 
retically, for an ideal Bose gas the temperature 
coefficient of the specific heat should have a 
finite discontinuity at the temperature where 
“Einstein condensation” begins. 

In every one of these cases the transition 
point is marked by a discontinuity which does 
not exist in the same form at any other tempera- 
ture; for even when a phase can be continued 
beyond its normal range of stability, the transi- 
tion in the extended range involves a finite heat 
effect. 

Whenever the thermodynamic functions have 
am essential singularity any computation by 
successive approximation is difficult because the 
convergence of approximation by analytic func- 
tions in such cases is notoriously slow. 

When the existing dearth of suitable mathe- 
matical methods is considered, it becomes a 
matter of interest to investigate models, however 
far removed from natural crystals, which yield to 
exact computation and exhibit transition points. 

It is known that no model which is infinite in 
one dimension only can have any transition.® ® 

The two-dimensional “Ising model,” originally 
intended as a model of a ferromagnetic,’ is known 
to be more properly representative of condensa- 
tion phenomena in the two-dimensional systems 
formed by the adsorption of gases on the sur- 
faces of crystals.* It is known that this model 
should have a transition.* Recently, the transi- 
tion point has been located by Kramers and 
Wannier.'® In the following, the partition func- 


2W. H. Keesom, Zeits. f. tech. Physik 15, 515 (1934); 
H. _ and R. O. Wilhelm, Rev. Mod. Phys. 7, 237 
(1935). 

3H. Klinkhardt, Ann. d. Physik [4], 84, 167 (1927). 

*W. H. and A. P. Keesom, Physica 2, 557 (1935); 
J. Satterly, Rev. Mod. Phys. 8, 347 (1938); K. Darrow, 
sbid. 12, 257 (1940). 

5 E. Montroll, J. Chem. Phys. 9, 706 (1941). 

*K. Herzfeld and M. Goeppert-Mayer, J. Chem. Phys. 
2, 38 (1934). 

7E. Ising, Zeits. f. Physik 31, 253 (1925). 

®R., Peierls, Proc. Camb. Phil. Soc. 32, 471 (1936). 

*R. Peierls, Proc. Camb. Phil. Soc. 32, 477 (1936). 

10H. A. Kramers and G. H. Wannier, Phys. Rev. 60, 
252, 263 (1941). 


tion and derived thermodynamic functions wil] 
be computed for a rectangular lattice with 
different interaction constants in two perpen- 
dicular directions; this generalization is of some 
interest and does not add any difficulties. 

The similar computation of the partition 
functions for hexagonal and honeycomb lattices 
involves but relatively simple additional con. 
siderations. A general form of the dual transfor- 
mation invented by Kramers and Wannier 
together with a rather obyious ‘“‘star-triangle” 
transformation are used in this connection, and 
the transition points can be computed from the 
transformations alone. While these were a great 
help to the discovery of the present results, the 
logical development of the latter is best relieved 
of such extraneous topics, and the subject of 
transformations will be reserved for later 
communication. 


OUTLINE OF METHOD 


As shown by Kramers and Wannier," the 
computation of the partition function can be 
reduced to an eigenwert problem. This method 
will be employed in the following, but it will be 
convenient to emphasize the abstract properties 
of relatively simple operators rather than their 
explicit representation by unwieldy matrices. 

For an introduction to the language we shall 
start with the well-known problem of the linear 
lattice and proceed from that to a rectangular 
lattice on a finite base. For symmetry we shall 
wrap the latter on a cylinder, only straight 
rather than in the screw arrangement preferred 
by Kramers and Wannier. 

The special properties of the operators in- 
volved in this problem allow their expansion as 
linear combinations of the generating basis 
elements of an algebra which can be decomposed 
into direct products of quaternion algebras. The 
representation of the operators in question can 
be reduced accordingly to a sum of direct 
products of two-dimensional representations, and 
the roots of the secular equation for the problem 
in hand are obtained as products of the roots of 
certain quadratic equations. To find all the roots 
requires complete reduction, which is best per- 
formed by the explicit construction of a trans- 
forming matrix, with valuable by-producis of 
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Fic. 1. Linear crystal model. Spin of entering atom O 
depends statistically on spin of atom in previous end 


position @. 


identities useful for the computation of averages 
pertaining to the crystal. This important but 
formidable undertaking will be reserved for a 
later communication. It so happens that the 
representations of maximal dimension, which 
contain the two largest roots, are identified with 
ease from simple general properties of the 
operators and their representative matrices. 
The largest roots whose eigenvectors satisfy 
certain special conditions can be found by a 
moderate elaboration of the procedure; these 
results will suffice for a qualitative investigation 
of the spectrum. To determine the thermo- 
dynamic properties of the model it suffices to 
compute the largest root of the secular equation 
as a function of the temperature. 

The passage to the limiting case of an infinite 
base involves merely the substitution of integrals 
for sums. The integrals are simplified by elliptic 
substitutions, whereby the symmetrical param- 
eter of Kramers and Wannier’ appears in the 
modulus of the elliptic functions. The modulus 
equals unity at the ‘Curie point’; the conse- 
quent logarithmic infinity of the specific heat 
confirms a conjecture made by Kramers and 
Wannier. Their conjecture regarding the varia- 
tion of the maximum specific heat with the size 
of the crystal base is also confirmed. 


THE LINEAR ISING MODEL 


As an introduction to our subsequent notation 
we shall consider a one-dimensional model which 
has been dealt with by many previous authors.” !° 
(See Fig. 1.) 

Consider a chain of atoms where each atom (k) 
possesses an internal coordinate un which may 
take the values +1. The interaction energy 
between each pair of successive atoms depends 
on whether their internal coordinates are alike 
or different: 


u(1,1)=u(—1, —1)=—x(1, —1) 


= —u(—1,1)=—J, 
or 


uu, wD) = —Ju ye, (1) 
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The partition function of this one-dimensional 
crystal equals 


eF kT = 0(T) (2) 
- ) 2 exp(— 2 u(u“, w*-Y)/RT). 


aM.) Oe MM was 


Following Kramers and Wannier,'® we express 
this sum in terms of a matrix whose elements are 


(u| V | go’) =e ew DET, (3) 


Then the summation over (u®, «++, u«-) in (2) 
is equivalent to matrix multiplication, and we 
obtain 


QAT)= LD (u| Vy) =rAnVO(1), (4) 


#1, BN 


where A», is the greatest characteristic number of 
the matrix V. Since all the elements of this 
matrix are positive, A» will be a simple root of 
the secular equation 


'V—A| =0, 
and the identity 
Am¥m(u) = (V, Ym(m)) (S) 


will be satisfied by a function (eigenvector) 
Ym(u) which does not change sign.“ In the 
present case we have 


(6) 
Am = 2 cosh(J/kT). 


It will be convenient to use the abbreviation 
J/kT=H. (7) 


The matrix V, written explicitly 


e, «4 
J = i eH ). (8) 
represents an operator on functions ¥(u), which 
has the effect 


(V, W(u)) =e*p(u) +e-#y(— 4). (9) 


This operator is therefore a function of the 
“complementary’’ operator C which replaces yu 
by —u 

(C, Wu) =¥(—x); 


C*=1, (10) 


1S. B. Frobenius, Preuss. Akad. Wiss. pp. 514-518 
(1909); R. Oldenberger, Duke Math. J. 6, 357 (1940). 
See also E. Montroll, reference 5. 
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_ Fic. 2. Cylindrical crystal model. Configuration of added 
tier-—O—O—O— depends on that of tier — @— @— @— 
in previous marginal position. 


and whose eigenfunctions are the (two) different 
powers of yu 
(C, 1)=1, (C, u)= —FH, (11) 


with the characteristic numbers +1 and —1, 
respectively. The linearized expression of V in 
terms of C is evident from (9), which is now 
written more compactly 


V=e%+e-7C. (12) 
The eigenfunctions of V are of course the same 


as those of C, given by (11); the corresponding 
characteristic numbers are 


Am =A, =e%+e-4=2 cosh H 


- : (13) 
A_ =e4 —e-#¥ =2 sinh H. 


We shall deal frequently with products of many 
operators of the type (12). It is therefore a 
matter of great interest to linearize the logarithm 
of such an operator. Following Kramers and 
Wannier,”® we introduce the function of H 


H*=}3 log coth H=tanh (e”), ~— (14) 
and obtain 
V=e" (1+e-°"C) =e” (1+(tanh H*)C), 
which may be written 
V =(e#/cosh H*) exp (H*C) 
=(2 sinh 2H)‘ exp (H*C). (15) 
The first, less symmetrical form serves to keep 


a record of the sign when H is not a positive 
real number. For example, 


(—H)*=H*44mi 


16 
cosh((— H)*) = +7 sinh H*. (16) 


For reference, we tabulate here the relations . 
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given by Kramers and Wannier. 


sinh 2H sinh 2H*=cosh 2H tanh 2H* 
=tanh 2H cosh 2H*=1. 






(17) 


Alternatively, we can express the hyperbolic 
functions in terms of the gudermannian angle 





ed u=2 tan—(e") —34=2 tan— (tanh 4x), (18) 






whereby 





sinh u cot(gd u) =coth u sin(gd u) 
=cosh u cos(gd u)=1. (19) 







Then the relation between H and H* may be 
written in the form 


gd 2H+gd 2H*=}r. (20) 









THE RECTANGULAR ISING MODEL 






We consider parallel chains of the same type 
as before. We shall build up these chains simy. 
taneously, tier by tier, adding one atom to each 
chain in one step (see Fig. 2). With the last atom 
in the jth chain is associated the variable ,, 
capable of the values +1. A complete set of 
values assigned to the m variables describes a 
configuration (u) = (u1, --+, un) of the last tier of 
atoms. The operator which describes the addition 
of a new tier with the interaction energy 












ui((u), (u’)) = -5 Jum =—-kTHD um; (21) 





is now 





Vi=I (e*+e-*C)) 
_ = (2 sinh 2H)? exp(H*B), (20) 
whereby 








B=¥ C, (23) 


jel 






and the individual operators C;, ---, C, have 
the effect 


(Ci, V(m1, "ee Hn)) 
=(u1, °° 


Next, let us assume a similar interaction betwee 
adjacent atoms in a tier, only with an inde 
pendent value J’ for the pairwise energy d 
interaction. For symmetry, let the mth atom k 
neighbor to the first; then with the conventia 






*» Mj-1) ~My, Mj+l) ** "5 Mn). (24) 
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that the atoms in a tier are numbered modulo n, 
Mjtn= His 
the total tierwise interaction energy will be 


tie(May ***» Hn) 


=—J' YD win = —kTH’ 2D ug (25) 
j=l j=l 

The effect of this interaction is simply to multi- 
ply the general term of the partition function, 
represented by one of the 2" vector components 
Y(u1, ***, Mn), by the appropriate factor 
exp(—we(u1, ***, Hn)/RT); the corresponding 
operator V: has a diagonal matrix in this 
representation. It can be constructed from the 
simple operators s;, ---, S, Which multiply y by 
its first, ---, mth argument, as follows 


(Sj, (ma, ++, on) =M(ur, +++, on) (26) 
A =r SiSit1y (27) 
V2=exp(H’A). (28) 


We now imagine that we build the crystal in 
alternate steps. We first add a new tier of atoms 
(see Fig. 3), next introduce interaction among 
atoms in the same tier, after that we add another 
tier, etc. The alternate modifications of the 
partition function are represented by the product 


+++ Va Vi Va Vi Ve Vi 


with alternating factors; the addition of one tier 
of atoms with interaction both ways is repre- 
sented by (V2Vi)=V. The eigenwert problem 
which yields the partition function of the crystal 
is accordingly 


W=(V,¥)=(V2Vi, y) 
= (2 sinh 217)"(exp(/7’A) exp(H*B), y). (29) 


Viz., in the more explicit matrix notation used by 
previous authors” 


AV (ur Hn) =exp(17 Lv ams) 
1 


xX 2 exp( > ums ) ¥(ur'* + + pn’). 
wi? spy’ = +1 1 
® See H. A. Kramers and G. H. Wannier (reference 10). 
E. Montroll (reference 5); E. N. Lassettre and J. P. Howe, 
J. Chem. Phys. 9, 747 (1941). 
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This 2*-dimensional problem we shall reduce to 
certain sets of 2-dimensional problems. So far, 
we have expressed A and B rather simply by the 
generating elements s;, Ci, «++, Sn, C, of a certain 
abstract algebra, which is easily recognized as a 
direct product of m quaternion algebras. How- 
ever, the expansion (27) of A is not linear. 
Following a preliminary study of the basis group, 
we shall therefore construct a new basis for a 
part of the present algebra, which allows linear 
expansions of both A and B in terms of ‘its 
generating elements. 


QUATERNION ALGEBRA 


We have introduced two sets of linear oper- 


ators $;, -*+, 5S, and C,, ---, C, by the definitions 
(Sj, V(u1, 7? eo Se Mn)) 

=pw(u, 7. oe oe Mn)» (30) 
(C;, ¥(u1, Le eee Mn)) 

=y(u1, tty yy ee, Mn). 


We shall see that the group generated by these 
operators forms a complete basis for the algebra 
of linear operators in the 2*-dimensional vector 


space of the functions ¥(u:, +--, ua). Their 
abstract properties are 
sPf=CP=1, siCp= —Cys;; 
Sk =SeSj, CiCe= CC; (31) 
$;Ch=Cysj; (jk). 


The algebra based on one of the subgroups 
(1, s;, C;, s;C;) is a simple quaternion algebra. 
The group of its basis elements includes (—1, 
—$;, —C;, C;s;) as well; but this completion of the 
group adds no new independent element for the 
purpose of constructing a basis. If we disregard 
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Fic. 3. Two-step extension of a two-dimensional crystal. 
(a) A new tier of atoms O is added (V;); their configuration 
depends on that of atoms @ in previous marginal position. 
(b) Interaction energy between marginal atoms O is 
introduced (V2), which modifies the distribution of con- 
figurations in this tier of atoms, 
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the sign, each element is its own inverse. Except 
for the invariant element (—1) and the unit 
proper, each element together with its negative 
constitutes a class. This algebra has just one 
irreducible representation 


1 0 1 0O 

Dit)=() i) Dd=(5 _,): 
0 1 ,. 4a 

picy=(; 5): Duoco=(", 4): 


Our algebra of operators in 2*-dimensional 
space can be described accordingly as the direct 
product of m simple quaternion algebras, which 
is still referred to as a (general) quaternion 
algebra. Each one of its 4" independent basis 
elements can be expressed uniquely as a product 
of m factors, one from each set (1, s;, Cj, s;C;). 
Still, each element furnishes its own inverse and, 
except for (—1) and the unit, each element with 
its negative constitutes a class, which immedi- 
ately characterizes all possible representations of 
the algebra: The characters of all basis elements 
with the sole exception of the unit must vanish. 
This of course is a consequence of the direct 
product construction. Another is the recognition 
that the algebra has only one irreducible repre- 
sentation 


D.XD:X---XD.=D," 


which is of 2" dimensions; our definitions (30) 
describe an explicit version of it. 

With the customary restriction that the ele- 
ments of the basis group are to be represented 
by unitary matrices, the representations of the 
4" independent basis elements must be orthogonal 
to each other (in the hermitian sense) on the 2” 
(or more) matrix elements as a vector basis. 
For the vector product in question is nothing but 


F D(Qr)eD(Qs)re 
rel =Tr. (D(Q:0:*)) =Tr. (D(Qs)), 


and if Q, and Q are independent elements of the 
basis group, then their quotient Q; is an element 
different from the unit and (—1), so that its 
character vanishes. 

Since orthogonality guarantees linear inde- 
pendence, the irreducible 2*-dimensional repre- 
sentation of our 4" independent basis elements 
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will suffice for the construction of any 2" x9n 
matrix by linear combination; our quaternion 
algebra is equivalent to the complete matrix algebrg 
of 2" dimensions. The choice between the two 
forms is a matter of convenience ; the quaternion 
description seems more natural to the problem 
in hand. One of the main results obtained by 
Kramers and Wannier will be derived below; the 
comparison with their derivation is instructive, 
In the following, the entire basis group will not 
be needed. Since our operator A involves only 
products of pairs of the generating elements 
$1, ***, S, and B does not involve any of these 
but only C, ---, C,, we shall have to deal only 
with that subgroup of the basis which consists 
of elements containing an even number of factors 
s;. This ‘‘even’’ subgroup, which comprises half 
the previous basis, may be characterized equally 
well by the observation that it consists of those 
basis elements which commute with the operator 


C=C,C2---C, (32) 

















which reverses all spins at once 


(C, W(ur, +++, wn) =¥(—m, °° 


The representation (30) is now reducible because 
the 2"—' even functions satisfying 


(C, y= 


are transformed among themselves as are the 
odd functions satisfying 


The corresponding algebras are both of the same 
type as before, with (m—1) pairs of generating 
elements 





| — pn). (32a) 












S1Sny S25ny ***y Sn—15n; Ci, C2, = Coxs 






and either algebra is equivalent to the complete 
matrix algebra of 2"~—' dimensions. The previously 
independent operator C, is now expressible by 
the others because either C=1 or C=-l, 
according to whether we are dealing with even or 
odd functions. This choice of sign determines the 
only difference between the two versions of the 
algebra in question. Moreover, the operator 
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is the same in both versions, but obviously 


B=2 C; 
j=l 
is not; here the interpretation of C, makes a 
difference. 

The arbitrary elimination of one of the sets of 
operators which enter symmetrically into the 
problem causes much inconvenience. The sym- 
metry can be retained if we start with the two 
series of operators 


$1S2, S283, ***) SnS1; Ci, C2, eis Ga (33) 


required to satisfy the identities 


S$ SoSoS53°* *$nS1 = S17So"- ° -S2=1 


(ts 


t= 1; (odd functions). 


(even functions) (34) 


C,C2:-:-C,=C= 


The commutation rules are simply described 
when the two sets are arranged in the interlacing 
cyclic sequence 


SnS15 Ci, 51Se, C2, iit Sa—1Sa Cos (S51, . -). (35) 


Here each element anticommutes with its neigh- 
bors and commutes with all other elements in 
the sequence. 


Dual Transformation 


The fundamental algebraic reason for the dual 
relation found by Kramers and Wannier is now 
readily recognized: In the representation which 
belongs to even functions, the two sets (33) are 
alike because the starting point of the sequence 
(35) is immaterial. Neither does it matter if the 
sequence is reversed, so that the general rules of 
computation remain unchanged if each member 
of the sequence (35) is replaced by the corre- 
sponding member of the sequence 


Cry SeSet1, Ceti, SktiSkt2, °° * 
Accordingly, the correspondence 
(36) 


S$ jp: Cet j, CpSee (5-1) Sk 45 


describes an automorphism of the whole algebra. 


Under this automorphism we have in particular 


A-B 
B-A 


or, compactly 
f(A, B)—f(B, A), 


all in the algebra of even functions. Similarity 
can be inferred from the additional observation 
that the characters of all representations are 
determined by the commutation rules, or equally 
well from the theorem that every automorphism 
of*a complete matrix algebra is an internal 
automorphism. When we return to our original 
mixed algebra of even and odd functions, our 
result takes the form 


(1+C)f(A, B)~(1+C)f(B, A). 


Application to the operator V given by (29) 
shows that the two operators 


(sinh 2H)-"?(1+C) V(A’, 1) 

=2-"(1+(C) exp(H’A) exp(H*B), 
(sinh 2H7’)"?(1+(C) V’(H*, H’*) 

=2-"(1+(C) exp(H’B) exp(H*A) 


(37) 


are similar, so that the part of the spectrum of 
V(A*, H’*) which belongs to even functions can 
be obtained from the corresponding part of the 
spectrum of V(H’, H) by the relation 


\,(H*, H’*) 


= (sinh 2H sinh 2H’)-"\,(H’, H). (38) 


This result applies in particular to the largest 
characteristic numbers of the operators in ques- 
tion because both can be represented by matrices 
with all positive elements so that, by Frobenius’ 
theorem," the corresponding eigenvectors cannot 
be odd. 

In the “screw” construction of a crystal 
preferred by Kramers and Wannier,’® the addi- 
tion of one atom is accomplished by the operator 
(M in their notation) 


exp(Hs,_15,)(sinh 2H)! exp(H*C,)T 

= (sinh 2H)! exp(Hs,_15,)T exp(H*C)), 
where T is a permutation operator which turns 
the screw through the angle 2x/n. They indicate 


in explicit matrix form a transformation which 


effects one of the automorphisms (36), namely, 
Sj 3417 Cn_j, C55 n— 5418 0-j 


whereby naturally 
TT, 
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The equivalent of (38) follows. In consequence, 
if the crystal is disordered at high temperatures 
(small HY), a transition must occur at the 
temperature where the equality 


H= H* =} log cot 1/8 (39a) 


is satisfied. 

Unfortunately, our simplified derivation fails 
to exhibit one of the most important properties 
of the “dual” transformation: When applied to 
large two-dimensional («*) crystal models, it 
always converts order into disorder and vice 
versa; this result is readily obtained by detailed 
inspection of the matrix representation. The 
consequent determination of the transition point 
for the rectangular crystal by the generalization 
of (39a) 


H'=H*; (kT=kT.) (39b) 


will be verified in the following by other means. 

On the other hand, the given results are still 
but special applications of a _ transformation 
which applies to all orientable graphs, not just 
rectangular nets. Except for its ‘‘order-convert- 
ing’’ property, this transformation is best ob- 
tained and discussed in simple topological terms. 
The properties of the ‘‘dual’’ and certain other 
transformations will be given due attention in a 
separate treatise; it would be too much of a 
digression to carry the subject further at this 


point. 
Fundamental Vector Systems 


We shall give brief consideration to the vector 
basis of our quaternion algebra, as described by 
the operator formulation (30). The systems of 
eigenvector which belong to the two sets of 
operators (5;, «++, Sn) and (Ci, ---, C,) are easily 
identified. The representations (30) of the former 
are already diagonal, so that their eigenvectors 


are 6 functions 
*,» Mn —pn’)) 


(Sj, 5(ui— mr’, ia Mj— Bj’, ia 


=u j'5(ui— mi, roe Mi— Hs, ain Hn — Un’). (40) 


The eigenvectors of (C:, ---, C,) are the different 
products of wi, ***, wa, with the normalizing 
factor 2-*/? 


(Cj, (ur + +n) =(—) (ur +a). (41) 
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The explicit representation of a transformation 
which leads from one of these systems to the 
other is obtained from the expansion of the 
identity 








*» Mn — pn’) 
=2-"(1-yi'ai)***(1-+Mn’un) (42) 


This product description of a 6 function is oftep 
more convenient. The two systems (40) and (41) 
yield diagonal representations of the operators 
A and B, respectively. The construction of the 
latter is a little simpler. With a slight variation 
of the notation we find 


5(ui— py’, “= 










(B, (uismie- - mix) = (nm — 2k) (usrmie- + * pix), (43a) 






so that any homogeneous function of (1, - - -, g,) 
is an eigenvector of B. With no repeated factors 
allowed, the number of independent homogene. 
ous functions of degree & is the binomial coeff. 







, n , 
cient (7). Accordingly, the secular polynomial 


1S 





|B—d| =T] (n—2k—2)). 


k=0 


(43b) 






Whereas B measures the degree of a function 
V(u1, ***, Mn), the operator A simply counts 
alternations of sign in the configuration (;, --:, 
un). The analog of (43a) is now 








Jet 


(4, I (1+u,) I (1—p,;)--- IT, (1+u)) 





= (n—4k) II (1+4;) IT (1—y,)--- 





xT (+m). (44a) 





The alternations of sign are located arbitrarily 
behind j1, j2, «+ *j2; but their total number must 
be even = 2k. On the other hand, the result (44a) 
is equally valid when yu; is replaced throughout 
by —u;; so that the spectrum of A repeats in 
duplicate that part of the spectrum of B which 
belongs to even functions, and its secular 
polynomial is . 










[n/2] . 
|A-A|= il (n—4k —y)2(2#), 


k=0 





(4b) 





By an obvious transformation, one-half of the 
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paired characteristic numbers of A may be 
assigned to even functions and the other, 
identical half to odd functions. We note again 


that the three operators 
(1-C)A~(1+C)A~(1+C)B 


are similar, whereas (1—C)B has a different 


structure. 
Some Important Elements of the Basis 


Asa preliminary to the simultaneous reduction 
of the representation of the operators A and B, 
we shall consider certain elements of the basis 
group which will recur frequently in the process 
and introduce a special notation for them. We 


have already encountered 
P..= ame Ce Poa+1=SeSe+1- 
In general we define P.» as the product 


Pp =SaCa+1Ca+2° 7s Cy—150; 


(a,b=1,2,-+-, 2m). (45) 


These operators are elements of the ‘‘even”’ basis 
and commute with C. They satisfy the recurrence 
formula 


Po041= ParPwPo,041= ParseCose41. (45a) 


The period of this recurrence formula is 2n; 
increasing either index by has the effect of a 
multiplication by (—C) 


Papin =Paino= —CPa= —P»C. (46) 


The product (45) obviously may be written in 
reverse order when (b—a)=n, and by (46), this 
result is easily extended to all (a, 6). Accordingly, 


Pip =SeCa41" + Cy_i8e8eCo-1+ ++ Ca4iSa=1. (47) 


More generally, we thus obtain the “rule of the 
wild index”’ 


PacP c= PaaP oa, PacPaa= PoeP oa: (48) 


The analogous result for indices congruent 
(mod=n) is obtained immediately by combina- 


tion with (46) 
PicPs, ct+n = PoaPr, d+-n= PoaaP rin, dy 
(49) 
PocPorn.d= PecPoin.a= PrePo, arn: 


While the period of the recurrence formula (45a) 
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is 2m, it is evident from (46) that the period of 
the commutation rules will be just ; for the 
factor C commutes with all the operators P,». 
For this reason we shall need a modified Kro- 
necker symbol with the period n, as follows 


D(0)=D(n)=1, D(m) sin(mr/n)=0. (50) 


Indices which differ by 0 or » (mod =2n) will be 
called congruent; (in the former case we call 
them equal). Only the congruence of correspond- 
ing indices matters ; no important identity results 
from the accident that a front index coincides 
with a rear index. The distinction between front 
and rear indices represents the distinction be- 
tween the opposite directions of rotation (1, 2, 3, 
+++, m) and (1, m, m—1, ---, 2) around the base 
of the “‘crystal.”’ 

The effect of simultaneous congruence is 
evident from (47) together with (46) ; one obtains 
either the unit (1) or the invariant basis element 
(—C). More generally, a product of the form 


Poa Pw: . Pree 


with any number of factors equals one of the 
invariant elements +1, +C if both front and 
rear indices congruent to each of the possible 
values 1, 2, ---, ” occur an even number (0, 2, 4, 
-++) of times. This rule can be obtained from 
(48) by simple manipulation; compare the 
derivation of exchange rules below. In addition, 
thanks to the relation 


PyyP 22° + + Pan=(—)"CiC2:--Cra=(—)"C, (51) 


an invariant element also results if both front 
and rear indices congruent to each of the numbers 
1, 2, --+, m occur an odd number of times. A 
count of the group elements shows that there 
can be no further identities of this sort; accord- 
ingly, if a product satisfies neither of the given 
conditions, it must be a basis element different 
from +1 and +C. 

The commutation rules for the operators P,» 
can be broken up into rules which govern the 
exchange of indices between adjacent factors in 
a product. These in turn can be obtained by the 
rule of the wild index from the commutation 
rules which involve P..= —C,. We first consider 
the product 


PaPun= —PuG= —SaCas Ae? * Cr—150Cp. 
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When a=), the commutative character of the 
multiplication is trivial. Otherwise, there is only 
one factor s, in P,»; so that this operator anti- 
commutes with C,= — Ps». Accordingly, 


PaPwt(—)??PePar=0. 
In view of (48) this may be written equally well 


PacPret+ (—)?°PrcPac= 0. 


However, the result does not depend on the, 


presence of a wild index, because 


PacP a= PacPaaP otP oa = PrePraP aaP va 
= —(—)?@-) P,P, aPra. 


The same reasoning applies to the exchange of 
rear indices, and we obtain the general rule 


PocPra= —(—)?° PrcPoa= — (—)? PaaPre 
=(—)P(e)+D(e-4) P,P... (52) 


Accordingly, the exchange of indices between 
adjacent factors in a product takes place with 
change of sign unless the trivial exchange of 
congruent (wild) indices is involved. 

The product P..Psa is seen to be commutative 
when neither: (a—0) nor (c—d) divides n, and in 
the trivial case when both do. When one and 
only one congruence exists, the factors anti- 
commute. 

For future reference we give here the direct 
operational definition of P.». In the case 


cot(a—4)x/n> cot(b—4)x/n 


we have 
(Pav, (ui, ***, Mn) = amo (Ma, °°, May 
Matty ***y —Mb—-1) Mb, ***, Mn); (53a) 
(sin((b—a—4)x/n) >0) 
(Pav, W(u1, °° +, Mn) = — bao (— mi, ***, 
fay Matt, ***, Mb—-1, — Mb, ***, —Hn); (S3b) 


(sin((—a—4)x/n) <0). 


When cot(a—})x/n= cot(b—})x/n, substanti- 
ally the same formulas apply, only the changes 


of sign are arranged (—y1, --*, —s-1, 

Mb, ***, May Matt, ***, —Mn) in the case (a) 
and 
(u1, ** *, Mb—1, —~ Mb, ***, “May Mat+l1, ***» btn) 


in the case (b). 
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The given formulas describe the effect of P,, 
in terms of the vector basis (40). The effect of 
this same operator on a vector which belongs 
to the fundamental system (41) can be described 
conveniently with the aid of functions which 
introduce a change of sign once in a period 
(1, 2, ---, m). Individually, these functions haye 
the period 2. We shall use the standard notation 











—-1; x<0 
sgnx=) 0; x=0 
| 1; x>0. 
Then the function 
sgn(sin(jx/n)) 






changes sign at 7=0 and j=” and vanishes for 
these values of 7, while the function 






sgn(sin(j—})x/n) 





has everywhere the absolute value 1, with 
changes of sign between 0 and 1, and between » 
and +1. With this notation we obtain 







(Pav, uisui2* * - wir) = sgn[sin((b—a—})x/n) 





k 
XII (sin(jn—a@ — 3) 4/n sin(j,—b+3)4/n) J 


h=1 





X Mabtoiirs* “Mik. (54) 





This formula can be verified by direct count of 
the sign changes caused by the operator 


Ca41Ca+2° ns Cr-1 







acting on each factor of the product 





MoMjifti2* * * Mik- 






THE BASIS OF A REDUCED REPRESENTATION 





Some Important Commutators 





The operators 





A=) ssiui=>d, Pett; B=) Cj=—-L Pu 
1 1 1 1 





are invariant against the dihedral (rotation 
reflection) group of transformations 


(j=1,2,---,m). (55) 






Mj Pej; 






From the operators P., we can form 2x linear 
combinations which exhibit the same symmetry, 
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most simply as follows: 


Au=)d, Pacten 


a=1 


shereby 
aww Ay= —B, A,=A. 


In view of (46) we have 


Anin= —CA,,= —A,C. 


From the commutators 
i Pon Pra ]= PacPra— PraP ac 


which vanish unless exactly one congruence 
between the indices is satisfied, we can form 
n—1 independent linear combinations of dihedral 
symmetry, which may be written most compactly 


Ga=} > (( Paz, Pasa z]—[P x, Poo) 


a=l1 
n 


=} } a (Poof osm, 27 Pal’ acta): 


a=1 


(S8) 


While the definition of G,, Go, will be 
extended conventionally to the complete period 
m=1, 2, ---, 2m, they satisfy the identities 


Gim=—Gn; Go=G,=0 


so that there are only m— 1 independent operators 
included in the set. The analog of (57) 


Gintn= — CGm= —GmC 


(58a) 


(59) 


is valid. The commutators of A, As, «++, Aon 
with each other necessarily can be expressed in 
terms of the set (58). By definition 


[Am, Am ]=0. 
By (57) we have also 
[Am, Amin ]=0. 
If we exclude these trivial cases, the congruences 
a+k=b+m 


will not be satisfied simultaneously when we 
compute 


a=); 


[A:, A,, |= > | Pactts P,b4m | 


a, b=1 


= i [ Poe, Pz c4m—t|+2>, | Pos Po-ath al 
c=1 


a=1 


whence 


[A ky Am ]=4Gi_m, (60) 
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With the convention (58a) this is also valid in 
the trivial case k=m. Conversely, the commu- 
tators of A;, As, «+--+, Ao», with their own commu- 
tators Gi, Ge, -+-, Gas can be expressed in 
terms of the former as follows: 


20 Gm, A, |= > ([ PesPe+m, Zz P 4,54 | 
a, b= 1 
—[ PrP ocrar Pr-x,» }). 
Again if we exclude the trivial case m=0, the 


congruences b=a, b=a+m will not be satisfied 
simultaneously. Hence, 


2[Gm, Ajot. Cl PesP esa. zy Poort | 


a=] 


+ PesPain 2 Pesnctiia | 
— [PrP s,c+m, Peso] 
—(P &Pacia, Pe-rincia))- 
Here 
[PesP ote, 2, Poors) ™ —2Pac+sPesP etm, s 


= —2P.. osal'actal cta.ett = — 2P aim, ark 


EP Pas zy Poasnctaia.! 
=2P42P sim, oP o+m,et+bim = 2P¢ o+b+m; 


etc., and we obtain 


2[Gn, A,J=> (—2Porm ett + 2P ao, c+bim 
a=1 
+ 2P 2, 0+m— 2Po—n+m, 0) 


=>" (4P, o+bim—4Po,a+k—m); 


a=l1 
or, simply 
Cc A, |=2Akim—2Ar—m (61) 


which is also valid in the trivial case m=0. It 
is easily shown from (60) and (61) that the 
commutators G;, ---, Ga_; commute with each 
other. We have generally 


[Gm, (Ax, Ao] ]=[[Gm, Ax], Ao]+[Ax, (Gm, Ao] ]. 
Using this together with (60) and (61) we find 
[Gm Ge ]=[(2Ab+m—2A nm), Ao] 
+[Ax, (2Am—2A_m) ] 
= 8Gi.4.m — 8Gk—m 


+8Gi-m—8Grim=0. (61a) 








A Set of Generating Basis Elements 


The commutation rules (60) and (61) suggest 
the introduction of the operators Xo, X1, ---, 


fm Yi, ee, | ae Zi, e's Bout defined as follows: 
2n 
X,=(4n)-! > ~cos((a—b)rx/n) Pw, 
a, b= 1 
2n 
Y,=(4n)-! > sin((a—b)rx/n) Pw, (62) 
a, b= 1 


Z,=(—1/8n) 5 sin((a—b)rx/n) 
— PR A—-P Ped. 


Again these definitions will apply conventionally 


for the entire period r=1, 2, ---, 2m, whereby 
X-=X,;, 
Y_,=-—Y,; Yo=Y,=0, (62a) 


Z,=—-Z,; Zo=Z,=0. 


Comparison of (62) with (56) and (58) yields 
directly : 
X,=(2n)— ¥ cos(mrx/n)Am, 
m=1 
Y,=—(2n)-1¥ sin(mrx/n)Am, (63a) 
m=1 


2n 


Z,= (t/2n) > sin(mrx/n)Gm. 


m=1 


Conversely, from the orthogonal properties of 
the coefficients 


A= (X,cos(mrx/n) — Y,sin(mrx/n)), 
Pr as (63b) 
Gn=—t > Z, sin(mrx/n). 


r=l 
In particular, A and B admit the expansions 
2n 
B= —Ay= —> A,-= —Xo—2X; 
1 
is —2Xn-1—- Xn, 


A=A\= > (X, cos(rx/n) — Y, sin(rx/n)) 


= Xo+2(X, cos(x/n) (63c) 
— Y; sin(x/n))+--- 
+2(X,_1 cos((m—1)x/n) 
— Y,_1 sin((m—1)x/n))—Xa. 
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Conversely, the operators (62) can be con. 
structed from A and B by means of (60), (61), 
and (63a). 

The relations (57) and (59) now become 


(1+(—)°C)X,=(1+(—)°C)Y, 

=(1+(-)'C)Z,=0 (64) 
so that the operators (62) belong to the algebra 
of even functions when r is odd, and vice versa; 
they annihilate even functions for even 7 and 
odd functions for odd r. We could deal with 
each set independently, but the two have many 
properties in common, and it is often convenient 
to keep them together. 

The commutation rules for the operators (62) 
are easily obtained from (60) and (61). Thanks 
to the orthogonal properties of the trigonometric 
coefficients, we find 

[X,, X.J=LY,, Y, |=[Z,, Z,|=0, 

[X,, YeJ=[¥r, Z.]=[2r, X.]=0; (65a) 
cos(ra/n) ~cos(sr/n). 

The only commutators which do not vanish are 


[X,, ¥Y,]= —2iZ,, 
[Y,,Z,]=—2iX,, (r=1,2,---,m—1)  (65b) 
[Z,, X,]=—2i¥,. 
We shall show next that the operators (X,, F,, 
Z,), which commute with the members of other 


sets (X,, Y,, Z,), anticommute among them- 
selves. For this purpose we compute 


2n 
(X,4+1Y, 2— (4n)-? > e(atb—e—d)rrifnP Pi, 


a, be 1 


by the definition (62). In the given expansion 
we may interchange (a, 6) and (c, d) independ- 
ently without disturbing the coefficient. Taking 
the mean of the four alternatives we obtain 


2n 
(8n)-? , e(atb—c—d)rri/n 


—_ X (PacPrat+ PaaPoet+ PrePaat+ PraPax) 
and by the exchange rule (52), this equals 


2n 
(8n)-? >. e(atb—c—d)rrij/n 


a, baw 1 «(1—(—)?) (1 —(—) 9) Pa Pra 


The only products which contribute to this sum 
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are those for which a=0 as well as c=d, whereby 
according to (46) and (47) 
Pi. = PaocPotn,ctn= 1; PacP a, c+n = PacPain,c= —C 


and we obtain 
(X,+i¥,)?= (8n)-? E Berend — (—)°C), 


a,c=1 


whence, observing (62a) 


X,?=4(1-—C) = Ro, 
Ro; (n even) 


(r=1,2,-+-,m—1). (67) 


(66) 


(X,+1Y,)?=0; 


The latter result implies that X,, Y, anticommute 


X,Y,=—Y,X, (67a) 
and that their squares are identical 
X= Y,?=R,. (67b) 


Now (65b) implies 
X,Y,= — Yx74,= —iZ,; 


and since the product in turn must anticommute 
with either factor, e.g., 


X( VX) =(X Y)X = —(YX)X, 
we have similarly 


Y,Z,= —Z,-Y,= —iX;; Z,X,= —X,Z,= —1Y,. 


These results together imply (67b) as well as 
its extension 


Z°=1X,Y,Z,=X,?= Y,2=R,. 


The operator R, thus defined satisfies the rela- 
tions appropriate to the unit of a quaternion 
algebra with the hermitian basis (R,, X,, Y,, Z,), 
viz., the real basis (R,, X,, Y,, iZ,) because 


R,X,= Y,Y,X,=1Y,2,=X,, 
etc., and 
R2=R,X,X,=XZP=R,. 


The multiplication table is accordingly 
R,=R?f=X?= Y?=Z,?, 
X,=R,X,=X,R,=tY,Z,= —iZ,Y,, 
Y,=R,Y,= Y,R,=1Z,X,= —iX,Z,, 
Z,=R,Z,=Z,R,=iX,Y,= —iY,X;. 


(68) 
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For r=0 and r=mn the operators Y, and Z, 
vanish ; in either case we obtain a commutative 
algebra with a degenerate basis (Ro, Xo) or 
(Rn, Xn); cf. (66). 

In the original quaternion algebra described 
by (31) the operators Ro, Ri, «++, Ra-1, Ra are 
projections. It remains to enumerate the dimen- 
sions of these projections [in the irreducible 
representation (30) ]. By (66), Ro and R, are 2"~ 
dimensional. By (64), the others have at most 
that many dimensions, because R, is contained 
in Ro or in (1—Rp») as r is even or odd. Further 
information can be obtained from suitable ex- 
plicit expansions of R, in terms of the basis 
elements P,,. A very compact expression is 
obtained by evaluating 


R,=}(X,?+ Y,*) 


directly from the expansions (62), which yields 


R,=2-*n-? > cos[ (a—b—c+d)rx/n |PacP oa. 

a,b,c, dat (69) 
To obtain an interesting modification of this 
formula we interchange the pairs of indices (a, }) 
and (c,d) independently. By (52) the mean of 
the four alternatives equals 


R,=2-n-? > (cos((a—b—c+d)rx/n) 


a,b, c,d 
— (—)?) cos((a—b+c¢—d)rx/n)) 
x (1 — ( — )P (eb) +D (e—4)) P, Pra. 


The anticommuting product terms in which 
either a=b but not c=d, or vice versa, cancel 
out. The terms for which neither congruence is 
satisfied yield 


2-5n-2 > sin((a—b)rx/n) sin((c—d)rx/n)PacP va. 
a,b, c,d 
The terms for which both congruences are 
satisfied yield either the identity or —(—)’C 
and add up to }(1—(—)’C). Accordingly, 
R,=(1—(—)'C) 
2n 


X(3+(8n)? > 


a, b, c, d= 1 


Xsin((c—d)rx/n)PacPea). (70) 


sin((a—b)rx/n) 


The dimensions (character) of the projection R, 
are simply enumerated from (70). When n>2, 
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we have —1+P,-Psa~1 unless a=b; c=d, so 
that the dimension of R, is one fourth of the 
total 2”. 

Dim(R,) = 2"-?; -s+,n—1); 


(n>2). (71) 


(r=1, 2, 


For the case n=2; r=1 we have from (63c) and 


(64) 


Bii+C)=—-4X1; (n=2). 


By (43a), the even eigenvectors of B are the two 
possible even homogeneous functions: (1) and 
(41me) ; the corresponding characteristic numbers 
of B are 2 and —2. Neither vanishes, so that in 
this case R; contains all of 1— Ro, and we have 
the “‘exceptional”’ result 


R,=1-—R; Dim(R:)=2; (n=2). (71a) 


Irreducible Representations 


In the previous section we have completed the 

abstract theory of the algebra generated by the 
operators A and B defined by (27) and (23), 
which are involved in the eigenwert problem 
(29). Apart from illustrative references to the 
representation (30), our main results have been 
derived from the abstract commutation rules 
(31). Now the task before us is to apply these 
results to the solution of problems such as (29), 
which can be formulated explicitly in terms of 
A and B. For this purpose we shall again think 
of the abstract numbers 5), s2, ***, Sn, Ci, Ca, 
+++, C, as linear operators in a 2"-dimensional 
vector space, defined by the matrix representa- 
tion (30). 
b. Essentially, the problem of reduction consists 
in finding the vector spaces which are invariant 
towards A and B, in the sense that neither A nor 
B operating on any member of an invariant set 
can yield anything but a linear combination of 
the vectors which belong to the same set. Since 
the set of operators (62) can be constructed from 
A and B and vice versa [cf. (63c) ], it is evident 
that the invariant sets in question must be 
invariant towards all the operators (62). In view 
of the commutation rules (65), we may specify 
for added convenience that the matrices repre- 
senting the operators (62) and their products 
shall be direct products of irreducible representa- 
tions of the factor groups (R,, X,, Y,, Z,). 


To answer the main questions dealt with jp 
the present communication it will suffice to find 
the largest solution to the problem (29). The 
identification of the corresponding irreducible 
representation of the operators (62) is an easier 
task than the complete reduction, which we 
shall leave aside for future attention. We shall 
go a little further and identify separately for the 
vector spaces of even and odd functions the 
representations of maximal dimensions. 

By (64) we may consider the sets of even and 
odd functions separately. The former are annihi- 
hilated by (R,, X,, Y,, Z,) for even r=(0, 2, 4 
+++), the latter for odd r=(1, 3, 5, +++). By 
(66), the space of even functions is formed by 
the projection 1—R,»=}3(1+C), that of odd 
functions by Ro=}3(1—C). 

In the representation (30) the operator of (29) 
has a matrix whose elements are all positive, 
Accordingly, by Frobenius’ theorem" the compo- 
nents of the “maximal” eigenvector must be all 
of one sign in this representation. It follows 
immediately that the vector in question cannot 
be odd. More precisely, it cannot be orthogonal 


to the even vector 
Xo(m1, liad, Mn) = const = * ee 


(normalized), which happens to be the only 
function of degree zero. By (43a), this is an 
eigenvector of B corresponding to the simple 
characteristic number 


(B, xo) =(B(1— Ro), xo) ="xo. (72) 
Since by (63c) and (64) 
—B=X,+Xeot- ++ +X, 
— B(1— Ro) =XitX3t- ++ +Xen-1 


. : (73) 
2X1 +2X3+::: +2X m1; 


(n= 2m) 
BO eS) le.) a 
+Xom41; (n=2m+1) 


and Xo, X;, ---, X, commute with each other, 
they commute with B. Accordingly, operation 
with X, on xo will either annihilate it or convert 
it into some function of the same degree. But % 
is the only function of degree zero. Hence, x08 
a common eigenvector of Xo, Xi, - --X, satisfying 


(X,, xo) = §-x0; (r=0, i. +++, mM). 


SO | 
nul 
we 


and 


The 


Hav 
Xi, 

whic 
oper 
satis 


then 
(X,, 


(Xs, 


so th 
mon 

t,isr 
vanis| 
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By (68) and (66) the operators X, all satisfy the 
identity 
X8—Xp=XA(X-—-1)(X-+1)=0, = (74) 
so that £0, £1, -**, &: are to be selected from the 
numbers 0, 1, —1. Since xo is an even vector 


we have . 
fo=fe=&s=---=0 


and by comparison with (73) 


fitist: ++ t+éni1=—2. 
There is just one way to satisfy this requirement 

fi=f3=---=fn1=—1, 

(X+, x0) = — Xo} 
(r=1, 3, «++, 20(m—1)/2]+1). 

Having identified one common eigenvector xo of 
X:, Xs, --*, we proceed to construct others 
which belong to the same invariant set by 


operations with Y;, Y3, ---. In general if x 
satisfies 


(75) 


(X,,x)=&x; (r=1,3,---), 
then by (68) and (65) 
(X,, (Yr, x)) =(X-¥,, x) =(— ¥-X,, x) 
= —£,(¥,, x), 
(Xs, (Vr, x)) = (Xs ¥r, x) =(VXe, x) 
=£,(Y,, x); 


so that the operation with Y, gives a new com- 
mon eigenvector; the sign of the corresponding 
t, is reversed and the others are unchanged. If é, 
vanishes, Y, will only annihilate x because by 
(68) 


(s#r) 


(¥,, x) =(V-X7?, x) = (Vr, & x). 


For this reason it is useless to operate with 
Yo, Y4, «++ on xo. It is equally useless to operate 
more than once with any one Y, because by (68) 
Y, is its own inverse in the subspace R, 


Y?7=X/=R,, 


so that the second operation at best undoes the 
effect of the first. Finally, since 
Z,=1X,Y, 


the operation with Z, can only give us such 
eigenfunctions of X, as we can obtain with Y, 
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alone. Accordingly, the vectors which we can 
obtain by operating with various combinations 
of (Yi, Ys, ---) on xoe—at most once apiece— 
form an invariant set. The vectors thus obtained 
can be distinguished from each other by their 
different sets of characteristic numbers; when 
necessary we shall designate them individually 


x(&1, &s, +++, Exm—13 (u)); 


In the case of odd n we have §,= —1 for every 
vector in the set; there is no Y,. For the rest, 
every combination occurs exactly once, as is 
expressed by the following formulation of the 
relations derived above 


(X,, x(E1, £3, +++, E2m—13 (u))) 
= E-x(£1, &3, +++, Eam—15 (m)), 
(Y,, x(&1, &3, +++, Eem—13 (u))) 
=x(£i, is Sa 
(Xn, x(E1, Es, +++, Env} (u))) 
= —x(f1, +++, En—2; (u)). 


A more explicit description of the functions x is 
not necessary. The relations (76) in conjunction 
with the verified construction of the one vector 


x(—1, =1, +++, = 1; (w)=2-"2 


(n=2m=n—1). 


(76) 
"ys Eom—13 (u)), 


(76a) 


suffice for definition. From the one given vector, 
the construction is in effect completed by means 
of the formal relations (68). 

The representation (76) is analogous to (30) ; 
but it has only m=[m/2] factors and the number 
of dimensions is 2!"/*], if we denote by [x] the 
greatest integer which does not exceed x. More- 
over, denoting by Q, a general function of the 
factor basis (R,, X,, Y,, Z,), the representation 
(76) is a direct product of the two-dimensional 
representations 


Dymax = D2(Q1) X D2(Qs) X + - - XK Do(Qn-1); 
(n= 2m) 
Drax = D2(Q1) X D2(Q3) X - - - 
X D2(Q,-2) X D_(Q;); 


—1 0 0 1 
Dix.=(; ‘): Da(¥.)=(, i 


Di(X,)=(1); D_(X,)=(—1) 


(n=2m+1) 
(77) 














supplemented by the one-dimensional repre- 
sentation D_(Q,) when 1 is odd. 

An invariant vector space may be constructed 
by the same procedure from any common 
eigenvector of (Xo, X1, --+, X,). The characters 
of the representation will be determined by the 
initial & and &, which are common to all 
members of an invariant set, together with the 


values (0 or 1) of &:*, &’, ---, ° As For the signs 
of £1, 2, «++, &.-1 can be changed independently 
by operation with Y;, Ys, ---, Yn-i, excepting 
those which vanish in the first place. To put it 
in geometrical language, the general type of the 
representation is determined by an intersection 
of projections 


PoP,:+-Pra 


where for r=1, 2, ---, n—1, P, is either R, or 
1—R,, while Ro is subdivided into }(Ro+Xo) 
and 3(Ro—Xo), and R, is similarly subdivided. 
The vector space constructed above from xo is 
the only intersection common to all of the 
projections R;, Rs, ---, Rtnj2}-1. For every such 
intersection must contain an even eigenvector 
of B corresponding to the characteristic number 
n when n is even or n—1—£, when 2 is odd, in 
either case an even function of degree less than 
2, and the only function that satisfies this 
specification is const. = xo. 

We have already pointed out that by Fro- 
benius’ theorem, that eigenvector which corre- 
sponds to the largest solution of (29) cannot be 
orthogonal to xo. Therefore, this eigenvector 
belongs to the invariant set constructed above, 
which contains xo. The expectation of a simple 
largest characteristic number is in accord with 
our finding that the invariant vector space which 
corresponds to the maximal representation (77) 
is unique. 

In connection with the considerations of 
Lassettre and Howe," it ia a matter of interest 
to identify not only the largest solution Amax of 


(29), but in addition the largest solution eda 
which belongs to an odd eigenvector. As it 
happens, this solution is the second largest in 
magnitude. We shall prefer a slight variation of 
the above procedure. From one common eigen- 


13 FE, N. Lassettre and J. P. Howe, reference 12. 
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vector of the operators 

r*= —X,cos(ra/n)+ Y, sin(rx/n); 
(r=0,1,---,m—1,m) (78) 
X,*=X, 







Xo*= —Xo; 


we shall construct others by operations with, 
say, Z1, -++, Zn-1, which anticommute with the 
corresponding X,*. Then instead of B we consider 









2n n—l 


A =>" SSjii= -> X*=Xo-> 2X,*—X,, 
j=1 


r=1 t=1 






RA =}(1—C)A (79) 
_ [Xo—2XF— AK F— ++ —2Xwa—X nj (n= Im) 
~ X—2XF—-2XF—---— 2X; (w= I2m41) 
(cf. (63c)]. By (68) the operators (78) satisfy 
X,"=R,=X?; X,%-X*=0. (80) 








The matrix which represents the operator 





exp(H’A) =] (cosh H’+5j5 5,1 sinh H’) 
jel 





in terms of the vector basis (41) which renders 
B diagonal consists of two 2*-'-dimensional 
matrices which transform even functions among 
themselves and odd functions among themselves, 
respectively. In either matrix all elements are 
positive (provided H’>0), so that Frobenius’ 
theorem applies in the vector spaces of even and 
odd functions separately. The modification by 
the factor exp(H*B) in (29) does not affect this 
conclusion because this factor is represented here 
by a diagonal matrix whose elements are all 
positive. 

By (44) the largest characteristic number n of 
A is double and belongs to the two configurations 
(+, +, -°°, +), (—-, —, °**, —), which are 
free from alternations of sign. When the corre- 
sponding even and odd eigenvectors 















xo" 


\- ae me ll (1+4,) 


xo* 
2° II (1—y;)) (81) 






are expanded in terms of the vector basis (41), 
their components are all positive, so that the 
representations which correspond to the charac- 









teristic numbers Amax and Amax of (29) must 
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belong to unique invariant vector spaces which 
contain xo* and x0, respectively. The former 
has been constructed above from xo. The odd 
vector xo is a common eigenvector of all X,* 
and satisfies 


nxo—* = (A, xo") 


= —> (X2q, x0") = — 


q=l 


7 + 
>> ei xo" 
q=1 


which requires 


o= —&*=1, 


te=1; (0<q<4n) 
t,=§,*= =1; 


so that the vector xo“* is included in the 
projection 


(n=2m), 


3(Rot+Xo0)ReRa- + + Ropnjay 


and in 3(R,—X,) when 1 is even. The ‘‘maximal 
odd” representation is therefore 


Dis = D,.(Qo) xX D2(Q2) > See X D2(Qn-1); 
(n=2m+1) 
Dmnax = Ds(Qo) X D2(Q2) X - + - 
X D2(Qn-2) XD_(Qn); ("= 2m) 


with the notation of (77). When is odd, Dinax 
has 2-/2 dimensions, the same number as 
Dax, and there can be no further representation 
with this number of dimensions, for there is only 
one pair of configurations with less than two 
alternations of sign. For even nm the representa- 


(82) 


tion Dax has only 2\-*)/? dimensions; it is still 
unique. However, one other representation with 
the same number of dimensions is compatible 
with the known spectrum (44) of A, namely, 


Dai = D_(Qo) X D2(Qz) X- - - 
x D2(Qn-2) Xx D.(Qn) ° 


The count of configurations with m (even) alter- 
nations yields one pair: (+ — + — --- + —) 
and (— + — +--+ — +); reasoning analogous 
to the above shows that exactly one invariant 
vector space belongs to (82a). 

The identification of the remaining invariant 
vector spaces involves more elaborate calcula- 
tions; the complete theory will be given in a 
later publication. However, in order to settle an 


(82a) 
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interesting question which will arise, we shall 
complete the discussion of those vector spaces 
which contain the m functions of degree one 


x1(27; 1, ***, Mn) =n) > exp(2rkwi/n)u,. (83a) 
k=l 

These functions may be classified according to 
the representations of the dihedral symmetry 
group (55). To the identical representation of 
(55) belongs x:(0; (u)) ; the corresponding repre- 
sentation of the algebra (68) is (82). For even n, 
xi("; (u)) belongs to another one-dimensional 
representation of (55) and to the representation 
(82a) of (68). The remaining functions (83a) 
belong pairwise to two-dimensional representa- 
tions of (55). To show that the pair of similar 
representations which belong to the pair 
x1(+2r; (u)) are described by 


D_(Qo) X Do(Q2) XI] D2(Qes); 


sr 


D_(Qo) X D_(Q,.) X Do(Q2-) X J] D2(Qz.); 


sar 


(m odd) (83b) 


(nm even) 


one may verify that the vectors in question are 
annihilated by Xe, 


(X er, x1(2r; (u))) = (Xery x1(— 27; (u))) =0; 


(0O<2r<n). (83c) 


The direct calculation by means of (62) and (54) 
is not too laborious; since X2, commutes with B, 
only the terms of first degree have to be com- 
puted. The remaining factors in (83b) are then 
determined as before by the relation 


(B, x1) = (n—2)x1. 
SOLUTION OF THE EIGENWERT PROBLEM 


It is a simple matter now to compute those 
solutions of the eigenwert problem (29) which 
belong to the representations (77), (82), and 
(83); we only have to solve a set of quadratic 
equations. 

As to the “physical” significance of the solu- 
tions, we recall that the characteristic number 
itself yields the partition function and with it 
the thermodynamic functions of the crystal. 
Moreover, as shown by Kramers and Wannier, 
the individual components (wu, «++, un) of the 
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eigenvector which belongs to Amax describe the 
statistical distribution of configurations of the last 
tier of atoms added to the crystal." On the 
other hand, if we know in addition the eigen- 
vector ¥’(u1, «++, un) which belongs to Amax for 
the adjoint operator V’= V,'V2’= V, V2, then the 
component products 


W(u1, -**, MnO’ (ua, + * +, Mn) 


describe the statistical distribution of configura- 
tions in the interior of the crystal. 

With the latter result in mind, we first trans- 
form V to self-adjoint form by operation with 
V;!; then 


V= Vi} VV, t= Vi} Ve V;! 
= (2 sinh 2H)” 
Xexp(3H’A) exp(H*B)exp(3H7’A) (84) 


has the same characteristic numbers as V, and 
the squares 


W (ur, -° +, Mn) =f (ur, - ++, Mn) 


of the components of its principal eigenvector 
will describe the distribution of configurations 
in the interior of the crystal. 

To solve the eigenwert problem we may 
transform V into a function of the operators 
Xo, X1, «++, Xn alone, or equally well into a 
function of Xo*, X,*, ---, X,*. The latter 
procedure involves slightly simpler relations; 
either result is readily converted into the other 
by means of the real orthogonal operator"! 


exp(= ((n=1)/2n)xiZ, ) 


r= 


=" (1—(1—sin(rx/2n))R,+cos(rx/2n)iZ,). 


r= 





For reference we observe here the somewhat 
more general rule 


exp(iaZ,)(X, cos 8+ Y, sin 8) exp(—iaZ,) 
=X, cos(2a+6)+ Y,sin(2a+ 8) (85) 


which is easily verified by means of the multi- 
plication table (68). The latter—and (85)—are 


“ This is not exactly the dual transformation; it takes 
a further transformation with II(1—R,+X,) to obtain 
(86). 
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also valid if X,, Y,,Z, are replaced throughout by 
X,*= —X, cos(rr/n)+ Y, sin(rx/n), 
Y,* =X, sin(rx/n)+ Y, cos(rx/n), (86) 
Z,*=—Z, 

(dual transformation). We now recall the 

expansions 

A= —X o*—2Xi*—---—-2X,1—-X," 

B= —Xo—2X1—- ++ —2Xn-1- Xn 


n—l 


= Xo*+2 > (cos(ra/n)X;,* 
= —sin(rx/n) Y,*) —X,.*. 
Substitution in (84) yields 
V =(2 sinh 2H)" exp((H* — H’)Xo*) 


n—l 
II (U,) exp(—(H*+H")X,.*) (87) 


r=1 


where 
U,=exp(— H’X,*) exp(—2H*X,) exp(— H’X,*) 
=(1—R,)+(cosh 2H’ cosh 2H* 

—sinh 277’ sinh 2H* cos(rx/n))R, 

— (sinh 2H’ cosh 2/]* 

—cosh 2H’ sinh 2H* cos(rx/n))X,* 

—sinh 2H* sin(rr/n) Y,*. (87a) 

By inspection of the factored form we note that 
transformation with (1—R,+7Z,) changes U, 
into U,-', which differs from U, in the signs of 


the coefficients of X,* and Y,*. Adding the two 
together we obtain 


U,+ U-"=2(1—R,-+(cosh 2H’ cosh 2//* 
—sinh 2H’ sinh 2H* cos(rx/n))R,) 
=2(1—R,)+2 cosh y,R, 
whence the characteristic numbers of U, be- 
longing to the subspace R, are e’ and e~*, of 
equal multiplicity. (Those which belong to 
(1—R,) are trivial, all equal to unity.) Accord- 
ingiy we may write 
U,=1—R,+ cosh y,R-—sinh y, cos 6,*X,* 
—sinh y, sin 6,* Y,* 


=exp(—vy,-(cos 6,*X,*+sin 6,* Y,*)) (88) 


where y, and 6,* can be computed from 2H’, 
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2H*, w-=rx/n by the rules of hyperbolic trigo- 
nometry, as indicated by Fig. 4. We note in 
particular the trigonometric relations'® 
cosh y=cosh 2H’ cosh 2H* 
—sinh 2H’ sinh 2H* cos w, (a) 
sinh y cos 6* =sinh 2H’ cosh 2H* 
—cosh 2H’ sinh 2H* cos w, (b) 
sin w/sinh y=sin 6*/sinh 2H* (89) 
=sin 6’/sinh 2H’, (c) 
cot 6*= (sinh 2H’ coth 2H* 
—cosh 2H’ cos w)/sin w. (d) 


With the aid of (85) and (86) we find immediately 
the transformations which reduce U, to diagonal 
form with regard to the eigenvectors of either 
X,* or X,; as a vector basis, 


U,=exp(— }6,*iZ,) exp(—v-X,*) exp(46,*7Z,), 
(90a) 
U,=exp(3(% — 


—5,*)iZ,) exp(—y-X,) exp(— 3(4 —w,— 6,*)iZ,), 
(r=1, 2, ---,m—1) (90b) 


and by substitution in (87) 


n—l n—l 
exp(X 48-2, V exp(- pm 18,%1Z,) 


r=1 r=1 


=(2 sinh 2H)" exp( (H*—1')Xe" 


n—l 
—> 7-X,*- (H'+H")X,"), (91a) 


r=1 


n—l 
wove 


r=1 
n—l 

- i)iZ,) V exp(X 3 (4 —w,— i)iz,) 
r=1 

=(2 sinh 2H)*"/? exp((H’-H1")X, 


n—l 


= elt (H'+H")X,). (91b) 


* The formulas of hyperbolic trigonometry are obtained 
from those of spherical trigonometry by the substitution 
of imaginary lengths for the sides. See F. Klein, Vorlesungen 
~~ eet Geometrie (Springer, Berlin, 1928), 
p. 195. 
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a D* 


2H’ 








2H* ~o 

Fic. 4. Hyperbolic triangle. Stereographic projection, 
conformal. Circles are represented by circles, geodetics by 
circles invariant towards inversion in the limiting circle 


Ce of the projection. See F. Klein, reference 15, pp. 
293-299. 


The spectrum of V is now known to the extent 
that the common eigenvectors x of (Xo, Xi, ---, 
X,n-1, X,) and the corresponding sets of char- 
acteristic numbers (&, £1, «++, &—1, &) are 
known. For every solution x to that problem V 
has an eigenvector y with characteristic number 
\ given by 


n—l 
y= (exp ap > ir—w,—4*)iZ,), x): 


T= 


(V,¥) =v, (92) 
log \= $n log(2 sinh 2H) + (H’ —H*)& 


n—l 
= ,¥ Vrér— (H’+H*)é,. 


r=1 


In the previous section we were able to show that 
the largest characteristic number of V belongs 
to the representation (77). In view of (76) and 
(76a) we must put {== =---=0. All the 
others may be taken negative, which calls for 
x=xo, and obviously this is the best choice. 
Accordingly, 
log Amax — $” log(2 sinh 2H) 
oe ‘+++ Yam-1; ("=2m) (93) 
Yityst ++ +Y2m-1+(H’+H") ; (n=2m+1). 
These two results can be combined if we adopt 
the natural conventions 


Y-r=r3 Yn =2H'+2H"*. (94) 
Then (93) takes the compact form 
log Amax — 3” log(2 sinh 2H) =4 ¥ ye,-1 
r= 


=} > cosh (cosh 2H’ cosh 2H* 
ral 


—sinh 2H’ sinh 2H* cos((2r—1)#/2n)). (95) 










































Fic. 5. Variation of 6* with w for the quadratic crystal, 


at various reduced temperatures (&//)7T=1/H. 


gd 2H 1/H 
0° oo 
15° 7.551 
30° — 3.641 
42° 2.472 
45° 2.269 
48° 2.089 
66°T 1.29 
77°t 0.92 


+ By mistake marked 60° on figure. 
¢ By mistake marked 75° on figure. 


The corresponding eigenvector can be described 
by either of the two formulas 


Vmax = (exp = > 4 (4 —wor-1 
0<2r<n 
= ii)iZn-), ), (96a) 
Wmax =f exp pa $05-1iZs-1) ’ w') . (96b) 
0<2r<n 
Here we have omitted operations with 


exp((const)Z:,); (r=1, 2, ---), because these 
have no effect on xo or any other even vector. 
The somewhat simpler form (96b) involves con- 
struction from the even 6 function xo* given by 
(81); one easily verifies that £,*, §;*, --- are all 
negative for this vector. 


The Principal Eigenvector 


The result (96) describes explicitly the distri- 
bution of configurations (¥”) of one row of atoms 
in a crystal. The terms of the description are 
quite unfamiliar. We should not regret this; for 
among the main objects of the development of 
this theory is the invention of more suitable 
methods for the description of such distributions. 
Even so, we must try to establish the connection 
with more familiar types of description in terms 
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of more tangible quantities. Moreover, the result 
(96) fails to describe the statistical correlation of 
configurations in different rows of atoms. Meth. 
ods of dealing with these two tasks have beep 
perfected to a considerable degree; but the 
algebraic apparatus involved is elaborate enough 
to make a separate publication advisable. 

One quite tangible result is readily obtained 
from (96b): The probability of a configuration 
free from alternations is 


2y(++-+++)=2¥(——----) 
= II (cos*($82,.)), 


Similarly, when 7 is even so that a configuration 
free from persistencies of sign is possible, the 
probability of such a configuration is 


2Y3(+ — + — + + —) =2p'(—+—4+---—+4) 
= J] (sin?(}82,_,)), 


Presumably, the probabilities of the 2"—4 re. 
maining configurations are intermediate between 
the given limits. One can show from the commu. 
tation rules that the probabilities of configura. 
tions with 2, 4, 6, --- alternations of sign can be 
computed successively from those terms in the 
expansion of (96b) which contain, respectively, 
i oo factors of the type (Z, sin }4,*), 
However, the labor involved increases rapidly 
with the number of steps. 

Of greater interest is the fact that the descrip. 
tion (96) brings out a striking qualitative differ. 
ence between those distributions which occur for 
n(H'—H*)>1 and those which occur for 
n(H*—H’')>1. 

The variation of 6* with w=rx/n is given 
analytically by (89d); but a qualitative inspec- 
tion of Fig. 4 is even more convenient. Let us 
keep the side OD*=2H’ of the triangle OD*D 
fixed; then the positions of the vertex D’ for 
w=n/n, w=3nr/n, +--+ are equidistant points on 
a circle of radius OD’ =2H*. Now if 2H*>2H,, 
then 6*=-2 for w=0, and 6* decreases from f to 
0 as w increases from 0 to x. On the other hand, 
if 2H* <2H’, then 6*=0 for w=0. As w increases, 
6* increases to a maximum value 


sin-"(sinh 2H*/sinh 2H") = bmax <4 





































given by (89c) for 6’=42; with further increase 
of w, 6* decreases again until 6*=0 for w=. 
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As regards the angle 3(*—46*—w) which occurs 
in (96a), the axiom 


i'+6*+w<r 
of hyperbolic geometry implies 
x—w—d*>0. 


A graphical representation of 5* as a function of 
w at various temperatures is given in Fig. 5 for 
a quadratic crystal (H’=H). As regards the 
variation of 6* with the temperature for a fixed 
w, we note that as the former increases from 0 
to «, 6* increases from 0 to r—w. For small 
values of w, nearly all of this increase takes place 
over a small interval of temperatures which includes 
the critical point. A detailed inspection of the 
operators Z, can be made to show that those for 
which r is small are the most effective in intro- 
ducing pairs of alternations far apart, which is 
just what it takes to convert order into disorder. 

Concerning the extreme cases of low and high 
temperatures, one sees easily that in the former 
case 6* and in the latter case r—w—6* will 
remain small for all w. This is neither new nor 
surprising; it means that the low temperature 
distribution consists mostly of configurations 
which resemble the perfectly ordered arrange- 
ment represented by xo*, while the high tempera- 
ture distribution resembles the completely ran- 
dom distribution described by xo. 


Propagation of Order 


The principal characteristic number (95) of 
the eigenwert problem (29) yields the partition 
function, and the principal eigenvector deter- 
mines the distribution of configurations. Never- 
theless, the remaining 2"—1 characteristic num- 
bers are of some interest. In general, these 
describe the propagation of order.'* In particular, 
if long-range order is present, at least one of the 
subsidiary characteristic numbers should be very 
nearly equal to the principal.'? Moreover, at the 
transition point where long-range order appears, 
one may expect a phenomenon analogous to the 
branching of multiple-valued analytic functions 
such as 2'/", whereby the order of the branch- 

'® F. Zernike, Physica 7, 565 (1938). 


'? See reference 13. See also E. Montroll, J. Chem. Phys. 
10, 61 (1942). 
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point bears a relation to the type of the singu- 
larity. A survey of the solutions which belong to 
the representations (77), (82), and (83) will 
answer the most important questions. 

To obtain all the solutions which belong to 
(77b), the angles 46,*, $6;*, --- in (96) may be 
replaced independently by 4$6,;*+42, 463*+42, 

++; the corresponding characteristic numbers 
are given by the formula 


log \— $n log(2 sinh 2/7) 
nen °° *+Vem—13 
pyibyst ++ Lb Yem—1 +37 2m413 


(n= 2m) (97) 
(n= 2m-+1) 


where the optional signs are independent and 
each combination occurs once. 

That part of the spectrum of V which belongs 
to the representation (82) is obtained similarly. 
The result is 


log \— $n log(2 sinh 2H) 


(H’ —H*) +yetyat +++ bY 2m-2+4 72m} 
n=2m 
= ( (98) 
(HI! —H*)Ayotyst +++ 472m} 


(n=2m-+1). 


For the largest of these we find the counterpart 
of (95) 


log AS. — 4 log(2 sinh 2H) 
n—l 


=(H'—H*) +32 72= 


r= 


= (x sgn(H’—H*) +5. v») (99) 


rel 
with the natural notation 


yo=cosh~'(cosh 2H’ cosh 2H* 


—sinh 2H’ sinh 2H*)=2|H’—H*|. (100) 


First of all, this differs from (95) in the substitu- 
tion of yo, 2, -** for v1, Ys, «+. Moreover, the 
sign of }7o is mandatory and changes at the 
critical point. This leads to a most remarkable 
limiting result for large values of m. The two 
sums 


Yityat +++ +Y2n-1=7(4/n)+7(3e/n)+---, 
Yor V2t + +Y2n-2= (0) +7(20/n)+--- 
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may be considered as different numerical quad- 
ratures of 


(n/2n) f ‘Velde 


by the trapeze rule. For periodic, analytic func- 
tions the approximation to the integral improves 
very rapidly (exponentially) with increasing 
number of intervals. Hence, with rapid con- 
vergence 


{ (H'>H*) - 


lim (Assax/Xmax) = (1 
exp(2H* — 2H’); (H'<H*). 


The distribution of the logarithms of the re- 
maining characteristic numbers which belong to 
(77) and (82) 


log Amax — 271, log Amax 

— 273, vey log Amax — 2¥1— 273, tee, 
log a log ed 

—2y4, +++, log Nmax— 272-274 °° 


becomes dense as n—«. These ‘“‘continua’’ re- 
main distinct from the two largest characteristic 
numbers as long as H’#H*. However, for 
H’ <H* the representations (83b) give rise to a 
continuum 


log A=log Amaxt+-Y0—Y2r; (H’<H*) (102) 


which contains it as a superior limit. Further 
complications due to representations which we 
have not investigated here are excluded by (43) 
and (63c). 

At the critical point H'=H*, we have yo=0; 
yr=O(r/n). In this exceptional case both “a 
and Amax itself are limits of ‘‘continua.”” The fol- 
lowing scheme summarizes the results 


Naz =A. > (continuum); (H’>H"*) 
Amax=Anax=lim (continuum); (H’=H*) (103) 
Nmax > Anse =lim (continuum); (/7’</#*). 


At temperatures below the critical point 
(H’>H*), we have ‘‘asymptotic degeneracy”’ of 
order 2, as a symptom of long-range order. 
Above the critical point, there is no degeneracy 
of the principal characteristic number. At the 
critical point we observe branching associated 
with an “asymptotic degeneracy” of infinite 


order. 


In regard to the ‘“‘propagation of order,” the 
mean distance to which a local disturbance in the 
crystal is propagated is inversely proportional] 
to log(Amax/A). We note that this ‘‘range”’ of the 
“short-range order’’ becomes infinite at the 
critical point. 

Since the functions (83a) belong to different 
representations of the dihedral group (55), —and 
the same is necessarily true fer their respective 
invariant vector spaces—, the result (102) ought 
to contain information about the anisotropy of 
the propagation of order. A tentative computa- 
tion leads to the implicit formula 


cosh 2H cosh 2H’—sinh 2H cosh(@ sin ¢) 
—sinh 2H’ cosh(8 cos ¢g)=0 (104) 


for the mean range (1/8) of the short-range 
correlation in the direction ¢, when H’ <H*. 


THERMODYNAMIC PROPERTIES OF A 
LARGE CRYSTAL 


To compute the partition function per atom 


A=Ax =lim(Amax)"” (105) 


n=O 


for an infinite crystal we replace the sum (95) by 
the integral 


1 r 
log A. =} log(2 sinh 2H) +—- f y(w)dw (106) 
rT 


0 


where 
cosh y(w) =cosh 2H’ cosh 2H* 
—sinh 2H’ sinh 2H* cos w. 


There are several ways to show that (106) 
actually describes a symmetrical function of H 
and H’. For example, with the aid of the useful 
identity 


9 


f log(2 cosh x—2 cos w)dw=2rx (107) 
0 
we can convert (106) into the double integral 


log(A/2) = ef f log(cosh 2H cosh 2H’ 
0 0 


—sinh 2H cos w—sinh 2H’ cos w’)dwdw’. (108) 


It seems rather likely that this result could be 
derived from direct algebraic and _ topological 
considerations without recourse to the operator 


eas 
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method used in the present work. Such a de- 
velopment might well amount to a great im- 
provement of the theory. 

A generalization of the expansion whose initial 
terms were given by Kramers and Wannier is 
easily obtained from (108). With the notation 

2x=tanh 2H/cosh 2H’ =sin g cos g’, 


; (109a) 
2x’ =tanh 2H’/cosh 2H=cos g sin g’ 


we expand the logarithm in powers of « and x’ 
and integrate term by term to obtain 


log A—} log(4 cosh 2H cosh 2H’) 


iff logit — 2x cos w— 2k’ cos w’)dwdw’ 
0 


=—4 > (2r+2s—1)!(r!)—*(s!)—*x?"x**. (109b) 


r+s>0 


Specialization to the case H=H’; x=x«’ of 
quadratic symmetry yields 


log \—log(2 cosh 2/7) 


= ef f log(1 — 4x cos w1 COS w2)dwidw. 


0 


5 (()'/)e 


=log(1 —x* —4«* — 29x® — 265x®— 2745" 
— 30773! —364315x!4—---) 


(109c) 


which confirms the result given previously to the 
order «x by Kramers and Wannier.”® 

These expansions converge for all values of H 
and H’ because 


| 2x cos w+ 2k’ cos w’| =2«+2«’=sin(g+g’) =1. 


The limit of convergence of the series is given by 
the transition point,!* and it converges even at 
the limit. Of course it is not very suitable for 
computation in the critical region; much better 
formulas for this purpose will be obtained. 

The partition function (106) yields directly 
the free energy F of the crystal; the energy U 
and the specific heat C are obtained by differ- 
entiation with regard to the temperature 7. 


* Kramers and Wannier might well have inferred this 
much from the uniform sign of the terms, which locates 
the nearest singularity on the real axis. 


For a crystal of N atoms 
F=U—TS=-—WNkT log i, 
U=F—T(dF/dT) = NkT°d(log \)/dT, 
C=dU/dT. 


(110) 


For the purpose of differentiation it will be 
convenient to consider A, given by (106), as a 
function of two independent variables H = J/kT 
and H’=J'/kT. In this notation we have 


U=—NJ(@ log }/8H) — NJ'(a log \/dH’) 
= — NkT(H(@ log \/dH) 
+ H’(d log \/dH’)) 
C= Nk(H?(@ log \/dH”) 
+2HH'(é log \/dHdH") 
+H"(# log \/dH"”)). 


(111a) 


(111b) 


The two terms of (111a) are separately the mean 
energies of interaction in the two perpendicular 
directions in the crystal (remember the interpre- 
tation of temperature as a statistical parameter). 

In differentiating the integral of (106) we 
consider H* as a function of H which satisfies 
[ef. (17) ] 

dH*/dH = —sinh 2H* = —1/sinh 2H. 


The following formulas, wherein y, 5’, and 6* are 
considered as functions of H’, H*, and w, are 
easily obtained by differentiation from the 
formulas (89); some of them are obvious by 
inspection of Fig. 4. 


d7/0H' =2 cos &*, 
d7/dH* =2 cos 6’, 
0°y/dH" =4 sin*® 6* coth y, 
0°y/dH* =4 sin? 8’ coth y, 
0°y/dH'dH* = —4 sin 5* sin 8’/sinh y. 


(112) 


Comparison with (106) yields for substitution in 
(111) 
8 log \/dH'’ = | cos 6* dw/r, 

ee (113a) 


8 log \/8H =cosh 2H* —sinh 2H f cos 6’ dw/m 


0 
and 


2 log \/dH”? = 2f sin? 6* coth y dw/r, 


0 
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3? log \/0 HAH’ =2 sinh 2H* 


xf (sin 6* sin 6’/sinh y)dw/r, 
0 


113b 
3 log \/0H*=2 sinh? 2H* \ 


4 (-1+c0th 2H f cos 5’ dw/r 
0 


+f sin? 5’ coth y du/r). 
0 


The reduction of these integrals by means of an 
uniformizing elliptic substitution for the hyper- 
bolic triangle of Fig. 4 is dealt with in the 
appendix. The qualitative behavior of the inte- 
grals is easily seen from Fig. 4: The integrals 
(113a) are continuous functions of H’ and H 
(or H*) for all values of these parameters, even 
for H’=H* (critical point). The three integrals 
(113b) are infinite at the critical point, otherwise 
finite. The singularity results from a conspiracy : 
In the case H’=H* we have y(0)=0 and at the 
same time 6'(0)=6*(0)=42, although at all 
other temperatures 
¥(w)27(0) =2| H’—H*| >0, 
sin 6’(0) =sin 6*(0) =sin 0=sin r=0. 

For the special case of quadratic symmetry 
H'=H the computation of the thermodynamic 
functions can be simplified considerably. The 
most convenient starting point is the double 
integral of (109c), which is converted into a 
single integral with the aid of (107). Using the 
notation 

k,=4x=2 sinh 2H/cosh? 2H, (114) 
ky” =+(1—k,?)!=2 tanh? 2H—1; [ky’’| =k,’ 


we obtain 


log(A/2 cosh 27) 


== [toga + 1a sin? g)'))dg. (115) 
2r Jo 


Differentiation under the integral sign yields for 
the energy U, 








2 
—NJ coth 21( 1+"i1"K:) (116) 
T 


LARS ONSAGER 


and for the specific heat C, 


d* log 


C= Nk 7 = Nk(H coth 2H)*(2/7) 





FY? 

X (2K1—2E:—(1—hki")(34+hi1"Ky)). (117) 
In these formulas K, and £; denote the complete 
elliptic integrals 

x/2 
K.=K(ki)= [ (1-ki'sin? ¢)-Mde, 
“70 
" (118a) 
E,= E(k) = fi (1—k;* sin? ¢) ‘doe. 
0 
The integral (115) cannot be expressed in closed 


form; but rapidly convergent series can be 
given. With the notation 


K,'=K(ky’), 
log Q=mryt= —7K,'/K,, (118b) 


G=1-*—3-24+-5-*—7-24.... 
= 0.915 965 594 (Catalan’s constant), 


one or the other of the following expansions (de- 
rived in the appendix) will be found suitable for 


computation: 


log \= 3 log(2 sinh 2H) — } log qi 


+ (—)"(2n—1) log(i—qi"“"), (119a) 


n= 1 


1 « 
log \= 4 log(2 sinh 2H) +(2/7)G+— > (—)" 


T n=(0 


(1+(2n+1)(wi/71) —exp[ — (2n+1)xi/7;}) 


(2n+1)? sinh*((m+4)i/71) 





119b) 


A singularity occurs for H=H*= } log cot 1/8, 
in which case k}x=1; Ki= ©; K,’=$; E,=1. 
The specific heat becomes infinite at this critical 
point; the energy is continuous because k,’’=0. 
The analytic nature of the singularity is evident 
from the approximate formulas 


1~log(4/k1’) ~log(28/| H—H-|), 
K,~log(4/k1’') ~log(2'/ | I) (120) 


C/Nk~(2/x) (log cét #/8)?(Ki—1—jr). 





oO - &@® 3 





CRYSTAL STATISTICS 


The critical data are 
H.=J/kT- =} log cot r/8 = 0.440 686 8, 
— F./ NkT-=log \.= } log 2+(2/x)G 
=0.929 695 40, 
— U,/NJ=2'=1.414 213 56, 
S./Nk=\log \-— 2! H.=0.306 470 
=log 1.358 621. 


(121) 


The critical temperature and energy were given 
already by Kramers and Wannier." In addition, 
their estimate of “‘about 2.5335” for \. was very 
close to the exact value 


Ae = 24e?/* = 2.533 737 28. 


The specific heat of a quadratic crystal is 
represented in Fig. 6 as a function of 1/H 
=(k/J)T; the result of the best approximate 
computation in the literature’ is indicated for 
comparison. 

Figure 7 represents the specific heat curve of a 
highly anisotropic crystal (J’/J=1/100), com- 
puted from Eq. (6.2) in the Appendix. The 
corresponding curves for the case of quadratic 
symmetry (J’/J=1) and for the linear chain 
(J'/J=0) are also shown. The abscissa is 
2/(H+H’')=2kT/(J+J’), so that the three 


curves have the same area. 


SOME PROPERTIES OF FINITE CRYSTALS 


The partition function of any finite crystal is 
an analytic function of the temperature. Perfect 
order and a sharp transition point occur only 
when the extent of the crystal is infinite in two 
directions. It is a matter of interest to study the 
approach to complete order and the sharpening 
of the transition point with increasing size 
(width) of the crystal. 


Boundary Tension 


First, let us compute the “‘boundary tension” 
between two regions of opposite order below the 
transition point. This we can do by a slight 
variation of the model. We simply chose a 
negative interaction energy between adjacent 
atoms in the same tier: 


n 
u=+J' Do ujmisr- 
1 
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Now if m is even, the reversal of J’ has the same 
trivial effect as a redefinition of the sign of 
every other y;. The interpretation of the result 
is a little different in that the type of order 
which occurs below the transition point is now a 
superlattice with the ideal structure indicated 
by Fig. 8. However, let us choose an odd value 
of n. Now a perfect alternation of signs around 
the polygon of m atoms is impossible; in one 
place at least the given superlattice must be 
adjacent to one of opposite order. As we build 
a long crystal on the odd polygonal base, we 
build it with one misfit ‘‘seam.”’ The free energy 
difference due to this seam will equal that of a 
boundary between regions of opposite order. 
The formula (93) is still valid, only now 


J'/kT = —H'>H*>0 
whence 
log Amax — 3” log(2 sinh 277) 
= t+73st +++ +Y2m-1— (| H’| —H*) 
= Vit y3st +++ + Y2m—1+ FY 2m41— 2(| HM’ | —H*) 
where ie 


cosh y,=cosh 2H’ cosh 2H* 
—sinh(—2H”’) sinh 2H* cos(x—w,). 


The replacement of w, by r—w, has no appreci- 
able effect on the sum; the modification of the 
partition function due to the seam is therefore 
practically equal to 


—2|H’|+2H*. 


This represents the effect of a lengthwise 
boundary for our original model as well, and 
we obtain 


a’ /kT = —log ds’ = 2(H’—H*), 
o’ =2J'—kT log coth(J/kT) 


(122a) 


for the boundary tension o’ (per atom) of a 
longitudinal boundary. The transverse boundary 


tension ¢@ is similarly given by 
o/kT = —log y=2(H—H"*). (122b) 


The boundary tensions both vanish at the 
critical point, as one should expect. The formulas 
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Fic. 6. Properties of quadratic crystal. — — — — — — 
Boundary tension ¢ between regions of opposite order. 


Specific heat C. ------------ Ap- 
proximate computation of C by Kramers and Wannier. 


for the boundary entropies and energies follow 
Se’ = —do’/dT =2k((H/sinh 2H)+H*), 
So= —do/dT =2k((H'/sinh 2H’)+H"*), 

up’ = o' + Ts,’ = (2/sinh 2H) J+2J’, 

us = o+7Ts,=2J+(2/sinh 2H’) J’. 


(123) 


The variation of o with the temperature for a 
quadratic crystal (J=J’) is represented graphi- 
cally in Fig. 6. 


Mean Ordered Length 


In a crystal of moderate width we may 
expect to find sections of rather well-defined 
order similar to that of an infinite crystal, with 
occasional transitions to opposite order. The 
result (122b) allows an estimate of the ‘‘mean 
ordered length’’ /) between such interruptions. 
By the general theory of fluctuations 


l=" = exp[2n(H—H’*) ] 


=(e" tanh H’)". (124) 


Evidently, when n is large enough we have 
l>>n. 
The estimate (124) is significant when this 
condition is fulfilled. 
The Specific Heat of a Finite Crystal 


Since the partition function of a finite crystal 
is an analytic function of the temperature, its 
specific heat will be finite at all temperatures. 
However, a maximum will occur near the 


transition poirit of the infinite crystal; this 
maximum will be increased and sharpened with 
increasing size of the crystal. 

Kramers and Wannier™ computed the specific 
heats at the transition point successively for 
n=1, 2, ---, 6 (screw arrangement). From the 
results they inferred the asymptotic rule that the 
height of the maximum increases linearly with 
log n. We shall verify this remarkable conjecture, 

The maximum of (C/N) for a finite crystal 
does not occur exactly at the asymptotic critica] 
point ; but the differences in location and magni- 
tude are only of the order n~ log n. This order 
of accuracy will suffice. 

To compute the specific heat we substitute the 
exact partition function (95) of a finite crystal 
in (111b), whereby the relations (112) are useful, 
The specialization to the critical case 


H'=H*; gd 2H+gd 2H’=}2; 
sinh 2H sinh 2H’ =1 


simplifies the computation in several respects: 
the relations 

6’ = 6* =5, 
are valid in this case. We write the result in the 
form 
C/Nk= —2H7? sinh? 2H’+ (H’+H sinh 2H’)? 


sinh }y=sinh 2H’ |sin }w}| 


1 n 
Xsinh 2H X— ¥ cosec[(r—4)x/n) 
N rm 


1 n 
+— > [(A’+H sinh 2H’)? 
N rend 


X (2 sin? 5e,_; coth y2,-1—cosech }72,_;) 

—4HH’ sinh 2H’ sin? 52,_, tanh }72,-1 
+ H? sinh 4H’ cos 52,-; ]. (125) 

Ls ] rn 


Lf I | 
Nk 
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Fic. 7. Specific heats for varying degrees of anisotropy. 
Filet, <n oo— J'/J=l 


(quadratic crystal); - ----- J’ =0 (linear chain). 
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The second summand, when considered as a 
periodic function of the variable we,_1 = (27 — 1) x/ 
an, has a differential quotient of bounded 
variation. Its sum may be replaced accordingly 
by the corresponding integral, with an error of 


the order n~*. 
We have purposely split off the first sum 


1 n 
Si=— L cosec((r—}3)x/n). 
N rant 


The corresponding integral diverges. For an 
asymptotic estimate of the sum we expand the 
cosecants in partial fractions; the result is easily 


arranged in the form 
1 1 1 
2n—1 2m+1 2n+3 





=(-+5+ + 
mma 3 


1 1 
on 
4n—1 4n+1 


2 > (—)™(v((m+1)n+ 4) —¥(mn+})) 


whereby 
¥(z) =I" (z)/T'(z) =log(s— 3) +O(z*). 


From the given asymptotic estimate of y(z) 
together with 


¥(}) = —2 log 2—Cz, 


Cg=0.577 215 665, (Euler’s constant) 


a +) = log _ (Wallis’ product) 


we obtain 
S,= (2/2) (log(8n/x) + Ce) +O(n-*). 


Substitution of the given estimates in (125) 
yields the asymptotic formula for the specific 
heat maximum 


C/Nk~(2/x)[sinh 2H(H’+4H sinh 2H’)? 
X (log n+log(8/x cosh 2H’) + Cp) 
+sinh 2H(H’ —H sinh 2H’)? 

—2(H' sinh 2H)? gd 2H’ 


—2(H sinh 2H’)? gd 2H]. (126) 


The result is not quite symmetrical with regard 


- > 


- + — 


> > 
- + - 


Fic. 8. Superlattice at low temperatures in crystal for 
which the interaction energy between neighbors of opposite 
spin is J>0 lengthwise and — J’ <0 transversely. 


to H and H’; neither was such symmetry to be 
expected. Interchange of the two leads to the 
replacement 


log(n/cosh 2H’)—log(n/cosh 2H). 
For the case of quadratic symmetry 
H=H' =} log cot r/8 


the coefficient of logm attains a maximum. 


We find 
C/Nk~ (2/2) (log cot 2/8)? 
X (log n+log(25/?/a) +Ce— in) 


= 0.4945 log n+0.1879. (127) 


The estimate indicated graphically by Kramers 
and Wannier'® may be described by the equation 


C/Nk=0.48 log n+0.21. 


APPENDIX 


We shall deal here with the evaluation of 
various integrals which occur in the text. Most 
of these can be reduced in straightforward 
fashion to complete elliptic integrals; only the 
partition function itself is of a type one step 
higher than the theta functions, and involves a 
little analysis which is not found in textbooks. 


1. General Notation 


One of the main reasons for setting this 
material apart is that the notation must be 
specified because there are several systems in 
common use. We shall adopt the notation used 
by Whitaker and Watson."® 


1” E. T. Whitaker and G. N. Watson, Modern Analysis 
(Cambridge University Press, 1927), fourth edition, 
Chapters 20-22. Referred to in the following as WW. 








Theta Functions 
33(z| r) = 0;(z) = d4(z+ 47) 


=> exp (n*xir) cos 2nz, 


—@ 


do(z| r) = Jo(z) = 31(2+ 437) 


= > exp[ ("+ 4)?xi7 ] cos(2n+1)z, 
sir (1.1) 
8,=08,(0); 3,”’=0,"(0); (r=2, 3, 4). 
Periods, Modulus, Comodulus 
2K=0;"; 2iK’=0;""7, 
k=sin 0=0,7/8;?, 
k’=cos 0=08,7/8;°, 
log g= wir = —2K'/K. 


(1.2) 


We shall only have to deal with ‘‘natural’”’ 
cases in which —ir is real (not negative); 
0<k=1; 0Sk’<1. The ratio 7 of the periods 
will be specified when necessary. Tables of log q 
are available.” 

Jacobian Elliptic Functions 
sn(#3°2) = (03/82)31(z)/d4(z), 
cn(03°z) = (04/82)02(z)/d4(z), 
dn(83"z) = (84/03) 03(z)/d4(z), 


am “u=sin—'(sn u) =cos~'(cn u). 


(1.3) 


We shall use Glaisher’s notation for the quotients 
of the Jacobian elliptic functions 


scu=snu/enu; ndu=1/dnu; etc. (1.4) 
The Functional Relations 
cn? u+sn? u=dn?u+k? sn?u=1 (1.5) 


are important for the uniformization of certain 
algebraic functions. 


Elliptic Integrals of the First Kind 


d(sn u)/du=cnudnu; etc. 
u= F(k,am w=f (1—k* sin? g)-'dy, (1.6) 
0 


K(k)=F(k, $x); K'(k)= K(k’). 


*¢ Four place tables of log g by 5’ intervals of the modular 
angle are given by E. Jahnke and F. Emde, Tables of 
Functions (B. G. Teubner, Leipzig and Berlin, 1933), 
second edition, p. 122. 
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Elliptic Integrals of the Second Kind 





These involve the following definitions ang 
relations. 


Z(u) =85-°04' (83-2u) /94(95-2u), 






x/2 
E=[ (1—k? sin? ¢)'dy 





(1.7) 





K 
-f dn? u du=(—82/'/d3od3")x/2, 
0 





E(u) =Z(u)+ (E/Kyu= { dn? v dv. 





Complete Elliptic Integral of the Third Kind 





K 
Kz(a)= [ k? sn acnadna sn? u 
0 





X(1—k*sn?asn*u)-'du. (1.8) 






2. Uniformizing Substitution 


The trigonometric functions of the angles and 
sides of a spherical triangle with two fixed parts 
can always be expressed in terms of single-valued 
functions by a suitable elliptic substitution. The 
same is true for hyperbolic (and plane) triangles, 
In any case the determination of the modulus js 
suggested by the sine proportion. 

For the triangle of Fig. 4, with two given sides 
2H’ and 2H*, the simplest uniformizing substi- 
tution depends on which side is longer. For 
temperatures below the critical point we choose 
the modulus k, the imaginary parameter ia, and 
the new variable u as follows 


k=ky=sin Oo=sinh 2H*/sinh 2H’ 
=1/sinh 2H sinh 277’ <1, (2.4 
am(ia) =am(2iKo'y, ko) = 2H’, 4 


am u=8’. 






















Then the trigonometric, viz., hyperbolic functions 
of 6*, 5’, 2H* and 2H’ are given by the formulas 





sin 6’=snu; sin 6*=ksn u 


cos 6’=cn u; cos 6*=dnu 


tan g’=sinh 2H’= —isn1 
an g’ =sin isn 1a (2.28) 





cot g=sinh 2H* = —1k sn ia 





sec g’=cosh 2H’ =cn ia 









cosec g=cosh 2H* =dn ia. 
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CRYSTAL STATISTICS 


We are using the abbreviation 

g’=gd 2H’; g=gd2H=}x—gd2H* (2.3) 
for the gudermannian angles, which can often 
be employed to advantage in numerical compu- 


tations. 
For temperatures above the critical point the 


starred and the primed quantities exchange roles: 
k=ky=sin Oo=sinh 2H’/sinh 2H* 
=sinh 2H’ sinh 2H= <1, 

(2.1b) 


am(ia) =am(2iKo'y, ko) = 21H", 


am u=6* 
and explicitly 
sin 6’=ksnu; sin 6*=snu 
cos 6’=dnu; cos é*=cnu 
tan g’=sinh 2H’ = —ik sn ia 
; oe (2.2b) 
cot g=sinh 2H* = —i sn ta 
sec g’=cosh 2H’ =dn ia 
cosec g=cosh 2H*=cn ia. 
The formulas for w and y are the same in the two 
cases. The abbreviation 


M=dn ia dn u—k cniacn u 
will be convenient ; then 
cosh y= (cn ia dn u—k dn ia cn u)/M 
sinh y= —ik” sn ta/M 
; ; (2.4) 
—cos w= (dn ia cn u—k cn ia dn u)/M 
sinw=k” sn u/M. 
Moreover, when w and y are considered as 
functions of u and a 
—dw/du=dy/da=k"/M, (2.5) 
0y7/duU=0w/da=tk"k sn ia sn u/M. , 


The conformal mapping indicated by (2.5) may 
be described by the functional relation 


cot $(w—ty) 
= (1+) sc }(u+ia) nd $(u+ia). 
3. Alternative Substitution 


(2.6) 


While the substitution given above is preferred 
for most purposes, there are certain advantages 
to be gained by a Landen transformation to 


modulus and periods given by 

ki’ =cos 9, =coth(H’+H"*) tanh| H’—H*| 

= |cos(g+g’) | /cos(g—g’) =(1—ko)/(1+ko), 
k,=sin 0,=2(sinh 2H sinh 2H’)*/ 
(1+sinh 2H sinh 2H’), (3.1) 

Ki=(1+ko)Ko, 
Ky’ _ 3(1 +ho) Ko’ = K,'/(1 +h’), 

71 =1K,'/K,=}19=}r. 
The parameter 


and the variable u; are given by 
am(ta,, ki) =i(H’+H"*), 
aaa ws peta | _ 
and the following relations hold 
am(K,—1a;) =}4—i|H’—H*|, 
am(K,—«) = }4r—}3|6’—6*|, 


am(u,+1a;) = }(4—w+i7). 


(3.4) 


Here the functions of the double parameter are 
of some interest 


sc(2ia,) = sc(4tyK,’) =i cos (g—g’), 


nc(2ia;) =nc(4iyK,’)=sin(g—g’), (3.5) 


dc(2ia;) = dc(4iyK,’) =sin(g+¢g’). 


This substitution shows up the analytic nature 
of the parameter through the critical point, and 
(2.6) is replaced by a simpler equivalent in (3.4). 
Moreover, formulas valid for all temperatures 
can be constructed by furnishing a sign for the 
comodulus 


k,'= |ky’"| . 


(3.6) 
ki" =coth(H’+H"*) tanh(H’ —H*). 


That advantage together with the previous usage 
of Kramers and Wannier"™ decided the choice of 
elliptic substitution in the text. 


4. Discussion of the Parameter 


The imaginary parameter ta given by (2.1) 
can be evaluated as a real elliptic integral of the 
first kind. For that purpose we use Jacobi’s 
imaginary transformation, described by the 
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relations 
ics(iu, k) sn(u, k’)=cn(iu, k) cn(u, k’) 
=dn(iu, k) cd(u,k’)=1 (4.1) 
and the parameter can be determined by the 
formula 
am(do, ko’) =am(2yKo’, ko’) =gd 2H’ =g’, 
ao= 2yKy' = F(cos Oo, g’). “_ 


In tables this elliptic integral is often denoted by 
F(90°— Qo, g’). Alternative formulas can be 
obtained by Landon transformations, whereby 
the “elliptic angle’”’ ry is invariant. From (3.4) 
and (3.5) we obtain, respectively, 


am(ai, cos @,;) =am(2yKy’, ki’) 
=gd(H’+H"*), (4.3) 

am(2a;, cos @;) =am(4yKy’, ky’) 
=}n—gt+g’. (4.4) 
Further Landen transformations are best de- 
scribed by the quotients of theta functions. The 


following are suitable for temperatures not too 
far removed from the critical point 


tan? 30.=sin 0, 
= (tan 2g tan 2g’)}/cos(g—g’), (4.5a) 
(cos 2)! sn(ae, cos O2) 
3:(2ey|—4/r) sinh| H’—H*| 
~ 94(2ey|—4/1) cosh(H’+H*) 
=|tan 3(g+g’—3x)|, (4.5b) 
3:(4xy| —4/r) 
34(4ry| — 4/7) 
=sin(g—g’)|cot(g+g’)|. (4.5c) 
For temperatures not too close to the critical 
point we may take 


k_,=sin O_,=tan? 36, 








(cos @2)! sn(2a2, cos O2) = 


(4.6a) 


and the theta series involved in the following 
formula (imaginary argument) will converge 
rapidly : 

8i(ery| 27) 

id 4(xry | 27) 

e?#: (H'>H*) 

= (4.6b) 

tanh Hl’; (H'<H*). 


eit sn(iyK_1’, k-1) = 


One or the other of the two formulas (4.5) and 
(4.6) will yield at least 3 significant figures per 
term of the theta series. When a table of log q is 
available (a table of complete elliptic integrals 
will do), it is not necessary to evaluate 42 or @_, 
explicitly. 

As regards the variation of the parameter q 
with the temperature, we note first of all that 


y =0/2Ko' =a;/2Ky' 


is analytic at the critical point. This is evident 
from (3.1) and (3.3); the elliptic integrals which 
determine a; and K;,’ are analytic functions. of 
cos? @;, which has no singularity there. In this 
sense the singularities of the various functions 
at the critical point are completely described by 
the degeneration of the period parallelogram 


t/r=K/K'=0; (H’=H"*). 
For extreme temperatures, with 
H=J/kT; H'=J'/kT 
the asymptotic estimates 


lim (2y) = H’/(H+A’) =J'/(J+J’), 

T=0 
lim (2y) =43 (4.7) 
T=02 


can be derived from (2.1). 

The parameter is an unsymmetrical function 
of H and H’; the nature of the asymmetry is 
evident from either of the formulas (4.4) and 
(4.5c), which involve the double parameter: 
Interchanging the roles of H and H’ replaces 
y by 4—y, and the parameter ia» is replaced by 
its elliptic complement i1(Ko’—ao). Thus (2.1a) 
and (2.1b) have the counterparts 


am(iK’ —ia, ko) =am(iKo'(1—2y), ko) 


=e (H’>H*) 4.4) 

2iH’*; (H'<H*). 

It is evident that if J’<J, we shall have 
y<4; a<}K’; (J’<J) (4.9) 


for all temperatures, and in the special case of 
quadratic symmetry 


y=}; a=}K’; (J=J'). (4.10) 


This is the reason why the results for that case 
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CRYSTAL STATISTICS 


can be expressed by the complete integrals K(k) 
and E(k) alone, as given in the text. 


5. Reduction of Integrals 


By the substitution (2.2), the integrals of (113) 
in the text are readily evaluated as complete 
elliptic integrals. For example, 


7” 2K 
f cos fdw= k’? cn u(dn ia dn u 
0 0 


—keniacnu)"du; (H’'>H"*). 


If we replace u by 2K—u, the sign of cnu is 
changed ; that of dn wu is preserved. Accordingly, 
that part of the integrand which is odd with 
regard to cn u may be omitted, which yields 


k cn ta cn? u k cn ia dn? ia sn? u 


————-=kcnta— 
1—k? sn* za sn* u 





1 —k?* sn* ia sn*u 

whence 

f cos 6’'dw = 2K(k cn ia—k" dsia Z(ia)); 
‘ (H’>H*). 


The other integrals (113) are treated similarly. 
The results are 


{ cos 6*dw = —2iK(k sn ia coth 2H* 

“ ~Z(ia) coth 2H’), 

y cos 6‘dw = — 2iK(k sn ia coth 2H’ 
—Z(ia) coth 2H*), 


f (sin 6* sin 6’/sinh y)dw 
0 
=2(K —E)(t/k sn ia), 
sinh? 2H { sin? 6* coth y dw 
0 


= sinh? 2H* f sin? 6’ coth y dw= — 2iKZ(ia). 
0 


The mixed notation is designed to throw the 
distinction between the cases (2.2a) and (2.2b) 
entirely into the choice of modulus and param- 
eter; so that the formulas (5.1) are equally 
valid for H’ <H* and for H’>H*. 

The Z function of imaginary argument which 
enters into (5.1) can be expressed in terms of the 
tabulated functions of real argument Z(a, k’) as 


follows: 
Z(iu, k) =dn(iu, k) sc(iu, k) 

—iZ(u, k’)—i(wu/2KK’). (5.2) 
For substitution in (113) it is convenient to make 
use of the function 
Z(iK’ —ia) = —dn iacsia—i(x/2K)—Z(ia) (5.3) 
together with Z(ia); this leads to more sym- 


metrical formulas. The connection between the 
corresponding functions of real argument is 
Z(a, k’')+Z(K’' —a, k’) 

=k"(—i) sciandia. (5.4) 
6. Derived Thermodynamic Functions 


In assembling the terms of the formula (113a) 
for the energy it is best to deal with the cases 
H'> H* and H' <H* separately. The results are 


U=—NJ'(@ log \/dH’) 
—NJ(d log \/dH), (6.1a) 
8 log \/0H’ =coth 2H’ (2y+(2K/x)Z(2yK’, k’)), 


(H’>H*) (6.1b) 
0 log A/dH =coth 2H 


X(1—2y+(2K/x)Z((1—2y)K’, k’)) 


valid for temperatures below the critical point; 
above the critical point we have instead 


8 log \/8H’ =coth 2H’ 
X (2y — (2K /w)Z((1—2y)K’, k’)) 


(H' <H*). (6.1c) 


0 log \/8H=coth 2H 
X(1—2y—(2K/x)Z(2yR’, k’)) 


The energy is continuous at the critical point 
because, when « is real [cf. (5.4) ] 


0SZ(u, k’) S1-—k. 


For the specific heat we obtain from (113b) 
and (5.1), in mixed notation valid for all temper- 
atures 


C 4 ; ; H' 2 
—=—{ —iKZ(éa)( ——_— 
Nk « sinh 2H’ 
: ; a |} 
—iKZ(R’—ia)( - 
sinh 2H 


+2(K B)( wean )uw’) (6.2) 
i sinh 2H’ a ai 
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Every term is positive. The specific heat is 
infinite at the critical point, where k= 1. This is 
caused by the factor K; the general case is 
substantially like the special case of quadratic 
symmetry treated in the text. 


7. The Partition Function 


To compute the partition function (106) in 
terms of 7 and y we first differentiate under the 
integral sign 


° ax OY Ow C*w 
ydo= f (-2=-, )au 
0 0 0a du  dadu 
Integration by parts yields with the aid of (2.5) 
Be OY Ow Pha. Ow 
f -—— du 
da du are da 
2K 
“S((a) +) 
-{ (= ia dn? u+k? cn? ia cn? “a 
— }du. 
1—k? sn? ia sn? u 


The integral is easily evaluated ; then considering 
that 








(0, u) =cosh—'!(1) =0 


for a=0, we have 


[ ydo=2K P(iu)du, (7.1) 
0 9 


&(u)=2 dn? u—k”?—2dnucsuZ(u). 


The function #(u) is even and periodic with the 
period 2K. Its singularities are poles at the points 


u=2mK+niK'; (m,n integers; n#0). 


Jacobi’s imaginary transformation (4.1), (5.2) 
yields 


(iu) = —k’?+dc(u, k’) ns(u, k’) 
X(2Z(u, k’)+(eu/KK’)). (7.2) 


Accordingly, the poles of ®(u) are simple; the 
residues are 


R(2mK +niK"') = R(niK')=(—)"nwi/K. (7.3) 


It is easily seen that #(u) satisfies the standard 


conditions for development in partial fractions 
gua function of cos(ru/K). The residues (7.3) 
together with the condition 


6(K)+(K+iK’) =k"—k?=0 


suffice to determine the expansion 


> 2 ) 
v= (Fe) (1-8 EO 


(g" cos(ru/K) —q*") ) , 
(1—2q" cos(xu/K)+q)/ 4) 
The integration required by (7.1) is readily 


performed, and the formula (106) of the text 
becomes 





rT 


F 
~peewlag ie $ log(2 sinh 2H) +— wg f x dw 


=} log(2 sinh 2H) —4y log g 
i— —q"ty 
-> (—)"n loe(- ae =). (7.5) 


n=l *y 





Specialization to the case y=} yields the formula 
(119a) of the text, where g:°=q. 

The series of (7.5) converges rapidly except in 
the immediate neighborhood of the critical point. 
To obtain an expansion suitable for computation 
in that region we note from (7.2) that 


&(iu) —(xu/2KK"')i cs iu dn iu 
=2dc(u, k’) ns(u, k’)Z(u, k”)—k" 


is a periodic function of u with the real period 
2K’. Its Fourier series is easily derived with the 
aid of the identity 


K(2 deans a Z(a) —k’) 
=f ¢en(u—2) sn(u+a) —cn(u—a) cn(u+a))du 


from those of the Jacobian elliptic functions." 
The result is 


ru 
(iu) ————i cs iu dn tu 
2KK’ 





=2(= yx 2 cos((2n+1)ru/K’') 
K’ <= sinh?((2n+1)"K/K’) 


2 aww.r reference 19, 22.6. 
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Moreover,” 
_ @ ' 
iu csiu dn in——_(s log(—ik! sn iu)) 
u 


= —log(k! sc(u, k’)) =log cot(ru/2K") 


© Qe-@nt)*K/K’ cos((2n+1)ru/K’) 


+2 (2n+1) sinh((2n+1)9rK/K’) 





Termwise integration of the Fourier series yields 
the alternative expansion for the partition 


function 


_ F/NkT =log d=} log(2 sinh 2H) 


— iy $(iu)du=y log(2 sinh 277’) 
© 


0 


vy 


+(}—y) log(2 sinh 2H)+2 } log cot (xz)dz 


1.2 (1+(4n+2)(xi/r) —exp[ — (4n+2)xi/r]) 
‘ (2n+1)? sinh?((2n+1)wi/r) 
Xsin((4n+2):y). 





T n=O 
(7.6) 
The formula (119b) of the text is obtained by 


specialization to the case H’=H; y=}, with 
7=27;, whereby 


7-24... 


x/4 
[ log cot x dx =1-*—3-° +5? — 
/9 


% WW, reference 19, 22.5, example 3. 
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8. Critical Data 


The transition temperature is given by the 
condition 


sinh 2H sinh 2H’ 
=sinh(2J/kT) sinh(2J’/kT) =1; 
Then 


(T'=T-). (8.1) 


r=0; K=ao; K'=}e 
(8.2) 
ry=a=gd 2H’ =}n—e¢d 2H. 
The formula (6.1) for the energy becomes 
U.= —(2/x)(NJ’' gd 2H’ coth 2H’ 
+NJ gd 2H coth 2H) 


and the formula (7.6) for the partition function 
simplifies to 


— F./ NkT =log \. =} log 2+(4—2y) log cot ry 


(8.3) 


v 
+2f log cotrzdz. (8.4) 
0 


The specific heat becomes infinite; the approxi- 
mate formula 


C/Nk~(4/x)(K(H+H’ sinh 2H)? sinh 2H’ 
— H? sinh? 2H’ gd 2H—2HH’ 


—H" sinh? 2H gd 2H’) (8.5) 


where 
K ~log(4/k’) ~}[log(47/ | T —T-|) 
—log(H coth 2H+H’ coth 2H") } 


is valid for temperatures near the critical point. 
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Continuous X-Ray Spectrum: 
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HE theoretical dependence of the continuous x-ray 

intensity from thin targets—relativity and retarda- 
tion of potential neglected—upon the angle of observation 
is given solely through the multiplicative factor 


Q(8) = (L,? sin? 8+ L,? cos? 8+ L,?), 


where @ is the angle between the direction of observation 
and that of the bombarding electrons; the quantities L,’, 
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L,’, L,? are calculable by means of complicated series ex. 
pansions.' The x axis is taken along the direction of the 
electron beam, while the y direction is perpendicular both 
to this and to the direction of observation. Through sym. 
metry, we have L,?=L/. 

Although the individual calculations of L,*, L,*(= L) 
can generally be carried through without an unreasonable 
amount of effort, the sum B=LP+Lf+LlLi=L2+213 
may be computed, usually, at the expense of appreciably 
less labor.* For this reason, comparison of experiment with 
theory may be facilitated somewhat if observation is made 
at an angle © for which the factor Q(@) is merely a multjple 
K of L, or 


L,* sin? 6+ L,?(1+ cos? @) = K(L,?+2L,?). 
Equating coefficients we obtain 
sin? @=K, 1+ cos? 0=2K, 
whence K =2/3; so that finally 
@ = 54°44’ or 125°16’. 

Hence it is suggested that—unless angular distribution js 
sought—future measurements upon the intensity of cop. 
tinuous x radiation from thin targets be made at an angle 
of 54°44’ with the (forward or reverse) direction of electron 


bombardment, provided the geometry of the experimental 
set-up is not rendered unwieldy as a result. 


1 R. Weinstock, Phys. Rev. 64, 276 (1943). 
2 A. Sommerfeld and A. W. Maue, Ann. d. Physik 23, 589 (1935), 
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PHYSICAL REVIEW 


HE 258th meeting of the American Physical 

Society was held at the California Institute 
of Technology, Pasadena, California, on Monday, 
December 27, 1943. Once again the meeting was 
very successful, in spite of the numerous handi- 
caps resulting from war conditions. The at- 
tendance varied from 50 to 70, and 40 members 
of the society and guests had luncheon together 
at the Athenaeum. Those in attendance came 
almost exclusively from Southern California, but 
limitations of travel made this fact almost 
inevitable. Much of the success of the meeting, 
and presumably the unexpected size of the 
attendance, was due to the four invited papers, 
which had been arranged by the combined efforts 
of Dr. Darrow, Dr. Millikan, and the local 
secretary. All of the four papers were important 
and were much appreciated by the audience. The 
number of contributed papers—fifteen—was far 
greater than had been anticipated, and resulted 
in a program almost too long for a single day. The 
morning session opened at 9:30 A.M. and ad- 
journed at 12:15 p.m. The afternoon session 


1. Analysis of a Multivibrator. S.C. SNowpon, California 
Institute of Technology (Introduced by W. H. Pickering).— 
A theoretical analysis of the action of a fundamental multi- 
vibrator has been undertaken along the lines of Van der 
Pol's analysis.! A more exact approximation to the transfer 
characteristic has been used and an exact solution of the 
equation has been found in the case where the shunt 
capacities of the tubes can be neglected. In this approxima- 
tion the period of oscillation has been calculated and 
checked experimentally. The limiting frequency of oscilla- 
tion in which the wave form is sinusoidal can be calculated 
as in Van der Pol's case. The results agree with experiment 
within three percent. 


1B. Van der Pol, Phil. Mag. 2, 978, 1926. 


2. A Voltage Regulator for High Voltages. W. H. 
PICKERING AND S. C. SNOWDEN, California Institute of 
Technology.—The standard degenerative voltage regulator 
circuit usually becomes difficult to use at voltages above a 
few kilovolts. This is particularly true in applications where 
the positive terminal is grounded. The new circuit avoids 
the difficulties by using a radiofrequency signal modulated 
by the variations in output voltage to transmit these 
variations across the potential difference to a tube in 
series with the high potential bus. The voltage limitations 
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opened at 1:30 p.m. and adjourned at 4:45 p.m. 
During the morning session two invited papers 
were presented, viz., 


Specificity of Intermolecular Interaction. Linus 
PAULING, California Institute of Technology. 

The Morphology of Nebulae and Clusters of 
Nebulae. Fritz Zwicky, California Institute of 
Technology. 


During the afternoon session the following two 
invited papers were presented: 


The Problem of Rendering Sea Water Potable. 
ALEXANDER GOETZ (in collaboration with S. J. 
Broadwell and D. C. Warner), California Insti- 
tute of Technology. 

Physical Methods in Hydrocarbon Research. 
Otto BEEcK, Shell Development Company, Emery- 
ville, California. 


Abstracts of the fifteen contributed papers 
follow. None of these was read by title. 


R. T. BirGe, 
Local Secretary for the Pacific Coast 


are accordingly removed and furthermore, either terminal 
of the output may be at ground potential. The sensitivity 
of the regulator depends as usual upon the total effective 
amplification of the circuit. Long-time stability depends 
on the stability of the amplitude or frequency of the radio- 
frequency signal. 


3. A Light Source for the Primary Standard of Wave- 
Length. Cuao-WancG Hsuen, California Institute of Tech- 
nology.—Despite some objections,! the red Cd _ line, 
6438.4696A, still remains the primary standard of wave- 
length for lack of a satisfactory substitute. In our labora- 
tory a light source of the Michelson type has been con- 
structed for this standard line. The interferometer spectro- 
grams show that the line is as narrow as, or even narrower 
than, that emitted from the Osira Cd lamp (Adam Hilger), 
and is very intense. Our lamp produces a satisfactory 
spectrogram with only thirty minutes’ exposure. Under the 
same conditions the Osira lamp gives a lighter spectrogram 
with six hours’ exposure. Furthermore, the intensity of our 
lamp can be varied over a wide range by changing the 
primary voltage of the transformer. The effective light- 
emitting area is 50X2.5 mm? which, if necessary, can be 
enlarged by means of a magnifying lens. The Doppler 
effect should disappear because the light is observed in a 
direction perpendicular to that of the current. Owing to 
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the high intensity, simultaneous observation of the stand- 
ard line and of lines from another source is possible either 
by using a half-reflecting mirror, or by placing the latter 
source in front of the spectrograph slit. These features 
eliminate the important objections against a light source 
for the standard Cd line. 


1W. F. Meggers, Rev. Mod. Phys. 14, 59 (1942). 


4. On the Elastic Vibrations of Cylindric Shells. PAuL 
S. Epstein, California Institute of Technology.—The prob- 
lem of the characteristic oscillations of curved shells with 
the usual border conditions (clamped edges and free edges) 
has remained hitherto unsolved. A method of treating 
elastic shells given by the author! permits of obtaining the 
solution, which was worked out for shells having the shape 
of circular cylinders. The result is particularly simple in 
the case of very thin shells. The three types of border condi- 
tions, namely, clamped edges, free edges, and supported 
edges, lead to the same characteristic frequencies. But the 
character of the oscillations, as evidenced by the distribu- 
tion of the nodal lines, is not the same for clamped as for 
free edges. 


1 Paul S. Epstein, J. Math. and Phys. 21, 198 (1942). 


5. Forbidden Lines of Call in a Stellar Spectrum. 
Pau W. MERRILL, Mount Wilson Observatory—A low 
dispersion spectrogram of the fourth-magnitude star v 
Sagittarii, type cApe, obtained by H. F. Weaver at Mount 
Wilson in 1942, showed the *S—?*P triplet of Ca II, \A8498, 
8542, 8662 in emission. Spectrograms with higher dispersion 
(21A/mm), obtained by the writer in 1943, showed in 
addition the forbidden doublet *S—*D, \7291.47, 47323.88. 
No previous observations of these forbidden lines either 
in the laboratory or in astronomical sources seem to be on 
record. The star is a spectroscopic binary of period 138 
days whose brighter component (the only one observed) 
has a range in radial velocity of 96 km/sec according to an 
orbit computed in 1914 by R. E. Wilson from the measured 
displacements of the dark lines. The bright lines of H and 
Ca II indicate little if any change in velocity. Hence they 
probably originate not in either component of the binary 
pair but in an enveloping gaseous cloud. 


6. Further Studies on the Least Penetrating Component 
of the Incoming Cosmic Rays. R. A. MILLIKAN, H. V. 
NEHER, AND W. H. PICKERING, California Institute of 
Technology—Although the war has prevented us from 
making such convincing tests on the existence of helium 
annihilation rays as we reported earlier! on silicon, oxygen, 
nitrogen, and carbon annihilation rays, we are able to 
present herewith heretofore-unpublished evidence obtained 
with the aid of high altitude electroscope measurements 
for the existence of approximately the energy to be com- 
puted from the helium annihilation hypothesis. 


! Millikan, Neher, and Pickering, Phys. Rev. 63, 234 (1943). 


7. The Average Specific Ionization of the Mesotrons 109 
Ft. Underground. Wayne E. Hazen, University of Cali. 
fornia.—Cloud-chamber photographs have been taken in a 
tunnel with 100 ft. of rock overhead. The chamber was 
operated in such a way that the ionization of cosmic-ray 
mesotrons could be obtained from the number of drops 
that condensed on positive ions. No magnetic field was 
used. A preliminary survey of the photographs has given 
a value of 60 ion pairs per cm in air at N.T.P. for ioniza- 
tion events in which less than 2500 ev of energy was lost 
by the cosmic ray per primary ionization (in other words; 
clusters of more than 80 drops were omitted). It was found 
that approximately 10 percent of the energy loss of a 
mesotron is attributable to collisions in which between 1200 
and 2500 ev is lost to an electron. There have been two 
determinations of the specific ionization of cosmic rays at 
sea level in this laboratory: one resulted in a value of 56, 
the other in a value of 46 ion pairs per cm of air at N.T.P. 
One expects an increase in the average specific ionization 
with depth since the average energy increases with depth, 
A rough calculation indicates that one might expect a 10 
percent increase in the average specific ionization at a 
depth of 100 ft. 


8. Friction of Binary Alloys with Thin Film Lubrication, 
S. Kyropou.os, California Institute of Technology.—An 
apparatus originally designed for investigation of lubricant 
characteristics' was converted to measure over-tempera- 
tures of alloys rubbing against steel. Equilibrium tempera- 
ture is usually reached in one-half hour; accuracy +0.3°C, 
As a basis of bearing alloys, first specimens were investi- 
gated covering the phase diagram Pb-Sb. Allowing for 
irregularities due to the crystallization process and solid 
solubility, the hardness of these alloys follows substantially 
a linear relationship. The overtemperature-composition 
curve is from 0-20 percent Sb almost its mirror image. 
From 20-100 percent Sb the curve goes through a minimum 
around 40 percent Sb, thus being substantially convex to 
the composition axis. This phenomenon may be inter- 
preted as an effect of the differential thermal expansion of 
the duplex structure, resulting in a shift toward more fluid 
lubrication. Apparatus and procedure followed permit 
study of precipitation hardening effects on friction and 
some insight into basic contributing factors such as plas- 
ticity and shearing off of particles and their reaction on 
various lubricants. 


1S. Kyropoulos, Refiner 18, 273, 320, 399 (1939); 19, 85 (1940). 


9. Present Status of Spectroscopic Methods of Analysis. 
M. F. Hascer, Applied Research Laboratories, Glendale, 
California.—Today’s versions of many instruments usually 
considered as tools of the physicist but now being used 
almost exclusively for the analysis of materials by the 
spectroscopic method, are described. The scope of the 
spectroscopic method is discussed as is its present day 
position. The future possibilities opened up by direct 
reading instruments, made practical by recent advances in 
vacuum-tube design, are considered. 
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10. Frequency-Band Character of Shock Vibrations. 
S, Leroy Brown, University of Texas.—A so-called shock 
vibration is one with a fundamental frequency but the 
amplitude is zero for more than half the cycle. Analysis of 
specific wave forms shows that, though the intensity of the 
harmonics decreases with frequency, there are bands of 
relatively greater intensity. For example a square-top 
wave of height A for a quarter of a cycle (height zero for 
three-quarters of the cycle) gives 1.4h/, 2h/2x, 1.4h/3x, 
0, 1.44/54, 2h/6m, 1.4h/7x, 0, 1.4h/9x, 2h/10x, 144/117, 
0, 144/132, 2h/14, etc., for the amplitudes of the har- 
monics. A vibration of 20 cycles per second, the wave form 
of which is half a sine curve for one-twentieth of the cycle 
(zero for nineteen-twentieths of the cycle), gives weak 
harmonics around 600 cycles per second, comparatively 
strong harmonics around 800 cycles per second, weak 
harmonics around 1000 cycles per second, etc. If the half- 
sine curve occupies a small fraction of a cycle, the ampli- 
tudes of the lower harmonics decrease very slightly with 


frequency. 


11. A High Sensitivity String Galvanometer. ARTHUR 
E. LocKENvITz, University of Texas.—A string galvanome- 
ter is described having a sensitivity as high as 10-* ampere 
per scale division and at this sensitivity a period of less than 
a tenth of a second. A new method of keeping the string 
under constant tension is used. A platinum wire 0.001 mm 
in diameter and about 7 cm long is fastened at the ends to 
rigid terminals but in such a way that the wire is very slack. 
Near one end of this platinum wire a second wire of the 
same diameter is tied so that it is at right angles to the first 
and its ends are fastened to fixed terminals like the first. 
This string assembly is in a magnetic field so that both 
wires are at right angles to the field. An auxiliary current 
of a few microamperes produces a force action in the second 
wire which causes the main string to be held under any 
desired constant tension. By means of a simple switching 
arrangement, the sensitivity and period can be varied over 
wide limits. The instrument is extremely stable and rigid. 


12. The Explosion of Ether-Air Mixtures by Electrical 
Sparks. A. M. ZareM, California Institute of Technology 
(Introduced by R. A. Millikan).—Sparks were passed be- 
tween spherical and pointed electrodes immersed in a 
combustible ether-air mixture. It was found that the energy 
in a spark required to explode the mixture decreased as the 
sparking potential was increased. At 4000 volts, a non- 
oscillatory spark current as low as 400 microamperes (peak 
value) exploded the mixture. At 2000 volts the peak current 
required for explosion was of the order of amperes. The 
shape of the electrodes affected the minimum energy neces- 
sary for explosion. No explosions were obtained using a 
prolonged corona discharge from a Wimshurst machine 
with voltages ranging up to 35,000 volts. Explosions were 
obtained, however, when a spark was passed between 
either pointed or spherical electrodes at 4000 volts. In 
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many cases the peak current in the spark which caused 
explosion at the lower voltage was not as great as the 
corona discharge current. 


13. Excitation of Electromagnetic Waves in Wave- 
Guides and Cavities. Joun Mites, California Institute of 
Technology (Introduced by R. A. Millikan).—The im- 
pedance of thin transverse wires in rectangular wave-guides 
and cavities has been evaluated by application of the vector 
potential. Both the resistive and the reactive components 
have been obtained, and one is thus enabled to place a 
driving or receiving element so as to match a given source 
or load for the maximum power transfer. The problem of 
short circuiting wires has also been calculated. The solu- 
tion of the analogous problems in a circular wave-guide is 
indicated. 


14. X-Ray Inspection of Spotwelds in Aluminum Alloys. 
R. C. McMaster, California Institute of Technology (In- 
troduced by R. A. Millikan).—This paper presents pre- 
liminary results of a study whose purpose is to develop 
practical industrial methods for the non-destructive testing 
of spotwelds in aluminum alloys. Procedures for x-ray 
inspection using industrial equipment have been developed, 
and correlations have been obtained between radiographs 
of spotwelds and spotweld structure, quality and strength. 
The scope and limitations of the method are indicated 
through the presentation of results of tests on thousands 
of industrially-made spotwelds. The results of preliminary 
tests employing fluoroscopic and special x-ray equipment 
are outlined. 


15. Sparking Potentials in Air in a Plane Parallel Gap 
at and below Atmospheric Pressure. Leon H. FIsHER, 
University of California.~—The streamer theory of the 
electric spark developed by Loeb and Meek" has indicated 
the importance of sparking potential studies in air at half 
an atmosphere pressure. Sparking potentials in a plane 
parallel gap in air have been measured as a function of plate 
distance and pressure. Below atmospheric pressure, the 
sparking potential curve as a function of distance at 
constant pressure shows a peculiar leveling off at large 
plate distances. This effect increases with decreasing 
pressure, and has been shown to be due to the effect of the 
metal walls of the large chamber used (50-cm diameter). 
Significant measurements at low pressures will require the 
use of larger chambers than have been used previously. 
The walls of the chamber do not influence the measure- 
ments at atmospheric pressure, and the sparking potential 
for a one-centimeter gap at atmospheric pressure at 22°C 
is 30.1 kv. The data at atmospheric pressure, when analyzed 
either by Townsend’s sparking equation or by Meek’s 
equation, are more self consistent than any other available 
data. Water vapor causes a small but definite scattering 
in the potentials at which sparks occur. 


1L. B. Loeb and J. M. Meek, The Mechanism of the Electric Spark 
(Stanford University Press, 1941). 





